
D. Herr mann G.Schneilhardt 

Schneider CPC 
Mathematik 



1 


Programm- 
Beispiele für den 
Anwender 


EE33 






Für Eva und Barbara 




D.Herrmann G.Schnellhardt 

Schneider CPC 
Mathematik 



Anwender 




ISBN 3-88 322-152-X 
1. Auflage 1985 

Alle Rechte, auch die der Übersetzung, Vorbehalten. Kein 
Teil des Werkes darf in irgendeiner Form (Druck, Fotokopie, 
Mikrofilm oder einem anderen Verfahren) ohne schriftliche 
Genehmigung des Verlages reproduziert oder unter Verwen¬ 
dung elektronischer Systeme verarbeitet, vervielfältigt oder 
verbreitet werden. 

Der Verlag übernimmt keine Gewähr für die Funktion einzel¬ 
ner Programme oder von Teilen derselben. Insbesondere 
übernimmt er keinerlei Haftung für eventuelle, aus dem 
Gebrauch resultierende, Folgeschäden. 

Schneider ist ein Warenzeichen der Schneider Rundfunk¬ 
werke GmbH & Co. 

Printed in Western Germany 
© Copyright 1985 by IWT-Verlag GmbH 
Vaterstetten bei München 

Holdenrieds Druck- und Verlags-GmbH, Füssen 
Umschlaggestaltung: Kaselow und Partner, München 



VORWORT 


Dieser Band enthält 39 mathematische Programme aus den Bereichen 

— Mehr-Register-Arithmetik 
— Zahlentheorie 
— Kombinatorik 
— Algebra 
— Geometrie 

— Numerische Mathematik 

Besonders nützlich werden Praktiker und Computerfans die hier neu vorgelegte 
Langzahl-Arithmetik finden, die die Grundrechenarten für Zahlen bis 255 
Stellen gestattet! 

Zahlreiche Anwendungen finden auch die hier angegebenen kombinatorischen 
Prozeduren: Wer wollte nicht schon mal alle Permutationen oder Kombina¬ 
tionen eines Problems durchspielen? Vielfältige Unterstützung findet auch das 
Rechnen mit Polynomen, Matrizen und komplexen Zahlen: Neben dem kom¬ 
plexen Hornerschema werden insbesondere Algorithmen zur Polynomdivision 
und Matrizeninversion gegeben. 

Anwendungsbezogen sind die Programme der Numerischen Mathematik. Hier 
werden Verfahren zur Lösung von linearen und nichtlinearen Gleichungen ge¬ 
geben, zum Eigenwertproblem von Matrizen und zur numerischen Integration 
und Differentiation. Alle Programme werden durch vollständige Beispiele und 
komplette Programmausdrucke erläutert. 


Dem Verlag danke ich für die Herausgabe des Bandes und für die stets freund¬ 
liche Zusammenarbeit. 


Anzing, im Juni 1984 
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ACHTUNG 


Bitte beachten Sie die folgenden Korrekturen! 


Seiten 21 und 24 und 29, hier bitte jeweils 
die Programmzeile 220 wie folgt ändern: 

220 IF M N THEN H$=A$:A$=B$:B$=H$:H=N:N=M:M=H 


Seite 82, hier ist das Ergebnis auf 6 Komma¬ 
stellen gerundet. 
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La machine arithmetique fait des effects 
qui approchent plus de la pensee 
que tout ce que font des animaux 

Pascal 


Daß die niedrigste aller Geistestätigkeiten 
die arithmetische ist, wird dadurch belegt, 
daß sie die einzige ist, welche auch durch 
eine Maschine ausgeführt werden kann 

Schopenhauer 


EINLEITUNG 


Mathematik und Computer 

Schon zu allen Zeiten hat es bemerkenswerte numerische Leistungen gegeben. 

Die berühmte Stieraufgabe, mit der Archimedesdie Mathematiker von Alexan¬ 
dria zur Verzweiflung trieb, hatte keine kleinere Lösung als 

10366482,7460514, 7358060 und 4149387. 

Ähnlich schwierig war die Lösung x = 1.3688081075 der Gleichung 
x3 + 2 x2 + lOx = 20 

die Leonardo von Pisa, Fibonacci genannt, um 1220 auf 10 Dezimalen genau 
berechnete. 

Van Ceulen (1539—1610) benötigte mehr als die Hälfte seines Lebens, um die 
Zahl pi auf 35 Dezimalen zu berechnen. Als Dank dafür hat man sie in seinen 
Grabstein in Leiden gemeißelt. Fermat (1601—1655) stellt in Briefen an seine 
armen Zeitgenossen mehrfach die Frage nach den Lösungen von 

x 2 i49 y 2 = i 

deren kleinste Werte x = 25801741449 und y = 2113761020 sind. 

In mühevoller Arbeit berechnete W. Shanks 1873 die ersten 707 Dezimalstel¬ 
len von pi. Ein Großteil dieser Mühe war jedoch vergeblich, da, wie sich später 
herausstellte, die 528ste Stelle bereits falsch war. 

An dieser Stelle muß auch der immense Rechenaufwand bedacht werden, den 
Ptolemaeus für die Erstellung seiner Sehnentafel Buergi und Brigg mit ihren 
Logarithmentafeln und Kepler mit den Rudolffinischen Tafeln hatte. 

Alle diese Probleme hätten in kurzer Zeit mit Hilfe eines Rechners gelöst wer¬ 
den können. Es erhebt sich die Frage, ob der Computer auch prinzipiell neue 
Fragestellungen in Angriff nehmen kann. 
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1976 kam es zu dem seltenen Ereignis, daß ein mathematischer Lehrsatz 
Schlagzeilen in der berühmten New York Times machte. Das sensationelle Er¬ 
eignis war, daß der Beweis des Vierfarbensatzes mit Hilfe von Computern ge¬ 
lungen war. 

K. Appel und W. Haken ernteten jedoch mit ihrem Beweis bei den meisten 
Mathematikern nur Kopfschütteln. Ein mathematischer Beweis war traditio¬ 
nell etwas, was man mit Papier und Bleistift nachvollziehen konnte. 

Wer aber sollte das 460-seitige Computerprogramm von Appel und Haken 
nachprüfen? Dazu kamen weitere 133 Seiten des lllionis Journal of Mathe- 
matics, auf denen die Autoren ihr Computerprogramm erläuterten. 

Das Vier-Farben-Problem, 1852 von Francis Guthrie aufgeworfen, besagt, 
daß jede ebene Landkarte mit 4 Farben so gefärbt werden kann, daß je 2 
benachbarte Länder verschiedenfarbig sind. Trotz erheblicher Anstrengungen 
widerstand die Vermutung über 120 Jahre lang den Beweisbemühungen der 
Mathematiker. 

Die Schwierigkeit war die, alle möglichen Lagemöglichkeiten der Länder in 
den Griff zu bekommen. Auf Grund der Vorarbeiten von Heesch (Hannover) 
gelang es Appel und Haken, alle möglichen Konfigurationen in einigen tau¬ 
send Fallunterscheidungen zu erfassen und vom Computer überprüfen zu 
lassen. 

Großen Aufschwung erfuhren die Teile der Mathematik, die sich mit endli¬ 
chen Mengen befassen, wie die Kombinatorik, Endliche Geometrie, Graphen- 
und Codierungstheorie, durch die Anwendung von Computern. Auch die 
Zahlentheorie erfuhr wesentliche Anregungen, so wurden z.B. sehr effektive 
Primzahltestverfahren entwickelt. Mit ihrer Hilfe konnte im September 1983 
die z.Z. größte bekannte Primzahl 
2 132 049 _i 

mit 39751 Stellen gefunden werden. Die Rechenzeit beträgt bei den hier be¬ 
nützen Lucas-Lehmer-Tests nur wenige Stunden. 

Konträr ist auch die Meinung der Mathematiker zu den stochastischen Prim¬ 
zahltests, die Primzahlen nur einer gewissen Wahrscheinlichkeit nachweisen. 
Beim Rabintest läßt sich die Irrtumswahrscheinlichkeit bei einigem Aufwand 
auf 

10 - 0 ° 

senken. Kann man dies als Beweis der Primzahleigenschaft ansehen? Oder ist 
nicht schon die Irrtumswahrscheinlichkeit eines Menschen größer? 

Bezeichnend ist die Geschichte von Rösser, Schönfeld und Yohe (Wisconsin), 
die bei ihrer Arbeit über die noch unbewiesene Riemannsche Hypothese einen 
Rechner benützten. Als "reine” Mathematiker mißtrauten sie dem Computer 
und überprüften in mühseliger Kleinarbeit dessen Betriebssystem. Tatsächlich 
fanden sie mehrere Fehler in der internen Rechnerlogik. Diese Fehler hatte 
bisher noch niemand bemerkt, obwohl der Computer bereits seit mehreren 
Jahren in Betrieb war. 
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Im Gegensatz zu den Primzahltests fehlt es bei der Primzahlfaktorisierung 
noch an effektiven Algorithmen. So weiß man seit langer Zeit, daß eine Zahl 
wie 

2256 + i 

zerlegbar ist, man kennt jedoch bis heute noch nicht ihre Primzahlzerlegung. 

Die Primfaktorisierung hat große Bedeutung in der Kryptologie gewonnen. 
Man kann nämlich damit Verschlüsselungsverfahren konstruieren und ihr Prin¬ 
zip offenlegen, ohne befürchten zu müssen, daß es "geknackt" wird (Verfah¬ 
ren von Rivest, Shamir und Adleman). Die Rechenzeit zur Faktorisierung 
einer lOOOstelligen Zahl beträgt z.Z. noch 4 Mrd. Jahre. Dies hat große poli¬ 
tische Bedeutung, falls es einmal zu einem generellen Atomwaffen-Sperrver- 
trag kommen sollte und die Codierung der Beobachtungssonden bekanntgege¬ 
ben werden müßte. 

Ein weiteres bekanntes bisher ungelöstes Problem der Zahlentheorie ist der 
sog. große Fermat-Satz. Er besagt, daß es keine ganzzahlige Lösung der Glei¬ 
chung =A 0 

x n + y n = z n 

für n>2 gibt. Gleichungen, von denen nur ganzzahlige Lösungen gesucht wer¬ 
den, nennt man Diophantisch. Vor kurzem (1983) konnte Faltings (Wupper¬ 
tal) zeigen, daß Diophantische Gleichungen vom Grad >3 höchstens endlich 
viele Lösungen haben können (Mordellsche Vermutung). Durch die Reduzie¬ 
rung auf endlich viele Fälle scheint das Problem nun mit dem Computer greif¬ 
bar. 

Im Falle der sog. Catalanschen Vermutung ist dies tatsächlich der Fall. Diese 
besagt, daß es keine zwei Potenzen >2 von natürlichen Zahlen — außer 8 und 
9 — gibt, deren Differenz 1 ist. Tijdemann konnte ebenfalls zeigen, daß das 
Problem höchstens endlich viele Lösungen hat. Durch umfangreiche Compu¬ 
terberechnungen glaubt er die Catalansche Vermutung mit an Sicherheit gren¬ 
zender Wahrscheinlichkeit bewiesen zu haben. Was ist von einem solchen Be¬ 
weis zu halten? 

Auch viele Fragestellungen der Kombinatorik, die infolge des großen Rechen¬ 
aufwandes bisher ungelöst waren, konnten seitdem beantwortet werden. Be¬ 
kannt ist das Eulerische Offiziersproblem (1782): Ist es möglich, 6 verschiede¬ 
ne Offiziere aus 6 verschiedenen Regimentern so in einem Quadrat aufzustel¬ 
len, daß in jeder Reihe und Spalte jeder Offiziersrang und jedes Regiment ver¬ 
treten ist? 

Eine solche Anordnung nennt man ein Paar von orthogonalen Lateinischen 
Quadraten. Euler vermutete, daß es keine orthogonalen Quadrate der Ord¬ 
nung 6, 10, 14, usw. gäbe. Erst in den Jahren 1959/60 konnten Bose und 
Shrikhande nach umfangreichen Rechnungen orthogonale Lateinische Qua¬ 
drate der Ordnung 10 und 14 angeben. Sie erhielten dafür den Spitznamen 
"the Euler Spoilers". 

Verwandt mit den orthogonalen Lateinischen Quadraten sind die Blockpläne 
der Endlichen Geometrie. Ausgangspunkt war das Schulmädchenproblem des 
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Geistlichen Kirkman (1851): Kann man 15 Schulmädchen an den 7 Wochen¬ 
tagen so in Dreierreihen spazieren führen, daß jedes Mädchenpaar an genau 
einem Tag zusammentrifft? 

Spezielle Blockpläne sind die sog. Steinersysteme S(t,k,v). Dabei müssen 
v Dinge zu je t Stück in k Reihen angeordnet werden können. Die Lösung des 
Schulmädchen-Problems stellt somit ein (2, 3, 15)-Steinersystem dar. Es ist 
im Rahmen des Buches nicht möglich, die Vielzahl der ungelösten Probleme 
im Zusammenhang mit den Steinersystemen aufzuzählen. Erwähnt sei, daß 
die kleinsten Steinersysteme, deren Existenz bisher unbewiesen sind, folgende 
sind: 

S (4,5, 17), S (4, 5,21), S (4, 5, 27) 

Durch umfangreiche Rechnung konnten Mendelsohn und Hung vor wenigen 
Jahren z.B. zeigen, daß das Steinersystem S (4, 5, 15) nicht existiert. Hier 
bleibt noch viel Forschungstätigkeit für künftige Mathematiker und viel 
Rechenarbeit für die nächste, noch schnellere Computergeration zu tun. 



"Aus Margarita Philosophica des Gregor Reisch, Freiburg 1503. 
Boethiusmit indischen Zahlen und Pythagoras als Linienrechner 
vorder Arithmetica sitzend" 
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Es macht mich traurig, daß gebildete Leute 
nicht einmal von der Existenz 
meines Fachgebietes wissen 

P. R. Haimos 


Die Furcht vor der Mathematik .. 
steht der Angst erheblich näher 
als die Ehrfurcht 

F. Auerbach 


Mathematik — die ungeliebte Wissenschaft 

Kein Zweifel, die Mathematik war schon immer eine unpopuläre Wissenschaft. 

Besonders dem Christentum waren die heidnischen Schriften eines Euklid- 
Diophant und Archimedes als heidnisch suspekt. Tertullian schreibt: 

"Nach Christus brauchen wir keinerlei Wissbegier mehr, nach den 
Evangelien sind keinerlei Forschungen mehr nötig". 

Als Folge der Verketzerung der Heiden, wurde die berühmte Mathematikerin 
Hypathia, Tochter des Akademievorstehers Theon, 415 von einer fanatischen 
Menge zu Tode gemartert. Dies, obwohl der Kirchenlehrer Synesios von 
Kyrene, Bischof von Ptolemais, ein Hörer ihrer Vorlesung war. Sogar die 
Mathematik hatte den Beigeschmack des Heidnischen. So lautete eines der 
Gesetze des Justinian-Kodex: "de Maleficis, mathematicis et ceteris similis" 
(Über Übeltätern, Mathematikern und dgl.): 

"Vollständig verboten ist die verdammenswerte Kunst der Mathe¬ 
matik" 

Die Mathematiker werden dort mit Giftmischern, Zauberern und Sterndeutern 
(Veneficis, Magis, Chaldeis) in einen Topf geworfen. Noch 1614 bezeichnete 
der Pater Caccini in einer Predigt, die gegen Galilei gerichtet war, die Mathe¬ 
matik als Teufelskunst und die Mathematiker als Urheber aller Ketzereien, die 
man aus allen Staaten vertreiben müsse. 

Er stützte sich dabei auf den heiligen Augustinus, der in De genesi ad literam 
geschrieben hatte: 

"Es besteht die Gefahr, daß die Mathematiker mit dem Teufel im 
Bunde, den Geist trüben und in die Bande der Hölle verstricken". 

Als 529 Kaiser Justinian die Athener Philosophenschule schloß, erlosch eine 
lange Tradition mathematischer Lehre. Noch bis 485 hatte an dieser Schule 
Proklos gelehrt, dem wir einen grundlegenden Euklid-Kommentar verdanken. 
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Es ist ein ganz wesentlicher Verdienst der arabischen Kultur ab 800 die wich¬ 
tigsten Werke des Apollonius, Archimedes, Heron, Euklid, Diophant, Ptole- 
maeus usw. in ihre Sprache übersetzt und somit der Nachwelt erhalten zu 
haben. Einige Übersetzungen aus dem Griechischen ins Lateinische führte 
auch Boethius (480 - 524) durch. Von ihm stammt auch die im Mittelalter 
übliche Einteilung der Künste in das Trivium (Grammatik, Rhetorik, Logik) 
und das Quadrivium (Arithmetik, Geometrie, Astronomie und Musik). 

Um eine mathematische Allgemeinbildung in Gang zu bringen, führte der 
Schotte Alkuin von York, im Auftrag Karls des Großen, 789 eine Schul¬ 
reform durch. Sie forderte, daß neben Psalmen, Noten und Grammatik ins¬ 
besondere auch der Computus gelehrt werden solle. Unter dem Computus 
Ecclesiaticus verstand man die Berechnung der beweglichen Kirchenfeiertage. 
Es war nämlich der Kirche ein Dorn im Auge, daß wegen der mangelhaften 
Rechenfertigkeit der Geistlichen das Osterfest zu verschiedenen Zeitpunkten 
gefeiert wurde. Alkuin selbst schrieb ein Rechenbuch zu Unterrichtszwecken. 


Der Rechenunterricht scheint — wie heute — nicht die besondere Freude der 
Schüler gewesen zu sein. Ein Mönch namens Strabo der Klosterschule Reiche¬ 
nau schreibt 822 in sein Tagebuch: 

"Zur Abwechslung und Unterhaltung lösten wir mathematische Rät¬ 
sel, welche Alkuin für den großen Karl gefertigt hatte. Viele vermoch¬ 
ten nicht allen diesen Rechnungen zu folgen, und bevor wir zur Geo¬ 
metrie übergingen, traten diejenigen aus, welche sich fortan den 
Studien der Medizin, Rechtswissenschaft und den Künsten der Male¬ 
rei und der Bildhauerei widmen wollten". 

Ein Mathematiker, Gerbert von Aurillac (940 - 1003) bestieg als Sylvester II 
999 den Stuhl Petri. Gerbert hatte einige Jugendjahre in Spanien verbracht 
und dort von den Arabern den Abakus (Rechenbrett) kennengelernt. Als Leh¬ 
rer der Domschule Reims verfaßte er ein Lehrbuch zum Rechnen am Abakus. 
Auch an den aus den Domschulen hervorgegangenen Universitäten Bologna 
(1088), Paris (1150), Salerno (1173) und Montpellier (1180) spielte die 
Mathematik keinerlei Rolle. 

Von der Universität Erfurt weiß man aus Quellen, daß die dort wirkenden 
Humanisten die Anstellung eines Mathematikers strikt ablehnten. Noch An¬ 
fang des 16. Jahrhunderts legten an der Universität Paris die Kandidaten für 
den Grad des Magisters der Künste kein Examen im Fach Geometrie ab. Sie 
mußten lediglich beeiden (!), Vorlesungen über die ersten 6 Bücher des Euklid 
gehört zu haben. 

Das Risiko einer wirklichen Prüfung wollte man nicht eingehen. Bezeichnend 
ist, daß der Lehrsatz des Pythagoras am Ende des 1. Kapitels den Namen 
"Magister Matheseus" trug und somit als Weisheit letzter Schluß galt. 

Welche Schwierigkeiten das Rechnen damals machte, wird aus der Rede 
Philipp Melanchtons vor der Universität Wittenberg (1517) deutlich. 
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. . Deshalb können ihre Anfangsgründe (der Rechenkünste) gar 
nicht dunkel und schwer sein; sie sind im Gegenteil so durchsichtig, 
daß Kinder sie begreifen können, weil ja alles so natürlich vor sich 
geht. Die Regeln des Vielfachen und Teilens allerdings erfordern viel 
mehr Fleiß, aber ihr Sinn wird schon bald von den aufmerksameren 
eingesehen werden. Übung und Anwendung erfordert diese Fertig¬ 
keit wie alle anderen." 



"Ein Vater meldet seinen Sohn bei einem Rechenmeister an. 
Holzschnitt von Hans Weiditz, 1535" 


Die Mathematik-Kenntnisse dieser Studenten aus den Lateinschulen sind ent¬ 
sprechend schlecht gewesen. Der Rechenunterricht wurde in der Regel von 
einem externen Rechenmeister übernommen, die dazu meist eigene Rechen¬ 
bücher verfaßten. Der bekannteste Rechenmeister war Adam Ries(e) in Erfurt 
und Annaberg. 

Es dauerte noch sehr lange, bis Mathematik ein reguläres Unterrichtsfach war. 
So wurde noch 1805 (!) im Jahresbericht eines lutherischen Gymnasiums in 
Essen vermerkt, daß an der Schule kein Rechenunterricht erteilt werde. Dies 
sei der Fall, weil ältere Schüler den Lehrern bekannt gemacht hätten, "daß in 
der Stadt jetzt ein fertiger Rechenmeister sei, der eine kürzere Methode habe". 

Oft wurde der Mathematik-Unterricht von nicht ausgebildeten Lehrkräften, 
z.B. von Theologen übernommen, die in Wartestellung auf eigene Pfarrpfründe 
waren. Ob es wohl diese Theologen waren, über die es in der Prüfungsordnung 
von 1820 für Preussen heißt: "Um das Eindringen untüchtiger Objekte in den 
Höheren Schulen Einhaltzu gebieten . . .". Auch in Bayern waren viele Mathe¬ 
matik-Lehrstühle im 19. Jahrhundert noch mit Theologen besetzt. 
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Daneben litt die mathematische Lehre noch unter anderen Schwierigkeiten. 
Mollweide in Leipzig erklärte es für unmöglich, neben der, für alle Fakultäten 
verbindlichen Mathesis pura, auch höhere Mathematik zu bringen, weil es da¬ 
bei zuviel Schreibens an der Tafel gäbe. 

Bezeichnend ist auch der Bericht von Neumann in Berlin von 1818: "Als ich 
mich beim Professor für Mathematik meldete, sagte dieser, ja ich habe die 
Vorlesung angezeigt, aber sie pflegt niemals zustande zu kommen. Ich verab- 





"Darstellung der Arithmetica 
Holzschnitte um 1500" 



Geometria aus dem Quadrivium, 
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redete mich mit 5 anderen, zu ihm zu gehen. Der Professor kam ins Audito¬ 
rium und schrieb ununterbrochen Formeln an die Tafel und sprach dabei kein 
Wort, bis die Zeit um war. Am 2. Tag kamen nur noch 2 Zuhörer. Der Profes¬ 
sor stellt sich wieder an die Tafel und zeichnete wieder ununterbrochen 
mathematische Formeln an diese. Am dritten Tag kam außer mir nur noch 
ein Zuhörer. Der Professor erschien, ging ans Katheder, wandte sich an uns 
und sagte: Sie sehen, meine Herren, es kommt kein Kolleg zustande". 



'Darstellung der Musica und Astronomia, Holzschnitte um 1 500" 
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1. Addition langer Zahlen 


Jeder, der schon länger mit dem Computer rechnet, ist schon einmal an die 
Grenze der Rechengenauigkeit gestoßen. Während bei Großrechenanlagen 
meist eine Langzahl-Arithmetik hardware-mäßig implementiert ist, ist dies bei 
Mikrocomputern nicht der Fall. 

Daher sollen im folgenden, Programme zu den Grundrechenarten für lange 
Zahlen angegeben werden. 

Diese Programme setzen voraus, daß die eingegebenen Zahlen kommafrei und 
positiv sind. Dies ist keine große Einschränkung, da die Vorzeichenregeln 
leicht von Hand ausgeführt werden können. Kommafreie Zahlen können 
durch entsprechendes Ausklammern von Zehnerpotenzen erhalten werden, 
z.B. 

12.3456789 = 123456789 * 10' 7 
0.000987654 = 987654 * 10' 9 

Alle Zahlen werden als Stringvariablen eingegeben und können daher bis zu 
255 Stellen haben, sie werden innerhalb des Programms in Felder zerlegt. Da¬ 
bei kann jedes Feld 1 bis 8 Ziffern umfassen. 

Umfaßt ein Feld mehrere Ziffern, so ist zu beachten, daß der Rechner führen¬ 
de Nullen nicht ausdruckt. Mit Hilfe der Stringfunktion kann man auch füh¬ 
rende Nullen ausdrucken lassen, wie es z.B. im Programm 7 geschieht. 

Wird der 1. Summand durch das Feld U(I), der zweite durch V( I) dargestellt, 
so kannbei ziffernweiser Addition das Addieren wie folgt beschrieben werden: 

90 K = 0 

100 FOR l=N TO 1 STEP-1 
110 W= U(l) + V(l) + K 
120 W(l) =W - INT(W/10) * 10 
130 K = INT(W/10) 

140 NEXT I 
150 W(0) = K 

Dabei ist N die Länge des größten Summanden, K der Übertrag. Die Summe 
ist im Feld W( I) gespeichert. Bei sehr großen Rechnungen, wie z.B. bei der 
Berechnung von 1000 Dezimalzahlen, zerlegt man die Summanden z.B. in 
Fünferblöcke U(l) bzw. V(l). 
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Die Fünferblock-Addition läuft wie folgt ab: 

90 K = 0 

100 FOR l=N TO 1 STEP -1 
100 W = U(l) + V(l) + K 
110 K = INT (W/10t 5) 

120 W(l) = W - K ♦ 10t5 
130 NEXT I 
140 W(0) - K 

Dabei ist N natürlich die Anzahl der Fünferblöcke. 


Zum folgenden Programm 

Als Beispiel werden im Programm die Zahlen: 
7830412596412348 

und 

914035812769267 
addiert. Hier ergibt sich die Summe 
8744448439181615 


Das hier gezeigte Programm arbeitet mit Einerblöcken, so daß auch Nullen 
richtig ausgegeben werden. 
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100 Rb.M Addition langer Zahlen 
110 : 

120 CLS:PRINT STRINOT(30,42) 

130 PRINT " Addition langer Zahlen" 

140 PRINT STRINGT(30,42) 

150 : 

160 DEFST'R A,B 

170 PRINT:INPUT "Eingabe 1.Summand";A4 
180 1'NPUT "Eingabe 2. Summand" ; B4 
190 : 

200 'Austauach und Auf f ui 1en der Stellen 

210 N-L.EN (A4) : M=LEN (B4) 

220 IF M>N TUEN SWAP A4,B4:SWAP N,M 

230 IF N~M TFIEN 280 

240 FOR 1-1 TO N-M 

250 B4="Ü"+B4 

260 NEXT I 

270 : 

280 ' Umwand1ung in Zah1en 
290 DEFINT I~k’,U~W 
300 DIN U (N) , V IN) ,W (N) 

310 FOR 1=1 TO N 

320 U(1)=VAL(MIOS(A4,1,1)) 

330 V(I)=VAL(NID4(B4,1,1)) 

340 NEXT I 
350 : 

360 ' st. e 11 enwei se Ad d i 1 1 on 

370 J=N:K=0 

380 WHILE J>=1 

390 W=U(J)+V(J)+K 

400 K=INT(W/10):W(J)=W NOD 10 

410 J=J-1 

420 WEND 

430 W(0)=K 

440 : 

450 PR INT:PR INT"Summe ="; 

460 K=l-K 

470 FOR I=K TO N 

480 PRINT W (I) ; 

490 NEXT I 
500 END 
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2. Subtraktion langer Zahlen 


Die Subtraktion langer Zahlen läuft als inverse Rechenoperation ähnlich wie 
die Addition ab. Da, wie anfangs erwähnt, auf das Rechnen mit Vorzeichen 
verzichtet wird, muß vorausgesetzt werden, daß der Minuend großer ist als der 
Subtrahend. Ist dies nicht der Fall, so vertauscht man Minuend mit Subtra¬ 
hend und versieht die Differenz mit einem negativen Vorzeichen. 

Wird der Minuend durch das Feld U(l), der Subtrahend durch V(l) dargestellt, 
so kann bei ziffernweiser Subtraktion das Subtrahieren wie folgt beschrieben 
werden: 


90 

K = 0 


100 

FOR l=N TO 1 STEP -1 

110 

W= U(l) - 

VO) + K 

120 

K = INT (W/10) 

130 

W(l) = W - 

K * 10 

140 

NEXT 1 



Dabei ist N die Länge des Minuenden, K der Übertrag. Die Differenz ist im 
Feld W(l) gespeichert. Übersteigt der Subtrahend den Minuend, so wird eine 
entsprechende Fehlermeldung ausgegeben. 

Bei sehr großen Rechnungen, wie z.B. bei der Berechnung von Fünferblöcke 
U(l) bzw. V(l). 

Die Fünferblock-Subtraktion läuft wie folgt ab: 


90 

K = 0 


100 

FOR l=N 

TO 1 STEP -1 

100 

W = U(l) 

- V(l) -K 

110 

IF W> = 

0THEN K = 0: GOTO 130 

120 

K = 1: W 

= W + 10t5 

130 

WO) = W 


140 

NEXT 1 



Dabei ist N natürlich die Anzahl der Fünferblöcke. 


Zum folgenden Programm 
Als Beispiel werden im Programm die Zahlen: 
784103516910456 

und 

9541217036428 

subtrahiert. Es ergibt sich die Differenz 
774562299874028 
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10 0 R E H B u b t. r a k t i o m 1 a n qer Zahle n 
1 10 : 

120 CLS:PRINT STR INGf(30,42) 

130 PR I NT " Bub t r a k 1 .1 on langer Zahlen" 

140 PR I NT STR INGf(30,42) 

150 : 

160 DEFST’R A,B 

170 PRINT:INPUT "Eingabe Minuend"; Af 
1BO INF'UT "Eingabe Subtrahend";Bf 
190 ; 

200 " A u s t a u s c h u n d A u f -f ü 1 len d e r S t e 11 e n 
210 N~LEN < Af) : M--LEN (B f) 

220 IF M: N TUEN SWAP Af , Bf:SWAP N , M 

230 IF N=M THEN 2S0 

740 FOR 1=1 TO N-M 

250 Bf= " 0 " +Bf 

260 NEXT I 

"v p, 

280 ' U mwand1unq in Zah1en 
290 DEFINT I~-K,U-W 
300 DIN LJ (N) , V (N) , W (N) 

310 FOR 1=1 TO N 

320 U <I)=VAL(MI DT(Af,1,1)) 

330 V ( T ) --VAL < MI Df (Bf , I , .1 ) ) 

'•40 NEXT I 
350 : 

360 'ste11enweise Subtraktion 

370 J=N:K=0 

380 WHILE J>0 

390 W=U( J) -V( J )+ K 

4 C > 0 K = IN T (W / .10) : W (J ) = W -- K * 10 

410 J=J-1 

420 WEND 

430 W(0)=K 

440 IF K=— 1 TUEN PRINT"Minuend 
< Subtrahend":END 

450 : 

460 PRINT: PRINT "Di-f-ferenz ="; 

470 K=0 

480 IF W (K) = O T H E N K=K +1 ; G O T0 480 
490 FOR I--K TO N 
500 PRINT W <1)5 
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3. Multiplikation langer Zahlen 


Nicht ganz so einfach wie die Addition und Subtraktion ist die Multiplikation 
langer Zahlen. 

Wird der größere Faktor stellenweise in das Feld U(I) und der zweite in V(I) 
zerlegt, so läßt sich die stellenweise Multiplikation wie folgt programmieren: 

100 FOR J = M TO 1 STEP —1 

110 IF V(J) = 0THEN W(J) = 0: GOTO 190 

120 K = 0 

130 FOR l=N TO 1 STEP-1 
140 T = U(I) * V(J) + W(l+J) + K 
150 W(l+J) =T - INTOV10) * 10 
160 K = INT(T/10) 

170 NEXT I 
180 W(J) = K 
190 NEXT J 

Dabei ist M die Stellenzahl des größeren Faktors, N die des kleineren Faktors. 
Im Feld W(l) wird die Ziffernfolge des Produkts gespeichert. 

Ähnlich, wie bei der Addition und Subtraktion, könnte man auch die Multi¬ 
plikation in Dreier-bzw. Viererblöcken durchführen. 


Zum folgenden Programm 

Das Programm berechnet nach dem oben angegebenen Verfahren das Produkt 
zweier langer Zahlen. 

Als 1. Beispiel wird die Zahl 
999999999999999 
quadriert. Es ergibt sich 

999999999999999800000000000001 

Im 2. Beispiel werden die beiden Primfaktoren 
59649589127497217 

und 

5704689200685129054721 
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der Zahl 


2128 + 1 

miteinander multipliziert. Diese Primzahlfaktorisierung gelang erst 1970 den 
Amerikanern Morrison und Brillhart, obwohl man schon sehr lange wußte, 
daß diese Zahl keine Primzahl ist. 

Das Programm liefert 

340 82366 92093 84634 63374 60743 17682 11457 
(vgl. Programmausdruck). 
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100 REM Multiplikation langer Zahlen 
110 : 

120 CLS:PRINT STRING*(35,42) 

130 PRINT " Multiplikation langer Zahlen" 
140 PRINT STRING*(35,42) 

150 : 

ISO DEFSTR A,B 

170 PRINT:INPUT "Eingabe 1.Faktor";A* 

180 IMPIJT "Eingabe 2. Faktor" ; B* 

190 : 

200 'Austausch und Auf-füllen der Stellen 
210 N--L.EN (A*) : M=LEN (B*) 

220 IF M >N TUEN SWAP A*,B*:SWAP N,M 

230 : 

240 Umwandlung in Zahlen 

250 DEFINT I-K,U-W 

260 DIM U(N),V(M),W(N+M) 

270 FÜR 1=1 T0 N 

280 U (I) = VAL. < MID*(A*,I , 1 > ) 

290 NEXT I 

300 FÜR 1=1 Tü M 

310 V(I)=VAL(MID*(B*,1,1)) 

320 NEXT I 

330 FÜR I=M+1 TD N+M 
340 W(I)=0 
350 NEXT I 
360 : 

370 'stellenweise Multiplikation 
380 J=M 

390 WHILE J >=0 

400 IF V(J)=0 THEN W(J)=0:G0T0 470 

410 I=N:K=0 

420 WHILE I>=0 

430 T=U(I)*V(J)+W(I+J)+K 

440 K=INT(T/10):W(I+J)=T MOD 10 

450 1=1-1 

460 WEND 

470 J=J-1 

480 WEND 

490 :: 

500 PRINT:PRINT"Produkt ="; 


29 






o 


(N 

Kl 




0 0 


Kl 

MD 




o 


■a 




o 


K) 




o 


MD 




o 


«t 




CO 


co 




cn 


KJ 




0' 


0' IN 




cn 


O ÜD 




cn 


CM <t 





H 






CM 





cn 

N 

cn th 





<a 






in 



* * 


cn 

O IN 

MD -i 


* * 



cn 


o 

* * 

Cn 0* 


C M CM 


CM 

* * 

CN Cn 

cn 

t-i 

MD CM 

Ul 

* * 

U- tN 


in cn 



* c * 

0- 0- 


CD <t 


o 

* Ql * 

£> Cn 

cn 

MD N 

Ki OD 

1— 

*-« * 

0' Cn 


O CM 


o 

* JC * 

Cn Cn 


* w * iH 


0 

* Ti * 

O O 

cn 

cm o- 

CM MD 

■ 1 

* NI * 

0- O 


cn 03 


tH 

* * 

0- Cn 


co in 


+ 

* L * 

Cn Cn 

cn 

MD Cn 

OD In 

DZ 

* Ql * 

Cn Cn 


'f 


11 

* cn * 

Cn Cn 


O MD 


DZ 

* c * 

0' Cn 

cn -i 

in cn 

CM -i 


* m * 

cn cn 


in in 


z r 

* -i * 





LU + ... 

* * 

f- c- 

cn o 

c- c- 

O 10 

X Z ~ 

* c * 

L. L 


L L 


H M 

* o * 

Q 0 


0 0 


□ ~ 

* -H * 

-M +J 

Cn O 

+J +J 

<a «t 

O 1- 3 

* -P * 

XL XL 


XL XL 


1! 

* <S * 

nj Qi 


rrj ro 


~ DZ H 

* Ui * 

Ll Ll 

9 

0 

Ix Ll 

Kl In 

1 II Z H 

* -rf * 

■ ■ 


■ • 


— 1—1 >-H 

*rt * 

-< C M 

ll 

*-> CM 

II 

3 ir H 

* Q. * 


o 



KixD 

* 'H * 

Ql Ql 

+J 

01 Oi 

XI 

Ix O ÜJ Z 

* X< * 

XI X! 


XI XD 

XL 

H lL Z LU 

* .-I * 

(5 Üi 

Zj O 

ns oi 

3 MD 


* 3 * 

cn cn 

TJ 

cn cn 

TJ 

o o o o o 

* Z * 

C C 

0 

c c 

0 

cm ro «t m -ß 

* * 

■H -r-t 

L O 

•i-t r-l 

l <a 

in in in ui in 

* * 

LU LU 

LL 

LU UJ 

X 


30 



4. Division langer Zahlen 


Wesentlich komplizierter als Addition und Multiplikation ist die Division lan¬ 
ger Zahlen. 

Das Programm führt eine stellenweise Division durch, entsprechend dem 
schriftlichen Rechnen. Zuerst wird geprüft, wie oft der Divisor in den Dividen¬ 
den hineingeht. Sodann wird der Divisor mit dem entsprechenden Vielfachen 
multipliziert. Diese Multiplikation einer langen Zahl mit einem einstelligen 
Faktor wird bei Programm 5 erklärt. Das sich ergebende Produkt muß noch 
vom Dividenden subtrahiert werden. Der genaue Algorithmus kann dem fol¬ 
genden Programm entnommen werden. 


Zum folgenden Programm 

Als Beispiel wird die schon aus dem Programm 3 bekannte Zahl 
2128 + i 

durch den Primfaktor 

59649589127497217 

dividiert. Das Programm liefert den ganzzahligen Quotienten 
5704689200685129054721 
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100 REM Division langer Zahlen 
110 : 

120 CLS:PRINT STRING*<30,42) 

130 PRINT" Division langer Zahlen" 

140 PRINT STRING*(30,42):PRINT 
150 : 

160 READ A*:PRINT"Dividend= ";A*:PRINT 
170 READ B*:PRINT"Divisor= " ; B* 

ISO : 

190 M=LEN(A*):N=LEN(B*):M=M-N 
200 IF M >=0 TUEN 220 

210 F'R I NT " Di vi dend < Divisor sein":END 
220 DIM U ( N+M ) , V (N> , W ( N ) , LI ( M ) 

230 : 

240 FOR 1=1 Tü N+M 

250 U(I)=VAL(MID*(A*,I,1)) 

260 NEXT I 

270 FÜR 1=1 TD N 

280 9(I)=9AL(MID*(B*„I,1)) 

290 NEXT I 
300 : 

310 D=INT(10/(V(1)+1)):U(0)=0 
320 IF D=1 THEN 450 
330 K=0 

340 FOR I=N+M TO 1 STEP -1 

350 U=U(I)*D+K:K=INT(U/10):U(I)=U-K*10 

360 NEXT I 

370 U(0)=K 

380 : 

390 K=0 

400 FOR I=N TO 1 STEP -1 

410 9=9 (I) *D+K: K= I NT (9/ 10) : V (I ) =9--K* 10 

420 NEXT I 
430 9(0)=K 
440 : 

450 J=0 

460 IF U(J)=9(1) THEN Q=9:G0T0 490 
470 IF 9(1)=0 THEN PRINT"führ-Null":END 
480 Q=INT((U(J)*10+U(J+1))/9(1)) 

490 Q1=(U(J)*10+U(J+l)-0*9(1))*10+U(J+2) 
500 IF 9(2)*0< =Q1 THEN 530 
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510 Q~Q—1:GOTO 490 
520 : 

530 K=0 

540 FÜR I=N TO 1 STEP -1 

550 W=V(I)*0+K:K=INT(W/10):W(I)=W-K*10 
560 NEXT I 
570 W(0)=K 
580 : 

590 K~0:F=0 

600 FÜR I--J+N TD J STEP -1 
610 U=U(I)-W(I-J)+K:K=INT(U/10) 

:U <I)=U-K*10 
620 NEXT I 

630 IF K=-1 TUEN F=1 
640 : 

650 Q(J)=Q 

660 IF F-0 TUEN 730 

670 Q(J)=Q(J)-1:K=0 

680 FÜR I=J+N TO J STEP -1 

690 U=U (I) + V (I - J ) +K: K= I NT (U/ 10) 

: U (I ) =U~K* 10 
700 NEXT I 
710 U (J-l) =U (J--1) +K 
720 : 

730 J=J+1 

740 IF J< =M THEN 460 
750 : 

760 K=0 

770 FÜR I=M+1 TG N+N 
780 U=U(I)+K*10:ü(I)=INT(U/D) 
:K=U-INT<U/D)*D 
790 NEXT I 
800 : 

810 REN Ergebnis 

820 PRINT:PRINT"Quotient=" :J=0 
830 IF Q(J)=0 THEN J=J+1:GOTO 830 
840 FOR I=J TO N 
850 PRINT 0(1); 

860 NEXT I:PRINT" 

870 FOR I=M+1 TO N+N 
880 PRINT U(I) ; 
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890 NEXT 
900 END 
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5. Potenzieren 


Eine weitere Anwendung der Mehr-Register-Arithmetik liefert das Potenzie¬ 
ren, da hier sehr schnell große Zahlen auftreten. 

Das Potenzieren wird im Programm als wiederholte Multiplikation durchge¬ 
führt. Eine andere Möglichkeit wäre das logarithmische Rechnen, jedoch ist 
die mehrfachgenaue Berechnung von Logarithmen noch aufwendiger. 

Teilt man die mit B zu multiplizierende Zahl in Achterblöcke R(l), so kann 
die Multiplikation wie folgt durchgeführt werden: 

100 U = 0 

110 FOR 1=1 TO R 

120 H = R(l) * B + U 

130 IF H<=-1E8THEN U = 0: GOTO 160 
140 U = INT (H/1E8) 

150 H = H - U * 1E8 
160 R(l) = H 
170 NEXT I 

Dabei ist R die Anzahl der Achterblöcke und U der jeweilige Überlauf. 

Soll die Zahl B zur N-ten Potenz erhoben werden, so muß das angegebene 
Programmstück noch in die Schleife 

FOR K=1 TO N.NEXT K 

eingebaut werden. 

Die benötigte Anzahl der Achterblöcke wird logarithmisch berechnet: Der 
Zehnerlogarithmus der gesuchten Potenz ßN ist 

N * LOG (Bl/LOG (10) 


Da der aufgerundete Zehnerlogarithmus einer Zahl die Stellenzahl angibt, 
kann daraus durch Division mit 8 die notwendige Registerzahl R berechnet 
werden. 
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Zum folgenden Programm 

Das Programm berechnet nach dem oben beschriebenen Verfahren beliebige 
Potenzen mit einstelliger Basis. 

Soll die Basis mehrstellig sein, so wird empfohlen die Registerzahl entspre¬ 
chend zu vermindern. 

Als Beispiel wird die hundertste Potenz von 2 berechnet. Es ergibt sich: 

1267650 60022822 94014967 032053756 
(vgl. Programmausdruck). 
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100 REM Berechnung grosser Potenzen 
1 10 : 

120 CLS:PRINT STRING*<35,42) 

130 PRINT" Berechnung hoher Potenzen" 

140 PRINT STRING*(35,42) 

150 : 

160 PRINT: INPLJT"Wel che Basis";B 
170 INPLJT"Zu welcher Potenz N 
1SO PR INT:PR I NT B;"hoch ";N;" = " 

200 DIM R(100) 

210 : 

220 'Zahl der Fuenferbloecke 
230 R -N*LC)G < B) / LOG < 10) 

240 R7.= INT(R/5) +1 
250 : 

260 ' Mul t; i p 1 i kat i onsschl ei f e 
270 11=0: R < 1 > =1 
280 FOR K=1 TO N 
290 FOR 1 = 1 TO R7. 

300 H=R(I)*B+U 

310 IF HO 100000 THEN ü=0: GOTO 330 

320 U=INT(H/100000 ):H=H-U*100000 

330 R(I)=H 

340 NEXT I 

350 NEXT K 

360 : 

370 'Ausgabe 

380 IF R (R7.) =0 THEN R7.=R7.~ 1: GOTO 380 
390 FÜR I=R7. TO 1 STEP -1 

395 R*=STR*(R(I)):R*=RIGHT*(R*,LEN(R*>-1) 

400 PR INT RIGHT*("0000"+R*,5); " " ; 

410 IF (R7.-I + 1) MOD 10 =0 THEN PRINT 
420 NEXT I 
430 END 
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B ö r e c h n u n q h o h e r P o t e n e n 
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6. FAKULTÄTS-BERECHNUNG 


Die Berechnung der Fakultäts-Funktion liefert ebenfalls eine Anwendung der 
Mehr-Register-Arithmetik. Die Fakultätsfunktion nl, definiert als das Produkt 
aller natürlichen Zahlen von 1 bis n, wächst nämlich noch schneller als eine 
Potenzfunktion. 

Die Fakultätsberechnung wird im Programm als wiederholte Multiplikation 
durchgeführt. Eine andere Möglichkeit wäre auch hier das logarithmische 
Rechnen, jedoch ist die mehrfachgenaue Berechnung von Logarithmen noch 
aufwendiger. Die Multiplikation erfolgt nach dem bei Programm 5 angegebe¬ 
nen Verfahren. 

Komplizierter ist nur die Berechnung der benötigten Register. Sie erfolgt 
ebenfalls logarithmisch. Da der Logarithmus eines Produkts gleich der Sum¬ 
me der Logarithmen ist, ergibt sich log(n!) aus 

100 FOR 1=2 TO N 
110 L=L+LOG(I) 

120 NEXT I 

Da die so ermittelte Zahl L, aufgerundet die Stellenzahl angibt, ist L/6 die 
benötigte Anzahl von Sechserblöcken. Da der Computer jedoch mit natürli¬ 
chen Logarithmen rechnet, muß die Zahl zuvor über den Faktor 

In (10) = .4342945 

auf den Zehnerlogarithmus umgerechnet werden. 


Zum folgenden Programm 

Das Programm berechnet die Fakultätsfunktion für dreistellige Zahlen in Fün¬ 
ferblöcken. Ist die Fakultät einer größeren Zahl gesucht, so wird aus Genauig¬ 
keitsgründen empfohlen mit Vierer- bzw. Dreierblöcken zu rechnen. Dazu 
müssen die im Programm auftretenden Sechsen entsprechend ersetzt wer¬ 
den. 

Als Programmbeispiel wird 

100 ! 

berechnet. Das Ergebnis kann dem Programmausdruck entnommen werden. 
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100 REM Fakultaet grosser Zahlen 
110 : 

120 CLS:PR INT STR INGT < 35,42) 

130 PRINT" Fakultaet grosser Zahlen " 

140 PR INT STR INGT(35,42) 

150 : 

160 PRINT:INPUT"Fakultaet welcher Zahl";N 
170 PRINT:PRINT N;"!=“ 

190 DIM R(100) 

200 ; 

210 'Zahl der Fuen-f er bl oecke 
220 L.= 1: R (1) = 1 
230 FÜR 1=2 TO N 
240 L=L+L0G(I) 

250 NEXT I 

260 L=L.*. 4342945: R7.=L/5+1 
270 : 

280 FÜR 1=2 TO R7. 

290 R(I)=0 
300 NEXT I 

310 : 

320 ' Mu 11 i p 1 i kat. i ansschielte 
330 L=1 

340 FOR I=N TO 2 STEP --1 
350 L=L+L_0G < I ) *. 434294575 

360 L7.=L/5+l 

370 U=0 

380 FOR J = 1 TO L7. 

390 H=R(J)*I+U 

400 IF H<-100000 TFIEN U=0:G0T0 420 

410 U=I NT(H/100000 ) :H=H-U*100000 

420 R(J)=H 

430 NEXT J 

440 NEXT I 

450 : 

460 'Ausgabe 

470 IF R (R7.)=0 TFIEN R7.=R7.-1: GOTO 470 
480 FOR I =R7. TO 1 STEP -1 

490 RT=STRT (RU)): RT=R I GH TT (RT , LEN (RT) -1) 

500 PR INT RI GHTT<"0000"+RT,5);" " ; 

510 IF (R7.-I + 1) MOD 8 =0 TUEN PRINT 
520 NEXT I 
530 END 
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(R7.-I + 1) MOD 8 =0 THEN PRINT 


o er- in o 
0 - 0 - o 
0-4 Kl O 
O M r-4 O 
O Kl N O 


N 0'- N O 
•43 0- IN O 
43 Ch CO O 
C-4 0' o o 
•43 m 04 o 


* * 

* * 

* * 

* * 

* * 

* * 

* * 

+ * 

* * 

* * 

* * 

* c * 

* tu * 

* ,-h * 

* JZ * 

* 03 * 

* r-4 * 

* * 

* L * 

* Qi * 

* C2 * 

* 0 * 

* L * 

* CT' * 

* * 

* +J * 

* :Hj * 

Li. 1— * -M * 

l-l X Q * r-t * 

LU Z * 3 * 

Z LU * -i * 

* 'Ti * 

o o o * u. * 

—I C-4 Kl * * 

in in in * * 



in 


0- 

«3 1 


CD 

r-H 

M- 



03 

C-4 

Cr 

Ti 


Ki 

in 

43 

» w « 


EN 

o-- 

Ki 

C ‘ 


O' 

03 

in 
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O 

K' 

C-4 

LJ 


•43 

■43 

■43 

« w « 


tH 

0- 

CO 

i^i 


CO 

C-4 

C4 

Q 

Q 

-43 

0- 

03 

0 

i’-'i 

C-4 

in 

r-H 

i_i 

rH 

in 

03 

in 

0 


tH 

■43 

43 
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r- 



in 

'S" 

rH 
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nj 


H 
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43 
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Ki 
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L 
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Ö- 

r-H 

Ql 
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03 

K< 
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1— 
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r—H 

in 

Ki 

43 

0 

0! 

•H 


er 

r-H 

2 

C-4 

■43 

tH 

r-4 


-0 

r-4 

'S" 

in 

-M 

C-4 

03 

0' 

03 

:rri 

1! 




-U 

— - 




r-H 

Ki 

43 

03 

r-H 

Ij 

O Kl 

O 

o 

r-H 

V 

O E> 

in 

q 

in 

•13 

rH '!^i 

r-H 

Ul 

r-4 

Ll 

o 

r- 

H 

03 
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7. e AUF 1000 STELLEN 


Mit der Zahl n ist die Eulersche Zahl 
e = 2.7182318 ... 
über die Gleichung 
e m = _i 

verbunden, dabei ist i = V - 1 • 

e ist insbesondere als Basis der natürlichen Logarithmen von Bedeutung 
y = In x => x = eY 
Ebenso stellt e den Grenzwert 


e = lim <1 + l)x 
x -*■ °° x 

dar. Mit diesem Grenzwert ist auch die Zinzeszinsformel verknüpft: Wird 
1 DM nicht jährlich oder monatlich, sondern kontinuierlich verzinst, so 
wächst der Zinseszins bei 100% nicht ins Unbeschränkte, sondern nähert sich 
dem Wert e DM. Die Reihenentwicklung der Exponentialfunktion liefert eine 
für numerische Zwecke brauchbare Formel für e 


_ 1 , 1 . 1 . 1 . 1 , 1 . 

1! 2! 3! 4! 5! 6! 


dabei stellt n! die Fakultätsfunktion dar. 

Erwähnenswert ist, daß auch e als Ergebnis eines Zufallsexperiments auftritt. 
Werden n Dinge zufällig angeordnet, so ist mit der Wahrscheinlichkeit 

1/e = 0.3678794. 

kein Gegenstand an seinem Ort. Solche Zufallsexperimente können mit Hilfe 
von Programm 12 erhalten werden. Der Versuch ist besser bekannt als das 
Problem der "vertauschten Briefe". Legt man n Briefe zufällig in n vorbereite¬ 
te Umschläge, so liegt mit der Wahrscheinlichkeit 

1 -1/e =0.6321205 

mindestens ein Brief im richtigen Umschlag. Im Band 1 der vorliegenden Pro¬ 
grammsammlung findet sich für n = 5 ein Grafikprogramm. 
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Zum folgenden Programm 

Das Programm berechnet die Dezimalbruchentwicklung von e mit Hilfe der 
angegebenen Reihenentwicklung. Die Kehrwerte der benötigten Fakultäten 
werden durch wiederholte Divisionen berechnet und fortlaufend aufsummiert. 

Das hier angegebene Programm läßt sich auf beliebige Stellenzahl ausdehnen. 
Die Variable N in Zeile 220 muß dazu so geändert werden, daß die Dezimal¬ 
stellen von 1/N! bei der gewünschten Genauigkeit vernachlässigt werden kön¬ 
nen. 
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100 REM e auf 1000 Stellen 
110 : 

120 CLS:PRINT STRING*(30,42) 

130 PRINT " e auf 1000 Stellen" 

140 PRINT STRING*(30,42) 

150 : 

160 DII V I A(202)Zahl der Fuenf erbl oecke 
170 FÜR 1 = 1 TO 202 
180 A(I)=0 
190 NEXT I 
200 : 

210 B = 100000 : A (0) = 1 

220 FÜR N-450 TO 1 STEP -1 

230 FOR 1=0 TO 201 

240 IF A(I)=0 THEN 280 

250 Q= I NT (A (I) /N > : R=A ( I) --Q*N 

260 A(I)=Q:A(I+1)=A(I+1)+B*R 

270 NEXT I 

280 A(0)=A(0)+1 

290 NEXT N 

300 : 

310 'Rundung 

320 A (201) =Ä (201) +1 NT (A (202.) / B+. 5) 

330 FOR 1=200 TO 1 STEP -1 

340 U=INT(A(I)/B)s A(I)=A(I)~B*U 

350 A(I-1)=A(I — 1)+U 

360 NEXT I 

370 : 

380 PEIIMT"e=" ; A (0) ; " . " 

390 FOR 1=1 TO 200 

400 A*=STR*(A(I)):A$=RIGHT*(A*,LEN< A* >-1) 

410 PR I NT RIGHT t- ( " 0000 " +A* , 5) ; " " ; 

420 IF I MOD 8 =0 THEN PRINT 
430 NEXT I 
440 END 


45 




****************************** 
e auf 1000 Stellen 
■*******■*■******■***•***•**■****■**■** 


71828 

18284 

59045 

66249 

77572 

47093 

62772 

40766 

30*^53 

52516 

64274 

27466 

81741 

35966 

29043 

59563 

07381 

32328 

82988 

07531 

95251 

93070 

21540 

89149 

76146 

06680 

82264 

74234 

54424 

37107 

55170 

27618 

38606 

75204 

49338 

26560 

71742 

64364 

22912 

31264 

f f /| er -j 

49669 

74639 

85002 

03510 

22517 

64866 

71751 

80122 

36490 

02456 

99960 

71749 

18641 

37605 

56454 

08460 

43587 

90914 

02547 

32649 

37577 

68259 

54657 

04140 

40868 

13658 

42489 

90958 

74336 

12838 

01761 

00965 

92466 

98058 

41272 

07192 

10263 

98035 

67741 

57062 

36086 

65511 

16056 

31994 

92993 

48212 

19923 

37104 

37136 

65791 

61003 

20703 

92931 

94953 

77431 

14315 

74683 

58961 

49749 

84728 



23536 

02874 

71352 

69995 

95749 

66967 

54759 

45713 

82178 

39193 

20030 

59921 

57290 

03342 

95260 

62794 

34907 

63233 

01901 

15738 

34187 

93488 

41675 

09244 

80016 

84774 

11853 

53907 

77449 

92069 

26133 

13845 

83000 

28864 

68800 

72705 

59039 

54818 

65870 

34792 

29434 

01033 

59931 

78084 

20832 

87375 

78427 

63076 

95587 

91367 

18709 

04258 

68953 

74226 

42957 

75390 

51037 

48265 

73137 

70140 

40475 

49435 

58471 

51052 

63327 

79621 

78661 

05465 

30091 

73296 

06216 

23923 

07785 

54983 

78516 

06389 

90530 

75586 

08405 

90325 

33337 

38466 

32017 

85929 

68436 

63182 

49742 

1 1443 

27646 

39803 

24716 

31157 

66078 

41698 

42363 

91701 

40066 

05718 

45316 
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8. LINEARE DIOPHANTISCHE GLEICHUNG 

Ein Standardverfahren zur Lösung einer linearen Diophantischen Gleichung 
ax + by = c 

gleichbedeutend mit der Kongruenz 
ax = c mod y 

ist der erweiterte Euklidsche Algorithmus. Mit seiner Hilfe kann der größte 
gemeinsame Teiler 

ggT (a,b) 

als Linearkombination von a und b bestimmt werden. 


Dies soll an einem Beispiel für a = 286 und b = 121 gezeigt werden: 


262 

121 =2 Rest 44 => 

286 = 2 • 121 +44 

121 

44 = 2 Rest 33 => 

121 = 2-44 + 33 

44 

33 = 1 Rest 11 => 

44 = 1 • 33 + 11 

33 

11=3 => 

33 = 3 ■ 11 


Somit gilt 

ggT (286, 121)= 11 
Durch Rückwärtsrechnen folgt 

11 = 44 - 33 

= 44 -(121 — 2 - 44) 

= 3-44-121 
= 3- (286-2 • 121) - 121 
= 3-286 -7-121 

Dies liefert die gesuchte Linearkombination 
3 - 286 - 7 - 121 = 11 
allgemein 

ax + by = ggT (a,b) 

Sind x Q , y 0 die Koeffizienten der Linearkombination, so hat die Diophanti- 
sche Gleichung 

ax + by = c 

die allgemeine Lösung: 


x = cx 0 + 


a 

ggT (a,b) * 


y = cy 0 - 


ggT (a,b) 


- 1 ; 


t ganzzahlig 
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Diese Lösung existiert nur, wenn c ein Vielfaches von ggT (a,b) ist; andern¬ 
falls ist die Gleichung nicht lösbar (vgl. [7]). 


Zum folgenden Programm 

Das Programm berechnet mit Hilfe des erweiterten Euklidschen Algorithmus 
nach dem oben angegebenen Verfahren eine spezielle Lösung der Gleichung 

ax + by = ggT (a,b). 

Damit wird die allgemeine Lösung bestimmt. 

Das im Programm angegebene Verfahren kann prinzipiell auch für Diophanti- 
sche Gleichungen der Form 

ax — by = c 

(b > C) angewandt werden. In diesem Fall ergibt sich für den ggT ein negativer 
Wert, im Gegensatz zur üblichen Definition; zwei aufeinander stoßende nega¬ 
tive Vorzeichen sind als positiv zu lesen. 


Als Beispiel wird die Diophantische Gleichung 
17x + 6y = 5 

gelöst. Das Programm liefert die Parameterlösung 

x = —5 —6t 
y = 15 +17t 

dabei durchläuft t alle ganzen Zahlen. Für jeden Wert von t ergibt sich wieder 
eine spezielle Lösung, z.B. für 

t = 0 : x = —5, y = 15 
t = 1 : x = —11, y = 32 
t = -1 : x = 1, y = —2 

usw. Zahlreiche Diophantische Gleichungen findet man in eingekleideter 
Form in den historischen Aufgaben des Anhangs. 
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100 REM Lineare Diophantisehe Gleichung 
1 10 : 

120 CLS:PRINT STRING$(30,42) 

130 PRINT" Lin.Diophant.Gl eichung" 

140 PRINT STRING#(30,42) 

150 : 

160 DEFINT A—C , X—Y 

170 PRINT"a*x + b*y = c" 

180 INPUT"Eingabe a,b,c";A,B,C 
190 : 

200 REM Initialisieren 
210 A0=A:A1=B:X0=1 
220 X1=0:Y0=0:Y1=1 
230 : 

240 REM Erweiteter Euklidscher Algorithmus 
250 Q=INT(A0/A1) 

260 IF Q=A0/A1 THEN 320 

270 A2=A0-Q*A1:X2=X0-Q*X1:Y2=Y0~Q*Y1 

280 A0=A1:A1=A2:X0=X1:X1=X2 

290 Y0=Y1:Y1=Y2:GOTO 250 

300 : 

310 REM 

320 PRINT"ggT(";A;",";B;")=";A2:PRINT 
330 F'RINT"Loesung der Gleichung"; 

340 PRINT" a*x+b*y=ggT(a,b)" 

350 IF Y2<0 THEN 370 

360 PRINT TAB(8)A;;X2;;B;; Y2; 

"=";A2:GOTO 380 

370 PRINT TAB(8)A;X2,"-";B; 

ABS(Y2>;"=";A2 

380 : 

390 IF C/A2<>INT(C/A2) THEN PRINT 
"existiert nicht":END 
400 PRINT 

410 PRINT"Al 1g.Loesung der Diophant. 

Gl eichung:" 

420 PRINT TAB(8)"x=";C*X2/A2;"—t*"B 
430 PRINT TAB(8)"y=";C*Y2/A2;"+t*";A; 

440 PRINT SPC(5)"t ganzzahlig" 

450 END 
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****************************** 

Lin.Diophant.Gl eichung 
****************************** 

a*x + b*y = c 
Eingabe a,b,c? 17,6,5 

ggT( 17 , 6 )= 1 

Lösung der Gleichung a*:;+b*y=ggT (a , b ) 
17 *-1+6* 3=1 

Al lg.Lösung der Diophant. 

Gl eichung: 

5 -t* 6 

y- 15 +t* 17 t ganzzahlig 
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9. KETTENBRUCH / BRUCHAPPROXIMATION 


Sucht man zu einem vorgegebenen Dezimalbruch einen Quotienten gleichen 
Werts — z.B. ein Windungszahlverhältnis beim Transformator — so kann dieser 
mit Hilfe eines Kettenbruchs bestimmt werden. 


Ein Kettenbruch hat folgende Form 


e= 2 + 


1 

1 + 1 
2 + 2 
3 + 3 
4 + 4 


Solche Kettenbrüche können für rationale, d.h. für Bruchzahlen durch den 
Euklidschen Algorithmus berechnet werden. 

Für 67/29 ergibt sich z.B. 

67 : 29 = 2 Rest 9 
29 : 9 = 3 Rest 2 
9 : 2 = 4 Rest 1 
2 : 1=2 


Somit folgt 


— = 2 + — = 2 + — 
29 29 29 


29 „ 2 „ 1 

r 3+ 9= 3+ 9 


9 , 1 
2 ^ + 2 


Setzt man diese Gleichungen rückwärts ineinander, so ergibt sich der Ketten¬ 
bruch 



1 

3+ 1 


4 


1 

+ — 
2 
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Der Euklidsche Algorithmus zeigt, daß der Kettenbruch einer rationalen Zahl 
abbricht, jedoch nicht der einer irrationalen Zahl. Ein Beweis dazu findet sich 
in [14], Kettenbrüche können somit durch folgendes Programm berechnet 
werden 

100 INPUT "DEZIMALBRUCH"; X 
110 PRINT INT (X) 

120 IF X=INT(X) THEN 150 
130 X=1/(X—INT(X)) 

140 GOTO 110 
150 END 

jedoch muß dabei noch die beschränkte Genauigkeit des Rechners beachtet 
werden. 


Zum folgenden Programm 

Das folgende Programm berechnet zu einem eingegebenen Dezimalbruch maxi¬ 
mal 8 Teilnenner nach dem oben angegebenen Programmstück. Dabei wird 
versucht Rechnerungenauigkeiten, so weit wie möglich, abzufangen. 

Der so berechnete Kettenbruch wird auf dem Bildschirm ausgegeben und die 
zugehörige rationale Zahl berechnet. 

Als Programmbeispiel wurde 
rr = 3.14159265 

eingegeben. Man erhält die Bruchapproximation 
2551690445 
80115557 

Der ausgegebene Kettenbruch kann dem Programmausdruck entnommen wer¬ 
den. 
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100 REM Kettenbruch/Bruchapproximation 
110 : 

120 CL.S: PRINT STRING# (35,42) 

130 PRINT" Kettenbruch/Bruchapproximation" 
140 PRINT STRING#(35,42) 

150 : 

160 PRINT:INPUT"Dezima1z ah1";Y 
170 DIM R(12) 

180 : 

210 K=8 

220 FOR 1=1 TG K 

230 R(I)=INT(Y+.000001) 

240 IF ABS (R(I)-YX. 000001 THEN 280 
250 Y=1/(Y—R ( I ) ) 

260 NEXT I:GOTO 290 
270 : 

280 K=I 

290 PRINT R(1);" + 1" 

300 FÜR 1=2 Tq'k-1 

310 PRINT TAB(I*6-7);" -" 

320 PRINT TAB<1*6-6);R<I);" + 1" 

330 NEXT I 

340 PRINT TAB <K*6—7);“-" 

350 PR INT TAB(K*6-6);R <K> 

36O : 

370 Z=R(K):N=1 

380 FOR J=K—1 TG 1 STEP -1 

390 H=Z:Z=R (J)*Z+N:N=H 

400 NEXT J 

410 PRINT Z 

420 PRINT"“";STRING#(LEN(STR# (Z) ),45); 

430 PRINT Z/N 
440 PRINT N 
450 END 
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Kettenbruch/Bruchapprox i mat. i on 
**************************** ******* 

Dezimalzahl? 3.14159265 
3 + 1 

7 + 1 


15 +1 

1 + 1 


288 + 1 


■+• 1 

4 + 1 


251690445 

=--= 3.14159265 

80115557 
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10. PERMUTATIONEN 


Da endliche n—elementige Mengen stets auf die Menge 



abgebildet werden können, arbeiten die folgenden Kombinatorik-Programme 
stets mit dieser Anfangsmenge der natürlichen Zahlen. 

Jede Anordnung der Menge 1,2, 3,.... n, bei der jedes Element genau einmal 
vorkommt, heißt Permutation. 

Permutationen treten vielfach in der Kombinatorik und Wahrscheinlichkeits¬ 
rechnung auf. Man benötigt sie insbesondere auch in Operations Research, 
wenn man z.B. auf der Suche nach der kürzesten Rundreise alle Anordnungen 
der zu besuchenden Städte durchmustern muß. 

Es gibt 

n! 

Permutationen einer n-elementigen Menge;dabei ist n! die in Programm 6 be¬ 
handelte Fakultätsfunktion. 

Die Fakultäten der Zahlen wachsen sehr schnell, so gilt 

6! = 720, 7! = 5040, 8! = 40320 
91 = 362880, 101 = 3628800 
11! = 39916800 

Somit gibt es für 10 Leute bei 10 Sitzplätzen 
3628800 

verschiedene Sitzordnungen. Spielt in einer Fußballmannschaft jeder auf jeder 
Position, so gibt es 

39916800 

verschiedene Mannschaftsaufstellungen. 


Zu folgendem Programm 

Das folgende Programm erzeugt nach einem Algorithmus von Ord-Smith die 
Permutation einer n—elementigen Menge in lexikographischer Anordnung. 

Dabei heißen zwei Zahlenfolgen aj und bj lexikographisch geordnet, wenn für 
die erste Stelle mit 


gilt 


aj *= bj 
aj<bj 


Der Programmausdruck zeigt die Permutation einer 4—elementigen Menge. 
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100 REM F'ermutat.i onen 
110 : 

120 CLS:PRINT STRINGS(25,42) 

130 PRINT“ F'ermutationen" 

140 PRINT STRING*(25,42) 

150 : 

160 INF'UT" Wi evi el e Elemente 1.. n" ;N: PRINT 
170 DIN X < N) ,Y(N) 

180 " 

190 REM Initialisieren 

200 FÜR 1=1 TO N 

210 X (I)=I 

220 NEXT I 

230 F=0:F-0:N1=N -1 

240 : 

250 REM wiederholen bis» Flag -f wieder =0 
260 FOR 1=1 TO N 

270 PRINT X(I); 

280 NEXT I:PRINT 

290 P=P+l 

300 GOSUB 350:' Unterprogramm Permutation 

310 IF F=1 THEN 250 

320 PRINT P;"Permutationen" 

330 END 
340 : 

350 REM Rekursive Erzeugung der Permutationen 
360 IF F=1 THEN 410 
370 F=1 

380 FOR K=1 TO N 
390 Y(K)=N 

400 NEXT K 

410 IF Y <N1)< >N THEN 450 
420 Y(N1)=N1 

430 H=X(N):X(N)=X(N1):X(N1)=H:GOTO 600 
440 : 

450 FOR J=1 TO NI 
460 K=N—J 

470 IF Y(K)< >K THEN 520 

480 Y(K)=N 

490 NEXT J 

500 F=0:K=1:GOTO 550 
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510 : 

520 M=Y(K) 

530 H=X(M):X(M)=X(K):X(K)=H 
540 Y(K)=M-1:K=K+1 
550 M=N 
560 : 

570 H=X <M> : X (M) =X <K> : X <K) =H 
580 M=M-1:K=K+1 
590 1F K<M TUEN 570 
600 RETURN 

******■*#■#•*•*•*•**■*•*■*•**•**•**■** 
Permutationen 
*******■»■»*■»** ■»***■*■***«■*** 
Wieviele Elemente .1. . n? 4 

12 3 4 

1 2 4 3 

13 2 4 

1 3 4.2 

14 2 3 

14 3 2 

2 1 3 4 

2 1 4 3 

2 3 1 4 

2 3 4 1 

2 4 1 3 

2 4 3 1 

3 12 4 

3 14 2 

3 2 14 

3 2 4 1 

3 4 1 2 

3 4 2 1 

4 12 3 

4 13 2 

4 2 13 

4 2 3 1 

4 3 12 

4 3 2 1 

24 Permutationen 
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11. ZUFALLSPERMUTATION 


Für viele Zwecke benötigt man zufällige Permutationen, z.B. für die bei Pro¬ 
gramm 8 erwähnten Monte-Carlo-Simulation der vertauschten Briefe. 

Als Anfangswert wird die Permutation P(I) der Zahlen 1, .n auf die natürli¬ 
che Reihenfolge gesetzt. 

Sodann werden n—1 ganzzahlige Zufallszahlen ausgelost und die Permutation 
mit entsprechendem Index ausgetauscht. Dies wird durch folgendes Pro¬ 
grammstück durchgeführt. 

260 FOR J=N TO 2 STEP —1 
270 R = INT (J * RND(D) + 1 
280 H = P(J) : P(J) =P(R) : P(R) = H 
290 NEXTJ 

Die so erhaltene Zufallspermutation wird ausgedruckt und bei Bedarf eine 
neue erzeugt. 


Zum folgenden Programm 

Nach Abfrage, wieviele Zahlen die Permutation umfassen soll und wieviele 
Permutationen erwünscht sind, werden nach dem angegebenen Verfahren mit 
Hilfedeseingebauten Zufallszahlengenerators die benötigten zufälligen Permu¬ 
tationen erzeugt. 

Der Programmausdruck zeigt 20 Zufallspermutationen der Zahlen 1 bis 9. 
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100 REM Zufallspermutation 

110 : 

120 CLS:PRINT STRINGS<25,42) 

130 PRINT" Zufal1spermutationen" 

140 PRINT STRING*(25,42) 

150 : 

160 INPUT"Wieviele Zahlen";N 

170 I NF'UT " Wi evi el e Permutat i onen " ; M: PR I NT 

180 DIN P(N) 

190 : 

200 REM Initialisieren 
210 1 = 1 
220 : 

230 RANDOMIZE TIMER 
240 WHILE I< =M 
250 FÜR J=1 TO N 
260 P<J)=J 

270 NEXT J 
280 : 

290 FOR J —N TO 2 STEP -1 

300 R=INT < J*F;ND> +1: 'Auslasen der Reihenfolge 

310 H=P<J)s P <J)=P(R):P <R)-H 

320 NEXT J 

330 FOR J = 1 TO N 

340 PRINT F' < J ) ; 

350 NEXT J:PRINT 
360 1=1+1 

370 WEMD 
380 END 
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************•*■**•*•*•*•*•**•*•*•#••* 
Zufal1spermutationen 
■****************•*■*•**■•*•**•** 

Wieviele Zahlen ? 9 
Wieviele Permutationen ? 20 


2 3 8 6 9 

4 3 5 6 7 

8 3 4 5 9 

5 4 2 6 7 

1 3 8 6 5 

8 5 6 7 9 

4 8 5 2 7 

8 3 9 1 2 

1 4 5 6 9 

4 5 8 3 7 

5 7 2 9 4 

1 9 2 4 5 

7 3 4 6 8 

7 3 8 6 9 

7 15 4 3 

4 6 7 1 9 

6 4 5 9 8 

8 4 7 5 9 

9 3 2 4 1 

9 8 5 2 6 


4 

1 

9 

7 

9 


6 

6 

6 

9 

5 

6 
8 


7 

7 


8 

1 

8 


6 

7 


8 

5 

4 

8 

7 

8 
4 


1 

9 

6 

1 

9 

1 

1 

4 

7 
9 

8 
7 

9 

r> 

i 


7 

2 

7 

ö 

4 

4 

£> 

8 

1 

1 

1 

1 


1 

6 

6 

1 
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12. KOMBINATIONEN 


Ordnet man nur k Elemente einer n-elementigen Menge ohne Wiederholung 
an, so erhält man eine k-Kombination oder auch Kombination k.-ter Klasse. 

Die Anzahl der k-Kombinationen einer n-elementigen Menge schreibt man 



und bezeichnet dies als Binomialkoeffizient. Diese tragen ihren Namen daher, 
weil sie als Koeffizienten in den binomischen Formeln, z.B. 



= a 3 + 3a 2 b + 3ab 2 + b 3 


auftreten. 


Die Binomialkoeffizienten 


0 , 1 , 2 , 


(n+1)-ten Zeile des Pascalschen Dreiecks: 


. ., n finden sich jeweils in der 


1 


1 

1 1 

1 2 1 
13 3 1 

1 4 6 4 1 

5 10 10 5 1 

6 15 20 15 20 1 


Die Binomialkoeffizienten können auch über die Fakultätsfunktion berechnet 
werden 

n! 


kl (n—k)! 

So kann beim Schafkopf (24 Karten, 3 Spieler) jeder Spieler 
= 735471 
Blätter haben. Beim Lottospiel gibt es 
= 13983816 

verschiedene Lottotips, beim neuen Mittwochslotto entsprechend 
= 12620256 
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Zum folgenden Programm 

Das Programm erzeugt nach einem Algorithmus von Kurtzberg alle k-Kombi- 
nationen (oder k-ten Klasse) einer n-elementigen Menge. 

Der Programmausdruck zeigt die 
4-Kombinationen 
einer 6-elementigen Menge. 




£yn ißewe 

^nno ivolgcartfnotc 
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(CetjQfcfccr nod) in VDcI|cf)tr |t»ad) nie 
geöiiicFc. t*urcf) petrmtt 2lpi<um 
» von fi.ev/?nicF/tf 2l|Fronomei 
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"Titelblatt der Arithmetik von Petrus Apianus, Ingolstadt 1527 
mit der ersten gedruckten Darstellung des Pascalschen Dreiecks" 
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100 REM Kombinatioen 
110 : 

120 CLS:PRINT STRING*(25,42) 

130 PRINT" Kombinationen" 

140 PRINT STRING*(25,42) 

150 : 

160 INPUT"Wieviele Elemente";N 
170 INPUT"Zu welcher Klasse";M 
180 DIM X(M) 

190 : 

200 REM Initialisieren 

210 FDR 1=1 TO M 

220 Xd) =N 

230 NEXT I 

240 L=0:PRINT 

250 : 

260 REM Bi nomi al koe-f -f i z i ent m aus n 
270 P=1 

280 IF M=0 THEN 340 
290 FÜR 1=1 TO M 
300 P=P*(N-I+l)/I 

310 NEXT I 

320 : 

330 WHILE L<P 

340 GOSTJB 440: 'Kombinationen 

350 FOR 1=1 TO M 

360 PRINT X(I); 

370 NEXT I:PRINT 

380 L=L+1 

390 WEND:PRINT 

400 PRINT P;"Kombinationen" 

410 END 
420 : 

430 REM Rekursion 
440 FOR K=1 TO M 
450 J=K—1:A=M~J:B=N—J 

460 IF X(A)<B THEN 530 

470 NEXT l< 

480 : 

490 FOR K=1 TO M 
500 X(K)=K 
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510 NEXT l<: GOTO 570 
520 : 

530 B~X(A) 

540 FÜR K=A TO M 
550 B=B+1:X(K)=B 

560 NEXT K 
570 RETURN 


******** *********** ***** * 
Kombinationen 
************************* 
Wieviele Elemente? 6 
Zu welcher Klasse? 4 

12 3 4 

12 3 5 

12 3 6 

12 4 5 

12 4 6 

12 5 6 

13 4 5 

1 3 4 6 

13 5 6 

14 5 6 

2 3 4 5 

2 3 4 6 

2 3 5 6 

2 4 5 6 

3 4 5 6 

15 Kombinationen 
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13. VARIATIONEN 


Ordnet man k Elemente einer n-elementigen Menge unter Berücksichtigung 
der Reihenfolge an, so erhält man eine k-Variation oder Variation zu k-ten 
Klasse. 

Es gibt 

n! 

(n-k)l 

k-Variationen einer n-elementigen Menge. 

Wollen von 20 Personen 8 Leute an einem bestimmten Tisch sitzen, so gibt 
es dort 

— = 5079110400 
12! 

verschiedene Sitzordnungen. Aus einem Gremium von 30 Personen gibt es 

^ = 24360 
27! 

verschiedene Möglichkeiten einen Präsidenten, dessen 1. und 2. Stellvertreter 
auszuwählen. 


Zum folgenden Programm 

Das Programm erzeugt alle k-Variationen oder Variationen zur k.-ten Klasse 
einer n-elementigen Menge. 

Der Programmausdruck zeigt alle 3 Variationen einer 4-elementigen Menge. 
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100 REM Variationen 
110 : 

120 CLS:PRINT STRING*<25,42) 

130 PRINT" Variationen" 

140 PRINT STRING*(25,42) 

150 : 

160 INPUT "Wi evi el e Elemente";!'! 

170 INPUT"zu welcher Klasse";K:PRINT 
180 DIM X(N),Y(K) 

190 : 

200 REM Initialisieren 

210 V=0:M=K 

220 IF K=N THEN M=M-1 

230 FÜR 1=1 Tü N 

240 X <I)=I 

250 NEXT I 

260 FÜR 1=1 TO M 

270 Y(I)=I 

280 NEXT I 

290 ; 

300 I=M 

310 FOR J=1 TO K 
320 PRINT X(J); 

330 NEXT J:PRINT 

340 V=V+1 

350 WHILE I>0 

360 IF Y(I)>=N THEN 400 

370 Y(I)=Y(I)+1 

380 H=X(I):X <I)=X <Y <I)):X <Y (I) )=H:GOTO 300 
390 : 

400 H=X(I):X(I)=X(Y<I)):X(Y(I))=H 

410 Y(I)=Y(I)— 1 

420 IF Y(I)>I THEN 400 

430 1=1-1 

440 WEND 

450 PRINT:PRINT V; "Variationen" 

460 END 
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******* ****** ************ 
Variationen 

************************* 
Wieviele Elemente ? 4 
zu welcher Klasse ? 3 



24 Variationen 
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14. 01-TUPEL 


Wählt man k Elemente aus einer n-elementigen Menge aus und läßt dabei be¬ 
liebige Wiederholung zu, so erhält man die 

k-tupel 

der Menge. Im Fall von k=2 spricht man auch von Paaren, bei k=3 von Tripel. 

Es gibt insgesamt 
n k 

k-Tupel einer n-elementigen Menge. 

Besondere Bedeutung haben die 8-Tupel der Binärzahlen 0 und 1, Byte ge¬ 
nannt, für Computer. Alle ganzen Zahlen werden mit Hilfe zweier Bytes dar¬ 
gestellt. Da das erste Byte das Vorzeichen angibt, können somit INTEGER- 
Zahlen mit 15 Bit dargestellt werden. Es sind dies die Zahl 0 bis 32767. 

Weitere Anwendungen finden die 01-Tupel, da mit ihrer Hilfe alle Teilmengen 
einer gegebenen Menge erzeugt werden können: Ist die i.-te Stelle des 01- 
Tupel Eins, so gehört das Element i zur entsprechenden Teilmenge, sonst 
nicht. 

Somit hat eine n-elementige Menge soviele Teilmengen, wie der 01-Tupel der 
Länge n gibt, nämlich 

2 n 


Wichtig sind 01-Tupel auch in der Boole'schen Algebra und Aussagenlogik, da 
mit ihrer Hilfe alle Wahrheitswerte von Aussageverknüpfungen ermittelt wer¬ 
den können (siehe Logikprogramme im Band 1 der vorliegenden Programm¬ 
sammlung). Allerdings wird im IBM-BASIC "wahr" nicht als 1, sondern als 
—1 codiert; dies kann aber im Programm durch eine Vorzeichenänderung er¬ 
reicht werden. 


Zum folgenden Programm 

Das Programm erzeugt nach einem Algorithmus von Berztiss alle 01-Tupel der 
gewünschten Länge. 

Läßt man in Zeile 270 statt X(K) das Element A(K) einer Menge — falls 
X(K) = 1 — ausdrucken, so erhält man alle Teilmengen dieser Menge. 

Der Programmausdruck zeigt alle 01-Tupel der Länge 5. 
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"Entwurf von Leibniz für eine Gedenkmünze, die das von ihm gefundene 
Binärsystem darstellt (1697)" 
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100 REM 01—Tupel 
110 : 

120 CLS:PRINT STRINGS(25,42) 

130 PRINT" 01—Tupel" 

140 PRINT STRING*<25,42) 

150 : 

160 PRINT:INPUT"Laenge der 01—tupel";N 
170 DIM X <N+1) 

180 : 

190 REM Initialisieren 
200 FOR K=1 TO N 
210 X(K)-0 
220 NEXT K 

230 J=INT(2"N+.00001):L=1 
240 : 

250 WHILE L< =J 
260 FOR K=1 TO N 

270 PRINT X <K); 

280 NEXT K:PRINT 
290 GOSUB 350:'Unterprogramm 
300 L=L+1 

310 WEND 

320 PRINT:PRINT J;" Ol-Tupel" 

330 END 
340 : 

350 REM rekursive Erzeugung der Ol-Tupel 
360 1=1 

370 IF X(DOO TUEN 390 
380 X < I) = 1:GOTO 400 
390 X(I)=0:1=1+1:GOTO 370 
400 RETURN 
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************************* 
01-Tupel 

************************* 
Länge der 01—tupel ? 5 

0 0 0 0 0 

1 0 0 0 0 

0 1 0 0 0 

1 1 0 0 O 

0 0 1 0 0 

1 O 1 0 0 

0 1 1 0 0 

1110 0 
0 0 0 1 0 

1 0 0 1 0 

0 1 0 1 0 

110 1 0 
0 0 1 1 o 

10 11 0 
0 111 0 
11110 
0 0 0 0 1 

1 0 0 0 1 

0 1 0 0 1 

1 1 0 0 1 

0 0 1 0 1 

1 0 1 0 1 

0 110 1 
1110 1 
0 0 0 1 1 

1 0 0 1 1 

0 10 11 
110 11 
0 0 111 
10 111 
0 1111 

32 01—T upel 
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15. PARTITIONEN 


Unter einer Partition der Zahl n versteht man die Zerlegung von n in Sum¬ 
manden, wobei es auf die Anordnung nicht ankommt. So hat die Zahl 6 
folgende 11 Partitionen: 

6 = 6 
= 5+1 
= 4 + 2 
= 4 + 1+1 
= 3 + 3 
= 3 + 2 + 1 
= 3 + 1 + 1 + 1 
= 2 + 2 + 2 
= 2 + 2 + 1 +1 
= 2 + 1 + 1 + 1 + 1 
= 1 + 1 + 1 + 1 + 1+1 

Solche Partitionen werden benötigt bei Geldwechsel- und Frankatur-Proble¬ 
men. 


Zum folgenden Programm 

Das Programm erzeugt sämtliche Partitionen einer Zahl nach einem Algorith¬ 
mus von McKay. 

Der Programmausdruck zeigt die 42 Partitionen der Zahl 10. 
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100 REM Partitionen 
110 : 

120 CLB:PRINT STRING*(25,42) 

130 PRINT" Partitionen" 

140 PRINT STRING*(25,42) 

150 : 

160 INPUT"Welche Zahl";N 
170 DIM P(N) 

180 PRINT"Die Partitionen von";N;" sind:" 
190 : 

200 REM Initialisieren 
210 P(1)=N:K=1:F=0:Z=0 
220 : 

230 REM Wiederholen bis Flag f=0 
240 GOSUB 330:'Unterprogramm 
250 FOR J=1 TD K 
260 PRINT P(J) ; 

270 NEXT J:PRINT 
280 Z = Z+1 

290 IF F=1 THEN 240 

300 PRINT Z;" Partitionen" 

310 END 
320 : 

330 REM Rekursion 
340 IF F=1 THEN 410 
350 F=1 

360 FOR 1=1 TO K 

370 IF P(I)=1 THEN 480 

380 NEXT I 

390 I=K:GOTO 480 

400 : 

410 L=K-1:K=I:P <I)=P(I)-1 
420 IF P(K)>L THEN 460 
430 L=L—P(K):K=K+1 
440 P(K)=P(K-1):GOTO 420 
450 : 

460 K=K+1:P <K)=L+1 
470 IF P(I)<>1 THEN I=K 
480 IF P(I)=1 THEN 1=1-1 
490 IF 1=0 THEN F=0 
500 RETURN 
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*** ****■*•**■*•*■* **'* **■**■#•*•*** 
Partitionen 

***************)(********* 
Welche Zahl ? 10 


Di e 
10 
9 
8 
8 
7 
7 


6 

6 

6 

6 

6 


o 

5 

5 

5 

4 

4 

4 

4 

4 

4 

4 

4 

4 


o 

1 


Partiti 


1 

'”1 

1 

1 

4 


1 

5 

4 


1 

4 

4 


1 


1 

*■> 


O 

1 

1 


1 

1 

1 

1 

'■? 

1 

1 

1 

1 

o 

1 

1 

9 

i 


i 


i 

i 


i 


i 

i 


2 

1 

1 

1 


1 


1 

1 


1 

1 

1 

1 

1 

1 

1 

n 

l 

l 

l 

i 

>-) 

l 

l 

l 

l 

l 


l 

l 

l 

l 


onen 


1 


1 

1 


1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

*9 

i 

i 

i 

i 

i 


von IO sind 


1 

1 1 


1 

1 

1 1 
1 1 1 

1 

1 1 
1 1 1 
1111 
1111 


42 Partitionen 
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16. KUBISCHE GLEICHUNG 


Die kubische Gleichung 

ax^ + bx^ + cx + d = 0 

kann mit Hilfe der Cardano-Formel gelöst werden. Die Cardano-Formel ist 
für das Rechnen von Hand äußerst umständlich, so daß man meist zu einem 
Iterationsverfahren greift (vgl. Programm 31). Sie liefert jedoch im Gegensatz 
zu den Iterationsverfahren sämtliche Lösungen der Gleichung. 


Dividiert man die kubische Gleichung durch a und führt die Substitution 


y = 


x + 


_b_ 

3a 


durch, so erhält man die sogenannte reduzierte Gleichung 
y + 3py + 2q = 0 


2b3 bc d 

27a3 3a2 a 


3ac — b2 
3a2 


Aus p und q berechnet sich die kubische Resultante 
D = q2 + p3 


Ihr Vorzeichen gibt Auskunft über die Art der Lösungen: 

Ist D > 0, so gibt es 1 reelle, 2 konjugiert komplexe Lösungen 

Ist D < 0, so gibt es 3 reelle Lösungen 

Ist D = 0, so gibt es mehrfache, reelle Wurzeln 

Mit Hilfe der 3. Wurzel 

TTir 

können nun die gesuchten Lösungen berechnet werden. Im Fall D<0 aller¬ 
dings, ist in der Cardano-Formel eine komplexe Wurzel zu ziehen. Deswegen 
wird hier eine trigonometrische Lösung durchgeführt. Die hierbei benötigte 
Umkehrfunktion arccos(x) der Cosinus-Funktion ist nicht in BASIC imple¬ 
mentiert. Sie muß daher über die Arctangens-Funktion ausgedrückt werden. 
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Es gilt 


arccos x 


\/"l — x2 

arctan- 

x 

Vl — x2 

arctan-+ n 

x 

7r 

2 


für 0 < x < 1 


für —1 < x < 0 


für x = 0 


Zum folgenden Programm 

Das Programm berechnet die Lösung der kubischen Gleichung nach dem an¬ 
gegebenen Verfahren von Cardano. 

Eine berühmte kubische Gleichung ist 
x3 + 2x2 + I0x — 20 = 0 

Sie wurde Leonardo von Pisa (ca. 1170 — 1250), genannt Fibonacci, von 
Johann von Palermo, einem Höfling von Friedrich II, zur Lösung gestellt. 
Fibonacci fand die Lösung 

x = 1.3688081075 
die auf 8(!) Dezimalen genau ist. 

Die Eingabe der Koeffizienten ins Programm liefert, neben der genannten Lö¬ 
sung, noch die beiden komplexen Lösungen 

-1.68440405 ± 3.43133135i (gerundet) 
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100 REM Kubische Gleichung 
110 : 

120 CLS:PRINT STRING*<25,42) 

130 F'RINT" Kubische Gleichung" 

140 PRINT STRING*<25,42) 

150 : 

160 PRINT"a*x•"•3+b** ■'"■2+c*x +d=0" 

170 PRINT"Eingabe der Koeffizienten a,b,c,d"; 

180 INPUT A,B,C,D 

190 IF A=0 TUEN PRINT" 

Kei ne kubisehe Gl eichung ":END 
200 A1=A:A=B/A1:B=C/A1:C=D/A1 
210 : 

220 'Parameter der reduzierten Gleichung 

230 P=A2/9+B/3 

240 Q : =A3 / 2 7 - A * B / 6+C / 2 

250 : 

260 'Kubische Resultante 
270 R=Q-'-2+F--3 
280 : 

290 'Fal1 unterscheidüng 
300 PRINT:PRINT"Loesung:" 

310 0N SGN <R)+2 GOTO 490,430,330 
320 : 

330 '1 reelle,2 komplexen Null stellen 
340 U=SQR (R) -Q: U=SGN < U) *ABS < U) -- < 1 / 3) 

350 V=~-SQR <R) -Q: V=S6N < 0) *ABS < V) < 1 /3) 

360 X 1 =4J+V—A/3 

370 'Real- und Imaginaertei1 

380 X2=-(U + V)/2-A/3:X3= <U-V>/2*SQR <3) 

390 PRINT":< 1 = " ; XI 
400 PRINT"X2=";X2;"+i*";X3 
410 PRINT"X2=";X2;"-i*";X3:END 
420 : 

430 'reelle Doppelloesung 
440 U=SGN <-Q) *ABS <Q) (1/3) 

450 Xl=LI*2--A/3: X2=-U-A/3 
460 PRINT"xl=";X1 
470 PRINT"X2=x3=";X2:END 
480 : 

490 '3 reelle Loesungen 

500 X =-Q/SQR (-F---3) 
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510 : 

520 'Definition des Arccosinus 
530 DEF FNARCC0S(X)=ATN(SQR(l~X"-2)/X) 
540 0N SGN(X)+2 GOTO 550,560,570 
550 PHI=FNARCC0S(X)+PI:GOTO 580 
560 PHI-FNARCCOS(X):GOTO 580 
570 PH I=P1/2 

580 X1=2*SQR <-P> *COS(PHI/3)-A/3 

590 X2=~2*SQR <-P)*COS<(PHl+PI)/3)-A/3 

600 X3=~2*SQR(-P)*C0S((PHI-PI)/3)~A/3 

610 PRINT"x1=";XI 

620 PRINT"x2=";X2 

630 PR I NT " 3= " ; X 3 

640 END 


Kubische Gleichung 
**■#•**•*■**•#*■*•*■*•*•*•*•*■*•*■*•*•*■*•*•* 
a * x 3+b * 2+c * + d=0 

Eingabe der Koeffizienten a,b,c,d? 1,2,10,- 

Loesung: 
xl= 1.368808 

X2=-l.684404 +i* 3.431331 
X2=-1.684404 -i* 3.431331 
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17. GLEICHUNG VIERTEN GRADES 


Ebenso wie die kubische Gleichung kann auch die Gleichung 4. Grades 
ax^ + bx3 + cx2 + dx + e = 0 

mit Hilfe von Wurzeltermen gelöst werden. Nach Ferrari kann die Lösung der 
Gleichung 4. Grades auf die kubische Gleichung 

8y3 _ 4cy2 + (2bd-8e)y + e(4c—b2) — d2 = 0 
zurückgeführt werden. Setzt man eine Lösung y in die quadratische Gleichung 

x + — (b+A)x + (y+ ———) = 0 
2 A 

mit 

A = ± \/8y + b2 - 4c 

ein, so liefert deren Lösung schließlich die gesuchten Lösungen der ursprüng¬ 
lichen Gleichung 4. Grades. 

Es ist klar, daß die Lösung nach Ferrari für das Rechnen von Hand zu kompli¬ 
ziert ist. Daher werden Gleichungen 4. Grades meist mit Iterationsverfahren 
angegegangen. Da jedoch Gleichungen 4. Grades nicht notwendig reelle Lösun¬ 
gen haben, kann es hier zu Konvergenzschwierigkeiten kommen. Da aber die 
Gleichungen vom 3. und 4. Grad häufig Vorkommen — insbesondere im Zu¬ 
sammenhang mit dem charakteristischen Polynom einer Matrix — ist das vor¬ 
liegende Programm sehr nützlich. 

Nachdem Ferrari 1545 seine Formel veröffentlicht hatte, bemühten sich die 
Mathematiker — insbesondere Euler und Lagrange — über 250 Jahre lang eine 
entsprechende Formel für die Gleichung 5. Grades zu finden. Erst 1824 zeigte 
Abel, daß Gleichungen von höherem als dem 4. Grad prinzipiell nicht mit 
Radikalen, d.h. mit Wurzeltermen zu lösen sind. Hier muß man stets zu Itera¬ 
tionsverfahren greifen. 


Zum folgenden Programm 

Das Programm löst die allgemeine Gleichung 4. Grades nach dem Verfahren 
von Ferrari. Die dabei auftretende kubische Gleichung wird nach der Cardano- 
Formel von Programm 17 gelöst. 

Als Beispiel wird die Gleichung 

x4 - 24x3 + i 50 x 2 _ 200 x - 375 = 0 

behandelt; sie ergibt sich in Programm 36 als charakteristisches Polynom der 
dort angegebenen Matrix. 

Das Programm liefert die reellen Lösungen 

—1, +5, +5, 15 
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100 REM Gleichung vierten Grades 
110 : 

120 CLB:PRINT STRING#<30,42) 

130 PRINT" Gleichung vierten Grades" 

140 PRINT STRING*<30,42) 

150 : 

160 PRINT"a*x A 4+b*x ""3+c*;■; "•2+d*x +e = 0" 

170 INPUT "Koeffizienten a,b,c,d,e";A,B,C,D,E 
ISO IF A=0 TUEN PRINT" 

Keine Gleichung 4.Grades":END 
190 PRINT:PRINT"Loesung:" 

200 : 

210 IF B=0 AND D=0 THEN 830: 'Biquadr.Gl eich. 
220 GOSUE^ 340: 'kubische Gleichung 

230 : 

240 quadratische Gleichung 
250 J=8*Y+B"'"2-4*C 

260 IF ABS <JX. 00000005 THEN J=0 
270 H1=SQR < J) :G=H1:H=<B*H1+J)/2 
280 I=Y*H1+B*Y-D:G0SUB 670 
290 : 

300 Hl=—Hl:H=<B*H1+J>/2:G=H1 
310 I=Y*H1+B*Y~D:GOSUB 670 
320 END 
330 : 

340 'kubische Gleichung 

350 Al=~C/2: Bl=B*D/4-E: Cl= <E* (4*C-B'"'2) -D""’2) /8 

360 p=-Al"-2/9+Bl/3 

370 Q=A 1 ••••3/27-A1 *B 1 /6+C1 /2 

380 'Kubische Resultante 

390 R=Q"-2+P--3 

400 : 

410 ON SGN <R)+2 GOTO 550,490,430 
420 : 

430 '1 reelle,2 komplexen Null stellen 
440 U=SQR < R) ~Q: U=SGN (U) * ABS < U) < 1 / 3) 

450 V=—SQR < R) -Q: V=SGN (V) *ABS (V) -- < 1 /3) 

460 Y=U+V—A1/3 
470 RETURN 
480 : 
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490 'reelle Doppel 1oesung 
500 U=SGN (-Q) * ABS (Q) (1 /3) 

510 : 

520 Y=2*U~Al/3 
530 RETURN 
540 : 

550 '3 reelle Loesungen 
570 X=-Q/SQR (-F--3) 

580 "Definition des Arccosinus 

590 DER FNARCCOS ( X ) — ATN (SQR (1 — X- v 2) / X ) 

600 0N SGN < X)+2 GOTO 610,630,620 
610 PHI—FNARCCOS(X)+PI:GOTO 640 
620 PHI=FNARCC0S(X):GOTO 640 
630 PH I=P1/2 

640 Y=2*SQR<-P)*C0S(PH 1/3)-A1/3 
650 RETURN 
660 : 

670 quadratische Gleichung 
680 D1=H""2-4*G*I 

690 IF G—0 AND 14=0 AND 1=0 TFIEN PRINT" 
Vierfache L.oesung x = " ; E"" (1 /4) : END 
700 IF D1<0 THEN 770 
710 : 

720 reelle Loesung 

730 PR INT (--H+SQR (Dl) ) / (2*G ) 

740 PRINT (-H-SQR(Dl))/(2*G) 

750 RETURN 
760 : 

770 komplexe Loesung 
780 D1=—D1 

790 PRINT -H/(2*G);" + i; SQR(Dl>/<2*G> 

800 PRINT —H/(2*G);i;SQR(Dl)/(2*G) 

810 RETURN 

820 : 

830 'biquadratisehe Gleichung 
840 K=C""2-4*A*E 
850 IF K<0 THEN 960 
860 Yl=(-C+SQR(C'-2-4*A*E> )/(2*A> 

870 Y2= (-C-SQR (C""2—4*A*E) ) / (2*A) 

880 IF Y1<0 THEN 900 

890 PRINT SQR(Y1):PRINT-SQR(Y1):GOTO 910 
900 PRINT"I*";SQR(-Y1) :PRINT"-1 *";SQR(-Y1) 
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910 IF Y2< 0 THEN 930 

920 PR INT SQR(Y2):PRINT-SQR(Y2):GOTO 940 
930 PRINT"I*";SQR(-Y2):PRINT"-I*";SQR(-Y2) 
940 END 
950 : 

960 RE=-C/(2*A):IM=SQR<-K)/<2*A> 

970 BT=SQR (RE--2+IM--2) 

980 R1—SQR((BT+RE)/2):I1-SQR((BT-RE)/2) 

990 BE=SQR (R1 -2+1 1 "-2) 

1000 R2=SQR((BE+Rl)/2 ):I2=SQR((BE-Rl)/2 ) 

1010 PRINT R2; " — I*"; 12 

1020 PRINT R2;"+I*";I2 

1030 PR INT-R2; " -1 * " ; 12 

1040 PRINT-R2; " +1 * " ; 12 

1050 END 


■#■*■*■** **■*■*•*■*■•*•** ■**•*•*■#•*■*•*■*** •**•*■■** 

Gleichung vierten Grades 
■***•*•#•**•*•■** ******************** 
a * 4 + b * ;■: 3+c * 2+d* ■+ e = 0 

Koef f i z i enten a,b,c,d,e? 1,-24,150,-200,-37 
Loesung: 

5 +i* 1.9301E—04 
5 ~i* 1.9301E—04 
-1 
15 
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18. POLYNOMBERECHNUNG 


Das folgende Programm bestimmt zu vorgegebenen Nullstellen das zugehörige 
Polynom. 

Die Koeffizienten a; (1 < i < n) eines normierten Polynoms n-ten Grades 

p(x) = x n + a n _i x n ~1 + a n _2 x n ~2 + a n _i x n— 1 + . .. + a 0 

sind nach dem Satz von Vieta elementarsymmetrische Funktionen der Null¬ 
stellen xj: 

xl +X 2 +X 3 +-+x n = -a n _i 

xix 2 + xix3 + xix4 + . . . + x n _i x n = a n -2 

X1X2X3 + X1X2X4 + xix 2 X 5 + . .. + x n _2 x n -1 x n = -a n -3 


xi • X 2 • X 3 • X 4 ■-• x n = (-l)n . a 0 


Zum folgenden Programm 

Das Programm berechnet die Koeffizienten des gesuchten Polynoms nach 
oben genannten Formeln iterativ innerhalb zweier verschachtelter Schleifen 

220 A(0) = 0 : A(N) = 1 

230 FOR l=N TO 1 STEP -1 

240 FOR J=H TO 1 STEP-1 

250 A(N-J) = A(N-J) - A(N—J+1) * X(N—1+1) 

260 NEXTJ 
270 NEXT I 

Als Programmbeispiel wird das Polynom mit den Nullstellen 
xi = 1, X 2 = 2, X 3 = 3, X 4 = 4, X 5 = 5 und xß = 6 
gesucht. Es ergibt sich das Polynom 

P(x) = x 6 - 21x5 + 175 X 4 _ 735 x 3 + 1624x2 - 1764x + 720. 

Das Programm funktioniert nur für reelle Nullstellen. Für komplexe Nullstellen 
müssen die entsprechenden quadratischen Polynome multipliziert werden, 
dies kann mittels Programm 20 erfolgen. 
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100 REM Polynomberechnung anhand 
110 ' gegebener Null stellen 

120 : 

130 CLS:PRINT STRING#(25,42) 

140 PRINT" Polynomberechnung" 

150 PRINT STRING#(25,42) 

160 : 

170 PRINT:INPUT "Wieviele Nul1stel1en";N 
180 DIM A(N),X(N) 

190 : 

200 FOR 1=1 TO N 

210 PRINT I;". Nul1stel1e"INPUT X(I> 
220 A(I)=0 
230 NEXT I 
240 : 

250 A(0)=0:A(N)=1 

260 FOR I=N TO 1 STEP -1 

270 FOR J=N TO 1 STEP --1 

280 A(N—J)=A(N—J)—A(N-J+l)*X(N-I+l) 

290 NEXT J 

300 NEXT I:PRINT 

310 : 

320 PRINT "Koef-f i 2 i enten des Polynoms:" 
330 FOR I=N TO 0 STEP -1 
340 PRINT"a(";I;")=";A(I) 

350 NEXT I 
360 END 
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************************* 
Polynomberechnung 
************************* 

Wieviele Null stellen? 6 

1 . Null stelle? 1 

2 . Null stelle? 2 

3 . Null stelle? 3 

4 . Null stelle? 4 

5 . Null stelle? 5 

6 . Null stelle? 6 

Koeffizienten des Polynoms 
a < 6 ) = 1 
a ( 5 ) =--21 
a < 4 )= 175 
a ( 3 ) ==-735 
a( 2 )= 1624 
a ( 1 ) =--1764 
a ( 0 >=■■ 720 
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19. KOMPLEXES HORNERSCHEMA 


Mit Hilfe des Hornerschemas kann der Wert des Polynoms 
f(x) = a n x n + a n _i x"-1 + a n —2 x n_2 + . .. + a Q . 

an der Stelle x ohne Potenzieren ausgewertet werden. Dies kann mit folgen¬ 
dem Programmstück geschehen 

100 F = A(N) 

110 FOR l=N—1 TO 0 STEP -1 
120 F=F * X +A(l) 

130 NEXT I 

Dabei ist N der Polynomgrad und F der gesuchte Funktionswert. Das Horner¬ 
schema ist auch im Betriebssystem vieler Rechner enthalten; mit ihm werden 
alle transzendenten Funktionen wie sin, arctan usw. über eine Polynomappro¬ 
ximation berechnet. Das angegebene Hornerschema kann auch für komple¬ 
xe Zahlen angewendet werden. Dazu muß die komplexe Addition und Multipli¬ 
kation geeignet sein. Sind die komplexen Zahlen als Zahlenpaare 

(RE(1), IM( 1)) (RE(2), IM(2)) 

gegeben, so kann die Addition wie folgt definiert werden 

500 REM KOMPLEXE ADDITION 
510 RE = RE(1) + RE(2) 

520 IM = IM(1) + IM(2) 

Entsprechend die Multiplikation 

600 REM KOMPLEXE MULTIPLIKATION 
610 RE = REO) * RE(2) — IM(1) * IM(2) 

620 IM = REO) * IM(2) + RE(2) * IM(1) 

Diese Rechenoperationen werden als Unterprogramm definiert und bei Bedarf 
angesprungen. Da jedoch in BASIC alle Variablen global sind, muß zuvor eine 
geeignete Parameterübergabe durchgeführt werden. 

Da hier die komplexen Zahlen als Zahlenpaare behandelt werden, müssen ent¬ 
sprechend auch die Koeffizienten des Polynoms und der Argumentwert als 
Zahlenpaare eingelesen werden. 
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Zum folgenden Programm 

Als Beispiel wird der Funktionswert des komplexen Polynoms 
p(z) = (1 + 1.5i)z3 — z2 + 3z + (1 + 2i) 
an der Stelle 

z = 1 -i 

berechnet. Die Koeffizienten werden im Programm mittels DATA-Werten ein¬ 
gelesen. 

Es ergibt sich der Polynomwert 
5—4i. 
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100 REM Komplexes HornerSchema 
110 : 

120 CLS:PRINT STRING*(25,42) 

130 PRINT" Komplexes Hornerschema" 
140 PRINT STRING*(25,42) 

150 : 

160 READ N:'Polynomgrad 
170 PRINT"Polynomgrad=";N 
180 DIM R(N),I(N),RE(2),IM(2) 

190 : 

200 PRINT"Koe-f -F i z i enten: " 

210 F0R K=N T0 0 STEP -1 


220 

READ 

RE, IM: ' 

Real- Imag.teil 

230 

GOSUB 470 


240 

R (K) 

=RE:I(K) 

= IM 

250 

NEXT K 



260 

PRINT" 

Argument 

— II „ 

!• 

270 

READ 

RE, IM 



280 GOSUB 470 
290 X=RE:Y=IM 
300 : 

310 REM Hornerschema 
320 RE <1)=R(N):IM(1)=I (N) 

330 FÜR K=N—1 TO 0 STEP -1 
340 RE(2)=X:IM(2)=Y 
350 GOSUB 560:'kompl.Mult. 
360 RE(1)=RE:IM(1)=IM 
370 RE(2)=R(K):IM(2)=I(K) 

380 GOSUB 510:'kompl.Add. 

390 RE(1)=RE:IM(1)=IM 
400 NEXT K 
410 : 

420 PRINT:PRINT"Polynomwert="; 
430 GOSUB 470 
440 END 
450 : 

460 REM Ausgabeprozedur 
470 IF IM>=0 THEN 

PRINT RE; " +i •* " ; IM: GOTO 490 
480 PRINT RE;"-i*";ABS(IM) 

490 RETURN 
500 : 
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510 REM Komplexe Addition 
520 RE=RE(1)+RE(2) 

530 JM=IM(1)+ IM(2) 

540 RETURN 
550 : 

560 REM Komplexe Multiplikation 
570 RE=RE(1)*RE(2)-IM<l)*IM<2> 
580 IM=RE<1)* IM(2)+1M(1)*RE(2) 
590 RETURN 
600 : 

610 DATA 3 

620 DATA 1,1.5,-1,0,3,0,1,2 
630 DATA 1,-1 


****** ****************** * 
Komplexes Hornerschema 
************************* 
Polynomgrad= 3 
Koe-f f i z i enten: 

1 +i* 1.5 
-1 +i* 0 
3 + i * 0 
1 + i * 2 

Argument= 1 -i* 1 
Polynomwert= 5 -i* 4 
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20. POLYNOM-MULTIPLIKATION 


Zwei Polynome vom Grad n und m 

a n x n + a n _ix n- 1 + a n _ 2 xn ~2 + . . . + a 0 

b m x m + b m _ix m - 1 + b m _ 2 x m— 2 + . . . + b 0 

werden multipliziert, indem man jeden Term des ersten Polynoms mit jedem 
des zweiten multipliziert. 

Dies kann mit Hilfe zweier verschachtelter Schleifen geschehen 

290 FOR 1=0 TO N 
300 FOR J=0 TO M 
310 K = I +J 

320 C(K) = C(K) + A(l) * B(J) 

330 NEXTJ 
340 NEXT I 

Das Verfahren läßt sich direkt auch auf Reihen übertragen; d.h. es liefert auch 
das Cauchy-Produkt zweier Reihenentwicklungen. 


Zum folgenden Programm 

Das Programm berechnet nach dem angegebenen Verfahren das Produkt 
zweier Polynome und gibt die entsprechenden Koeffizienten aus. 

Als Beispiel werden die Polynome 

p(x) = x^ — 5x3 + 3 x 2 + 6x — 2 - a 0 

q(x) = x2 — 1 

multipliziert. Eingabe der Koeffizienten von links nach rechts ins Programm 
liefert die Koeffizienten des Produktpolynoms 

1, —5,2, 11,-5,-6,2 

Das Ergebnis ist somit 

x 6 _ 5 x 5 + 2x4 + 11 x 3 — 5 x 2 — 6 x + 2 
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»Darstellung der Multiplikation 
3069-45 = 138105 
aus einer chinesischen Arithmetik 
des Jahres 1355«. 
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100 REM Polynom-Multiplikation 
110 : 

120 CLS:PRINT STRING*(25,42) 

130 PRINT" Polynom-Multiplikation" 
140 PRINT STRING*(25,42) 

150 : 

160 INPUT"Grad des 1.Polynoms";N 
170 DIM A(N) 

180 FÜR I=N T0 0 STEP -1 

190 PRINT"Koeffizient x<" ; I;")="; 

200 INPUT A(I) 

210 NEXT I:PRINT 

220 : 

230 INF'UT"Grad des 2. Pol ynoms" ; M 
240 DIM B(M),C(N+M) 

250 F0R I=M T0 0 STEP -1 

260 PRINT"Koeffizient x ( " ; I ; " )="; 

270 INF'UT B(I) 

280 NEXT I:PRINT 
290 : 

300 F0R 1=1 TO N+M 
310 C(I)=0 
320 NEXT I 
330 FOR 1=0 TO N 
340 FOR J=0 TO M 
350 K=I+J 

360 C(K)=C(K)+A(I)*B( 3 ) 

370 NEXT J 
380 NEXT I 
390 : 

400 PRINT"Koe-f f i z i enten d. Produkts:" 
410 FOR I=N+M TO 0 STEP -1 
420 PRINT Cd); 

430 NEXT I:PRINT 
440 END 
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************************* 
Pol ynom-liul tipl i kati on 
************************* 


Grad des l.Polynoms? 

4 

Koef fizient x( 

4 

) =? 

1 

Koef fizient x ( 

“7 

) =? 

—5 

Koeffizient x ( 

•~y 

) =? 


Koe-f f i z i ent x ( 

l 

) =? 

6 

Koeffizient x< 

0 

) =? 

_ •-} 

Grad des 2.Polynoms? 

n 

X. 

Koeffizient x( 


) =? 

l 

Koeffizient x( 

1 

) =? 

0 

Koeffizient x( 

0 

) =? 

-1 


Koef f i 2 i enten d . Produkts: 
1-5 2 11-5-6 2 
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21. POLYNOM-DIVISION 


Die Polynom-Division findet vielfache Anwendung beim Abdividieren von ge¬ 
fundenen Nullstellen, bei der Partialbruchzerlegung u.s.w. 

Die Polynom-Division kann in völliger Analogie zur gewöhnlichen Division 
durchgeführt werden: 

Es wird geprüft, wie oft der entsprechende Koeffizient des Divisors in den 
Dividenden hineingeht. Das entsprechende Vielfache des Divisors wird dann 
vom Koeffizienten des Dividenden subtrahiert. 


Zum folgenden Programm 

Das Programm dividiert zwei Polynome nach dem gewöhnlichen Divisionsver¬ 
fahren. 

Es wurde versucht Rundungsfehler aufzufangen, indem sehr kleine Reste Null 
gesetzt werden. Die Polynomgrade und Koeffizienten werden in Form von 
DATA-Werten eingelesen. 

Als Beispiel werden die Polynome 

2 x 7 + g x 6 _ 5x5 + 42x4 _ 24x3 + io x 2 + 21x - 8 
x3 + 5x2 — 3x + 4 

dividiert. Das Programm liefert den Quotienten 

2x4 _ x 3 + 6x2 + x - 7 

mit dem Rest 

24x2 - 4x + 20 
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100 REM F'olvnom-Di Vision 
110 : 

120 CLS:PRINT STRING*(25,42) 

130 PRINT" Polynom-Division" 

140 PRINT STRING*(25,42) 

150 : 

160 PRINT"Koe-f f i z i enten des Dividenden 

170 READ N:'rem Polynomgrad 

180 DIM ACN+l),B(N+1),Q(N+1>,R(N+1) 

190 FÜR I=N+1 TO 1 STEP -1 
200 READ A(I):PRINT A(I); 

210 NEXT I:PRINT:PRINT 
220 : 

230 PR INT “Koe-f -f i z i enten des Divisors" 
240 READ M 

250 FÜR I=M+1 TO 1 STEP -1 
260 READ B(I):PRINT B(I); 

270 NEXT I:PRINT:PRINT 
280 : 

290 IF N >—M TUEN 320 

300 F‘R I NT " Pol ynomgrad des Dividenden" 
310 F’R INT "< = der des Di vi sors" : END 
320 R=N-M 

330 FÜR I—N+1 TO M+l STEP -1 
340 Q <I-M)=A(I)/B <M+1) 

350 FÜR J=0 TO M 

360 A(I-J)=A(I-J)-Q(I-M)*B(M+1-J) 

370 IF ABS(A(I-J))>=ABS(Q(I-M)* 

B(M+1-J))*.0000001 THEN 390 
380 A(I-J)=0 

390 NEXT J 
400 NEXT I 
410 : 

420 F’RINT"Koe-f -f i z i enten d. Quotienten: " 
430 FOR I=R+1 TO 1 STEP -1 
440 PRINT Q(I); 

450 NEXT I:PRINT:PRINT 

460 : 

470 PR INT "Koe-f -f i 2 i enten d . Restpol ynoms 
480 IF A(R)=0 THEN R=R-1:G0T0 480 
490 FOR I=R TO 1 STEP -1 
500 R(I)=A(I) 
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510 

PRINT RU): 

520 

NEXT 

I:PRINT 

530 

END 


540 

: 


550 

DATA 

7 

560 

DATA 

'" 7 * Q _ R 40 

1 ' 1 u f l 

570 

DATA 

~r t 

580 

DATA 

1,5,-3,4 


************************* 

Pol ynom~Division 
************************* 
Koeffizienten des Dividenden 
2 9 -5 42 -24 10 21 -8 

Koeffizienten des Divisors 
1 5-3 4 

Koef f i z i en t. en d. Quotienten: 
2-1 6 1-7 

Koeffizienten d, Restpolynoms 
24 -4 20 
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22. MATRIZEIM-MULTIPLIKATION 


Matrizen sind rechteckige Anordnungen von reellen Zahlen. Eine Matrix hat 
die Ordnung (m, n), wenn sie m Zeilen und n Spalten besitzt. 

Zwei Matrizen A und B werden multipliziert, indem man jede Zeile von A mit 
jeder Spalte von B multipliziert. Dies zeigt, daß das Produkt von A und B nur 
existiert, wenn die Spaltenzahl von A mit der Zeilenzahl von B übereinstimmt. 
Ist A von der Ordnung (n, m) und B der Ordnung (p, q) mit m=p, so ist die 
Produktmatrix von der Ordnung (n, q). 

Die Matrizenmultiplikation kann mit Hilfe von folgendem Programmstück 
durchgeführt werden 

320 FOR 1=1 TO N 
330 FOR J=1 TO Q 
340 S=0 

350 FOR K=1 TO M 
360 S=S +A(I,K) * B(K,J) 

370 NEXT K 
380 C(I,J) = S 
390 NEXT J 
400 NEXT I 


Zum folgenden Programm 

Das folgende Programm berechnet das Produkt zweier Matrizen, die in Form 
von DATA-Werten eingelesen werden. Dabei wird geprüft, ob die Spaltenzahl 
der ersten Matrix mit der Zeilenzahl der zweiten übereinstimmt, andernfalls 
wird eine Fehlermeldung ausgegeben. 

Im Programmbeispiel werden die Matrizen 



miteinander multipliziert. Die Produktmatrix ist 



(vgl. Programmausdruck). 
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100 REM Matrizenmultiplikation 
110 : 

120 CLS:PRINT STRING*(25,42) 

130 PRINT" Matr i zenmul t.i p 1 i kati on " 
140 PRINT STRING*(25,42) 

150 : 

160 PRINT"Gegebene Matrizen:“ 

170 READ N,M:'Ordnung der 1.Matrix 
ISO DIM A(N,M) 

190 FQR 1=1 TO N 
200 FOR J=1 TO M 

210 READ A(I,J):PRINT A(I,J); 

220 NEXT J:PRINT 
230 NEXT I:PRINT 
240 : 

250 READ P,Q:'Ordnung der 2.Matrix 
260 DIM B(P,Q) 

270 IF M=P TUEN 320 

280 PRINT"Eingabefehler:Spaltenzahl" 
290 PRINT"der 1.Matrix muß mit der" 
300 PRINT"Zeilenzahl der 2.Matrix" 
310 PRINT"übereinstimmen":END 
320 FOR 1=1 TO P 
330 FOR J=1 TO Q 

340 READ B (I , J > :PR INT B(I,J); 

350 NEXT J:PRINT 
360 NEXT I:PRINT 

370 : 

380 PRINT“Produktmatrix:" 

390 FOR 1=1 TO N 
400 FOR J=1 TO Q 
410 S=0 

420 FOR K=1 TO M 

430 S=S+A(I,K)*B(K,J) 

440 NEXT K 

450 C(I,J)=S:PRINT S; 

460 NEXT J:PRINT 
470 NEXT I:PRINT 
480 END 
490 : 

500 DATA 3,4 
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510 

DATA 

8 , 5 , 2,3 

520 

DATA 

6 , 4 , 1 , 0 

530 

DATA 

1 , 2 , 2,4 

540 

: 


550 

DATA 

4,2 

560 

DATA 

1 —r> 

57 O 

DATA 

_ r ~i 

580 

DATA 

1,1 

590 

DATA 

1,2 


* ******** ********** ****** 
Matrizenmultiplikation 
********* **************** 


Gegebene Matrizen: 
8 5 2 3 

6 4 1 0 

12 2 4 


1 -2 
t 

i i 

1 2 


Produktmatrix: 

3 7 

-1 1 
3 14 
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23. ZWEIERPOTENZEN VON MATRIZEN 


Da für quadratische Matrizen die Zeilenzahl mit der Spaltenzahl überein¬ 
stimmt, können quadratische Matrizen mit sich selbst multipliziert werden. 

Solche Potenzen von Matrizen haben zahlreiche Anwendungen in der numeri¬ 
schen und Wirtschaftsmathematik, Wahrscheinlichkeitsrechnung und Graphen¬ 
theorie. 

Als Beispiel werden hier die (homogenen) Markowketten herausgegriffen: 

Hat ein System endlich viele Zustände und geht ein Zustand i, unabhängig 
von seinem vorhergehenden Zustand, in den Zustand j mit der Wahrschein¬ 
lichkeit 


Pij 

über, so stellt das System eine Markowkette (auch Markoff geschrieben) dar. 
Die Matrix P mit den Elementen 

Pij 

heißt eine stochastische Matrix. Solche Matrizen sind dadurch gekennzeichnet, 
daß ihre Zeilensumme den Wert 1 haben, da die Summe aller Übergangswahr¬ 
scheinlichkeiten natürlich den Wert 1 hat. 

Hat dieses System ein stationäres Verhalten, d.h. strebt es einer Grenzvertei¬ 
lung zu, so kann diese Grenzverteilung durch fortgesetztes Potentieren berech¬ 
net werden [10]. 

Multipliziert man die Matrix mit sich selbst und quadriert das jeweilige Ergeb¬ 
nis, so erhält man die Zweierpotenzen der Matrix. 


Zum folgenden Programm 


In einem bestimmten Land gebe es 3 Parteien A, B und C. Folgendes Wahlver¬ 
halten sei gegeben: 


80% aller 
10% C. 

A-Wähler 

wählen 

bei 

der 

nächsten 

Wahl 

wieder 

A, 

10% 

B 

und 

70% aller 

15% C. 

B-Wähler 

wählen 

bei 

der 

nächsten 

Wahl 

wieder 

B, 

15% 

A 

und 

60% aller 
20% B. 

C-Wähler 

wählen 

bei 

der 

nächsten 

Wahl 

wieder 

c. 

20% 

A 

und 
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Damit ergibt sich folgende stochastische Matrix 


0.8 

0.1 

0.1 

0.15 

0.7 

0.15 

0.2 

0.2 

0.6 


Wie man dem Programmausdruck entnimmt, nähern sich die Zeilen der Matri¬ 
zenpotenzen dem Vektor 

(0.461912194 0.307419108 0.230668699) 

Dieser stellt die Grenzverteilung dar; d.h. auf lange Sicht wählen 

46.2 % aller Wähler A 
30.7 % B 

23.1 % C 
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100 REM Zweierpotenz einer Matrix 
110 : 

120 CL5:PRINT STRING*(30,42) 

130 F'RINT“ Zweierpotenz einer Matrix" 
140 PRINT STRING*(30,42) 

150 : 

160 READ N:'Ordnung der Matrix 

170 READ P:'gewuenschte Zweierpotenz 

180 DIM A(N,N),B(N,N) 

190 : 

200 PR INT"Gegebene Matrix:" 

210 FOR 1=1 TO N 
220 FOR J=1 TO N 
230 READ A(I,J):PRINT A(I,J>; 

240 NEXT J:PRINT 
250 NEXT I:PRINT 
260 : 

270 P=INT(LOG(P)/LOG(2)+.5) 

280 FOR L=1 TO P 

290 PRINT 2""L; " . Potenz:" 

300 FOR 1=1 TO N 

310 FOR J=1 TO N 

320 S=0 

330 FOR K=1 TO N 

340 S=S+A(I,K)*A(K,J) 

350 NEXT K 

360 B(I,J)=S:PRINT S; 

370 NEXT J:PRINT 

380 NEXT I:PRINT 

390 : 

400 FOR 1=1 TO N 

410 FOR J=1 TO N 

420 A <I,J)=B(I,J) 

430 NEXT J 

440 NEXT I 

450 NEXT L 
460 END 
470 : 

480 DATA 3 

490 DATA 16 

500 DATA .8,. 1 , . 1 
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cn ui 


10 DATA .15,.7,.15 
20 DATA .2,.2,.6 


****************************** 
Zweierpotenz einer Matrix 
****************************** 
Gegebene Matrix: 

0.8 0.1 0.1 

0.15 0.7 0.15 

0.2 0.2 0 . 6 

2 . Potenz: 

0.675 0.17 
0.255 0.535 
0.31 0.28 

4 . Potenz: 

0.547025 
0.37365 
0.40775 


0. 155 

0.21 

0.41 


203875 

237975 

27495 


0.2491 0. 
0.388375 0. 
0.3173 0. 


8 . Potenz: 
0.475442597 
0.446546517 
0.453719452 

16 . Potenz: 
0.461912195 
0.461128662 
0.461337398 


0.297697678 
0.319420823 
0.312043548 


0.307419108 
0.307993905 
0.307836579 


0.226859726 
0.234032661 
0.234237002 


0.230668699 

0.230877434 

0.230826025 
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24. MITTELWERTE 


Neben dem arithmetischen Mittel A der Zahlen xj (1 < i < n) 
A = 1 (xi + X2 + X 3 + X 4 . . . + x n ) 

sind auch andere Mittelwerte von Bedeutung: 

das geometrische Mittel 

_ n/- 

G = V xi • X 2 • X 3 ■ X 4 •-■ x n 

das harmonische Mittel 

- = n( —+ — + — + —+ ) 

H xi X 2 X 3 X 4 x n 


das quadratische Mittel 


Q 


= J\ 


1 , 2 , 2 2 , 2. 
(x 1 + x 2 + x 3 + x 4 + .... +x n ) 


Das meistgenutzte arithmetische Mittel ist sehr empfindlich gegen Ausreißer; 
d.h. bei kleinen Zahlenreihen wird das arithmetische Mittel durch einen zufäl¬ 
lig großen Wert verfälscht. In diesem Fall sollte man lieber den Median (vgl. 
Programm 25) nehmen. 

Das geometrische Mittel ist stets bei prozentualen Größen zu nehmen. Wird 
z.B. eine Gehaltserhöhung von 3% erst zum 1. Juni ausgezahlt, so ist die effek¬ 
tive Gehaltserhöhung nur 

( ^ 1.037 _ i) . 100 %= 1.7 % 

Hier findet man häufig falsche Angaben, auch in angesehenen Tageszeitungen. 

Fährt ein Auto eine Strecke bergauf mit 30 km/h und bergab mit 90 km/h, so 
ist die Durchschnittgeschwindigkeit keineswegs 60 km/h, sondern nur das har¬ 
monische Mittel von 45 km/h. 


Für das arithmetische, geometrische und harmonische Mittel gilt die Unglei¬ 
chung von Cauchy 

A > G > H 
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Das quadratische Mittel hat große Bedeutung in der Statistik, da die Standard¬ 
abweichung einer Meßreihe gleich dem quadratischen Mittel der einzelnen Ab¬ 
weichungen ist. 


Zum folgenden Programm 

Das Programm berechnet für beliebige Zahlenreihen das arithmetische, geome¬ 
trische, harmonische und quadratische Mittel nach den oben angegebenen 
Formeln. 

Als Beispiel wurde die Zahlenreihe: 

19,22,20, 18, 19,23, 17,21,25, 17 
eingegeben. Es ergeben sich folgende gerundete Mittel: 

arithmetisches M.: 20.10 

geometrisches M.: 19.95 

harmonisches M.: 19.80 

quadratisches M.: 20.26 
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100 REM Mittelwerte 
110 : 

120 CLS:PRINT STRING*(25,42) 

130 PRINT" Mittelwerte " 

140 PRINT STRING*(25,42) 

150 : 

160 INPUT "Wieviele Zahlen ";N 
170 A=0:G=1:H=0:Q=0:DF=1 
180 : 

190 PRINT:PRINT"Eingabe der Zahlen:" 
200 FOR 1=1 TO N 

210 PRINT USING" ##. Zahl";I;:INPUT X 
220 A=A+X: G=G*X : Q=Q + X"'"2 
230 IF X=0 THEN DF=0:G0T0 250 
240 H=H+1/X 
250 NEX'T I: PRINT 
260 : 

270 A=A./N 
280 G=G’- (1/N) 

290 0=SQR(Q/N) 

300 IF HOO THEN H=N/H 
310 : 

32C> F : 'R INT "Ar i thmet i sch es Mittel™"; 

330 PRINT USING "#####.###";A 
340 PRINT" Geometrisches Mittel ™"5 
350 PRINT USING "#####.###";G 
360 IF DF=0 THEN 370 ELSE 380 
370 PRINT" Harmonisches Mittel 

nicht definiert":GOTO 400 
380 PRINT" Harmonisches Mittel™"; 

390 PRINT USING "#####.###";H 
400 PRINT " Quadratisches Mittel™"; 
410 PRINT USING "#####.###";Q 
420 END 
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************************* 
Mittel werte 

************************* 
Wieviele Zahlen ? 10 

Eingabe der Zahlen: 


1 . 

Zahl? 

19 

r ~y 

Zahl? 

'•? 


Zahl? 

20 

4. 

Zahl? 

18 

5 . 

Zahl? 

19 

6 . 

Zahl? 

r-yrr 

“7 

/ m 

Zahl? 

17 

8 . 

Zahl? 

21 

9. 

Zahl? 

25 

10. 

Zahl? 

17 


Arithmetisches 
Geometrisches 
Harmonisches 
Quadratisches 


Mittel= 

20.100 

M i 11: e 1 = 

19.948 

Mittel= 

19.801 

Mittel= 

20.256 
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25. STATISTISCHE MITTELWERTE 


Neben dem arithmetischen Mittel werden in der Statistik auch noch andere 
Mittelwerte betrachtet. 

Der Modus ist der häufigste Wert einer Stichprobe oder einer Meßwertreihe. 
Für Durchschnittslohnangaben ist sicher der Modus aussagekräftiger als der 
Mittelwert, der durch Bezieher sehr großer Einkommen verfälscht wird. 

Denkt man sich alle Stichprobenwerte der Größe nach geordnet, so stellt das 
in der Mitte liegende Element den Median dar. Ist die Zahl der Werte gerade, 
so wählt man den Mittelwert der beiden mittleren Elemente zum Median; 
z.B. beim Stichprobenumfang 100 ist der Median das Mittel aus dem 50. und 
51. Stichprobenwert. Der Median hat sehr große Bedeutung beim Testen auf 
vorgegebene Verteilungen (siehe z.B. [16]). Bei einer geringen Zahl von Stich¬ 
probenwerten ist das arithmetische Mittel besonders empfindlich gegen 
"Ausreißer". So verwendet die Stiftung Warentest als mittleren Preis eines 
Gerätes den Median, wenn weniger als 5 Geräte gekauft werden. 

Weitere Parameter, die die Verteilung von Zahlen kennzeichnen, sind die 
Varianz, die Schiefe und der Exzeß. 


Sind xj (1 < = i < = n) die Stichprobenwerte, so stellt 

(£xj>2 


n—1 

die Varianz dar. Sie ist ein Maß für die Streuung der Werte um ihren Mittel¬ 
wert. Die Quadratwurzel aus der Varianz wird Standardabweichung o ge¬ 
nannt. Im Falle einer Normalverteilung gibt letztere Zahl den Bereich um 
den Mittelwert an, in dem sich 65% aller Werte befinden. 


Die Schiefe wird berechnet aus 

Zxi3 £x;2 Zx; 


+2 



/a3 


Sie ist ein Maß für die Symmetrie der Verteilung. 


Der Exzeß 
e = 




n 


Sxj2 

n 



/o4 -3 
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(auch Wölbung genannt) ist ein Maß für die Steilheit der Kurve. Positiver 
Exzeß zeigt, daß die Verteilungskurve steiler als die Normalverteilung verläuft, 
negative Werte zeigen einen flacheren Verlauf. Für die Normalverteilung ver¬ 
schwindet sowohl Schiefe wie Exzeß. 


Zum folgenden Programm 

Gibt man die Stichprobenwerte in Form von DATA-Werten ein, so werden 
die oben genannten statistischen Kennwerte nach den angegebenen Formeln 
berechnet. Existieren zwei Werte mit gleich großer Häufigkeit, so wird der 
kleinere als Modus gewählt. 

Für die Stichprobe 

(18,20,21, 19,22,23, 17,21,18,24) 

können die sich ergebenden Maßzahlen dem Programmausdruck entnommen 
werden. 

Die Schiefe ist klein, da genau 5 Werte unter und 5 Werte über dem Mittel von 
20.3 liegen. 

Der Exzeß ist negativ, dies zeigt eine flache Verteilung an. 
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100 REM Statistische Mittelwerte 
1 10 : 

120 CLS:PRINT STRINGS(30,42) 

130 PRINT" Statistische Mittelwerte 
140 PRINT STRING*<30,42) 
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640 IF (H (I) =K) AND (IOL) THEN M=1 
650 NEXT I 

660 IF K=0 THEN PRINT"Modus nicht de-f i ni er t" : GOTO 690 
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Median= 20.5000 
Modus“ 18.0000 
Modus nicht eindeutig 






"Pythagoras zeigt an Glocken, gespannten 
Saiten und Flöten den Zusammenhang von 
Harmonie und ganzen Zahlen. Aus Theorica 
Musice, Mailand 1492" 
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26. PYTHAGOREISCHE ZAHLEN 


Ganzzahlige Lösungen ¥= 0 der Diophantischen Gleichung 
x2 + y 2 = z 2 

heißen pythagoreische Zahlen. Ein Dreieck ist genau dann rechtwinklig, wenn 
seine Seitenlängen pythagoreische Zahlen sind. 

Setze man zwei solcher rechtwinkeligen Dreiecke mit gleicher Kathete zusam¬ 
men, so erhält man ein Dreieck, dessen Höhe ganzzahlig ist. Die Seitenlängen 
eines solchen Dreiecks nennt man auch Heronsche Zahlen, da sie rationale 
Werte der Heronschen Flächenformel liefern. 

Bekannte pythagoreische Tripel sind 

(3,4,5) 

( 8 , 6 , 10 ) 

(5, 12, 13) 

(9, 12, 15) 

Pythagoreische Zahlen können durch die Formeln 
x = u 2 - v2 
y = 2uv 
z = u 2 + v 2 

erzeugt werden; dabei durchlaufen u und v alle natürlichen Zahlen. Die 
obengenannten Tripel erhält man für 

u = 2, v = 1 
u = 3, v = 1 
u = 3, v = 2 
u = 4, v = 2 

Wie man an den Beispielen sieht, sind auch stets alle Vielfachen eines solchen 
Tripels wieder pythagoreisch. Im folgenden Programm werden daher nur tei¬ 
lerfremde Tripel erzeugt. 


Zum folgenden Programm 

Das folgende Programm berechnet nach dem oben angegebenen Verfahren be¬ 
liebig viele pythagoreische Zahlen, wobei Vielfache anderer Tripel ausgeschlos¬ 
sen werden. 

Die ersten 20 vom Programm erzeugten Zahlentripel können dem Programm¬ 
ausdruck entnommen werden. 
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Pythagoreische Zahlen 
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27. KOORDINATEN BESONDERER DREIECKSPUNKTE 


Sind jdI, jd2, p3 die Ortsvektoren der Dreiecks-Eckpunkte, so ergeben sich die 
Ortsvektoren der besonderen Dreieckspunkte aus folgenden Formeln: 

Schwerpunkt 

_ £1 + £2+£3 


Inkreismittelpunkt 

aifii + a2P2 + a3P3 

i = - 

ai + a2 + a3 

Höhenschnittpunkt 

h = mi +t2£2 + t3£3 
t1t2t3 


Umkreismittelpunkt 

cisigi + C2S2£2 + C 3 S 3 £3 

2 si S2S3 


dabei stellen ai,a 2 , a 3 die Seitenlangen, 
s 1. s 2. s 3 die Sinuswerte, 
c 1 , C 2 , C 3 die Cosinuswerte, 
t 1- t 2- t 3 die Tangenswerte 

des Dreiecks dar. 


Zum folgenden Programm 

Das Programm berechnet die 4 besonderen Punkte eines Dreiecks nach den 
oben angegebenen Vektorgleichungen. 

Als Programmbeispiel wird das räumliche Dreieck 

(5 101-4), (-21613), (61-31-5) 
behandelt. Es ergeben sich folgende Werte, 

Schwerpunkt: 

(3| 1 I -2) 
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Inkreismittelpunkt: 

(4.60337933 | -0.50521699 | -3.60337933) 

Höhenschnittpunkt: 

(14.6 | 17.0666654 1 -13.56) 

Umkreismittelpunkt: 

(-2.8 I -7.03333264 | 3.8) 

Zur Kontrolle kann man den Satz anwenden,daßder Schwerpunkt die Strecke 
Umkreismittelpunkt — Höhenschnittpunkt im Verhältnis 1:2 teilt 

s = - (2 • u + h) 

3 

Dies ist, wie man leicht nachrechnet, erfüllt. 
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***********■**•***•****■****■**■*■#••*•*■*•*•*•*•**••**•#•* 


Koordinaten 

, y, z 

des 

1.Punkts? 

5,0,-4 

Koordinaten 

- iYi- 

des 

2.Punkts? 

-2,6,3 

Koordinaten 

x 1 V,z 

des 

3.Punkts? 

6,-3,-5 


Schwerpunkt: 

3 1 -2 

I n k r ei smi 11 e 1p un kt: 

4.60337933 -0.505216989 

Hoehenschnittpunkt: 

14.6000006 17.066669 

Umkreismittelpunkt 
-2.79999994 -7.03333456 


-3.60337933 
-13.6000003 
3.79999987 
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28. POLYGON-BERECHNUNG 


Das folgende Programm berechnet Flächeninhalt und Umfang eines Vielecks 
bei bekannten Koordinaten der Eckpunkte. Der Flächeninhalt wird nach der 
Formel, vgl. [4] 

A = [(xi + X 2 ) (yi - y2> + (x 2 + X 3 ) (y 2 - yß) + . .. 

+ (x n + xi) (y n — yi )]/2 

berechnet. 

Da der Abstand des Punktes (xj+i I yj+ 1 ) vom Punkt (xj I yj) nach Pythagoras 
gleich 

d = 

ist, ergibt sich der Umfang des Vielecks durch Aufsummieren der Abstände 
der einzelnen Eckpunkte. 

Der Schwerpunkt ergibt sich aus 
1 n 

x s = gß 2 UiYi+l - xj+iyj) (xj + xj+i) mit x n +i = xi 

1 n 

y s = 6 Ä?! (xivi+1 — x i+ 1 Vi) (yj + Yi+l) y n +1 =V 1 


Zum folgenden Programm 

Das Programm berechnet für ein beliebiges Vieleck Umfang und Flächeninhalt. 
Bei der Eingabe der Koordinaten muß darauf geachtet werden, daß der Um¬ 
laufsinn eingehalten wird, da die Formel nicht für "überschlagene" Vielecke 
gilt. 

Als Programmbeispiel wird das Sechseck 

(412), (10 10), (1219), (8111), (3 15). (01 3) 
gewählt. Es ergibt sich 

Flächeninhalt = 67 

Umfang = 35.5551423 

Schwerpunkt = (7.35323384 I 5.20895523) 
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"Visieren von Winkeln mit Hilfe des Jakobstabs. Holzschnitt aus dem 16. Jahrhundert" 
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100 REM Polygonberechnung 
1 10 : 

120 CLS: PRINT STRING* <25,42) 

130 PRINT" Polygonberechnung" 

140 PRINT STRING*(25,42) 

150 : 

160 INPUT"Wieviele Eckpunkte";N 
170 DIM X(N+l),Y<N+1> 

180 : 

190 PRINT"6ib Koordinaten der Eckpunkte" 

200 PRINT"im gleichen Umlaufssinn ein !" 

210 FOR 1=1 T0 N 

220 PRINT"X<";I;"),Y(";I;"):INPUT X(I),Y(I) 

230 NEXT I 

240 X(N+1)=X(1)s Y(N+l> =Y <1) 

250 : 

260 ' FTaechenberechnung 
270 F=0 

280 FOR 1=1 TO N 

290 F=F+ ( X (I ) + X (I +1) ) * (Y (I) Y (I +1) ) 

300 NEXT I 
310 F=ABS(F)/2 
320 : 

330 'Umfang 
340 U=0 

350 FOR 1=1 TO N 

360 U=U+SQR ( ( X (I +1) - X (I) ) •-•2+ (Y (I +1) - Y (I) ) A 2) 

370 NEXT I 
380 : 

390 'Schwerpunktskoordinaten 
400 XS=0:YS=0 
410 FOR 1=1 TO N 

420 XS=XS+(X(I)*Y(1 + 1)-X(1 + 1)*Y<I))*(X (I)+X(1 + 1)) 
430 YS=YS+(X <I>*Y<1 + 1>-X <1 + 1> *Y <I> > * <Y <I> +Y <I + l)> 
440 NEXT I 

450 XS=XS/<6*F):YS=YS/(6*F) 

460 : 

470 PRINT: F'R I NT "Fl aechen inhalt = ";F 
480 F’RINT'TJmf ang=" ; U 

490 PRINT"Schwerpunkt (";XS; "," ; YS;")" 

500 END 
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************************* 

Pol ygonberechnung 
************************* 
Wieviele Eckpunkte? 6 
Gib Koordinaten der Eckpunkte 


i m 

g1 eichen 

Umlaufseinn ein 

X ( 

1 

) 

,Y< 

1 

) ? 

4,2 

X ( 

ü? 

) 

,Y< 

9 

) ? 

10,0 

X ( 


) 

,Y( 

y % 

> ? 

12,9 

X ( 

4 

) 

,Y( 

4 

) ? 

8,11 

X ( 

5 

) 

,Y( 

5 

) ? 

^ 

X ( 

6 

) 

,Y( 

6 

) ? 

0 , 3 


F1 aechen i nha 11-= 67 
Umfang“- 35.5551423 

Schwerpunkt < 7.35323383 , 5.20895522 ) 
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29. DETERMINANTE 


Determinanten haben zahlreiche Anwendungen: Mit ihrer Hilfe kann geprüft 
werden, ob 

— Vektoren linear abhängig sind 

— ein lineares Gleichungssystem lösbar ist 

— zwei Polynome gemeinsame Nullstellen haben 

— ein Extremum einer Funktion mehrerer Variablen vorliegt 

usw. 


Die Resultante der Polynome 

f (x) = x3 + 2x2 + 3x + 5 
g(x) = 6x2 + 8x + 9 


Det (f,g) 


1 2 3 5 0 
0 12 3 5 
6 8 9 0 0 
0 6 8 9 0 
0 0 6 8 9 


z.B. ist genau dann Null, wenn f und g mindestens eine gemeinsame Null¬ 
stelle haben. Hier gilt Det (f,g) = 2873, somit haben f und g keine gemeinsa¬ 
men Nullstellen. Ersetzt man g durch die Ableitung f'von f, so gibt die De¬ 
terminante Auskunft, ob f mehrfache Nullstellen hat. 


Erwähnenswert ist auch die Cramersche Regel, die erlaubt einzelne Unbekann¬ 
te eines Gleichungssystems zu finden, ohne das ganze System lösen zu müssen. 
Ein Programm dazu findet sich in [15]. 


Zum folgenden Programm 

Das Programm berechnet die Determinante durch fortgesetztes Entwickeln 
nach der 1. Spalte. 

Dazu werden die Zeilen 2 bis n mit dem Element an ^ 0 multipliziert. Da¬ 
durch vervielfacht sich zunächst der Wert der Determinante um 
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Der entsprechende Faktor wird sodann abdividiert. Subtrahiert man von jeder 
der Zeilen 2 bis n das 

an -fache von Zeile 1 
so wird in der 1. Spalte der Vektor 

(1 1010 I.... 10) 

erzeugt. 

Entwickelt man nun nach der 1. Spalte, so verbleibt eine Determinante der 
Ordnung n—1. 

Setzt man das Verfahren fort, so verringert sich die Ordnung der Determinan¬ 
te bei jedem Schritt, bis sie nur noch zweireihig ist. Im letzteren Fall kann der 
Wert leicht ermittelt werden. 

Die Determinante wird im Programm in Form von DATA-Werten eingelesen. 

Als Beispiel wird die Determinante 

1 2 3 4 5 6 

2 2 3 4 5 6 

3 3 3 4 5 6 

4 4 4 4 5 6 

5 5 5 5 5 6 

6 6 6 6 6 6 

berechnet. Das Programm liefert den Wert —6. 
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100 REM Determinante 
110 : 

120 CLS:PRINT STRING$(25,42) 

130 PRINT" Determinante" 

140 PRINT STRINGS(25,42) 

150 : 

160 READ N:'Ordnung der Determinante 
170 DIM A(N,N) 

180 : 

190 'Einlesen der Derterminante 
200 FÜR 1=1 TO N 
210 FOR J=1 TO N 

220 READ A ( I , J) :PRINT A<I,J); 

230 NEXT J:PRINT 
240 NEXT I:PRINT 

250 : 

260 REM Berechnung der Determinante 
270 GOSUB 310 

2B0 PRINT"Determinante=";D 
290 END 
300 : 

310 REM Unterprogramm Determinante 

320 D=1:M=N 

330 WHILE M>2 

340 P=0:REM Pivotsuche 

350 IF ABS(A(1,1))>.00001 THEN 490 

360 FOR K=2 TO N 

370 IF ABS < A <1,K))< =ABS(A(1,K—1)) 

THEN 390 

380 F‘=A (1 , K > : L=K 

390 NEXT K 

400 IF P=0 THEN D=0:G0T0 640 
410 : 

420 REM Spaltentausch 

430 FOR J=1 TO M 

440 H=A(J,L):A(J,L)=A(J,1):A(J,1)=H 

450 NEXT J 

460 D=-D 

470 : 

480 REM Entwicklung nach der 1.Spalte 

490 FOR J=2 TO M 

500 FOR K=2 TO M 
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510 A(J,K)=A(1,1)*A(J,K)-A(J,1)*A(1,K) 

520 NEXT K 

530 NEXT J 

540 T=SBN (A (1,1)) *ABS (A (1,1) (M-2) ) : 

D=D/T 

550 FOR J=2 TO M 

560 FOR K=2 TO M 

570 A(J-1,K-1)=A(J,K) 

580 NEXT K 

590 NEXT J 

600 M=M-1 

610 WEND 

620 'Zweireihige Restdeterminante 
630 D=D*(A(1,1> *A(2,2)-A(1,2)*A(2,1)) 

640 RETURN 
650 : 

660 DATA 6 

670 DATA 1,2,3,4,5,6 
680 DATA 2,2,3,4,5,6 
690 DATA 3,3,3,4,5,6 
700 DATA 4,4,4,4,5,6 
710 DATA 5,5,5,5,5,6 
720 DATA 6,6,6,6,6,6 


************************* 
Determinante 

■*•*** ******** ** *********** 

1 2 3 4 5 6 

2 2 3 4 5 6 

3 3 3 4 5 6 

4 4 4 4 5 6 

5 5 5 5 5 6 

6 6 6 6 6 6 

Determinante=-6 
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30. MÜLLER-ITERATION 


Die Iteration ist eine Programmiermethode, bei der durch zahlreiche Iterations¬ 
schritte ein gesuchter Punkt immer näher eingekreist und schließlich mit der 
vorgegebenen Genauigkeit bestimmt wird. 

Iterationsverfahren werden hauptsächlich zur Lösung von nichtlinearen Glei¬ 
chungen genutzt, da z.B. Nullstellen von Polynomen vom Grad >4 nicht mehr 
formelmäßig berechnet werden können. Iterationen finden jedoch auch An¬ 
wendung bei der Lösung von großen linearen Gleichungssystemen (z.B. beim 
Gauß-Seidel-Verfahren). 


Ein einfaches Verfahren ist die sog. Fixpunkt-Iteration 
x n+1 = F (x n ) 

Dabei wird der Graph der Funktion F(x) mit der Geraden y = x zum Schnitt 
gebracht. 


Bei der Regula falsi 


*n+1 = 


xn - f(x n ) 


x n -x n -1 

f(x n ) - f(x n —i) 


wird die Nullstelle der Verbindungsgerade zweier Graphenpunkte von f(x) be¬ 
stimmt. 


Bessere Konvergenz erzielt das Newton-Verfahren 

f(x n ) 

x n +1 = x n - ——- 
f (x n ) 

bei dem jeweils die Nullstelle der Tangente bestimmt wird. Die Newton-Itera¬ 
tion läßt sich auf nichtlineare Systeme und somit auch auf komplexe Glei¬ 
chungen ausdehnen (siehe Programm 32). Programme zu den erwähnten Itera¬ 
tionsverfahren finden sich z.B. in [15] und [17], 

Die Müller-Iteration (1956) (in der englischsprachigen Literatur Müller ge¬ 
nannt) kann als Verallgemeinerung der Regula falsi aufgefaßt werden. Im Ge¬ 
gensatz zu dieser wird durch 3 Kurvenpunkte eine Parabel gelegt und deren 
Nullstelle bestimmt. Die Bestimmung der Parabelgleichung erfolgt mittels 
Interpolation. 
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Sind xi, X 2 und X 3 die 3 Kurvenpunkte bzw. die Startwerte, so berechnet 
man den nächsten Iterationswert aus 

x = X 3 + I (x 3 - X 2 ) 

Der Faktor I berechnet sich aus dem Quotienten 
— 2 Df(x 3 ) 

G ± N /G 2 - 4 DCf(x 3 ) 

Die Konstanten ergeben sich aus 

G = L2 f(x 1 ) - D2 f( x 2> + (L+D) f(x 3 ) 

C = L(Lf(xi) — D f(x 2 ) + f(x 3 )) 
mit 

L= x 3 ~ x 2 
x 2 - xi 


x 2 - X 1 

Das Vorzeichen im Nenner von I wird aus numerischen Gründen so gewählt, 
daß der Nenner möglichst groß wird. 

Wegen der im Nenner auftretenden Wurzel könnte das Verfahren auch ins 
Komplexe gehen. Da die komplexe Newton-Iteration effektiver ist, ist das 
hier angegebene Müller-Verfahren auf das Reelle beschränkt. 


Zum folgenden Programm 

Die zu lösende nichtlineare Gleichung ist im Programm in Zeile 170 zu defi¬ 
nieren. Der Startwert, die Abbruchgenauigkeit und die max. Zahl der Iteratio¬ 
nen werden über INPUT eingegeben. 

Liegt der Startwert in der Nähe der gesuchten Nullstelle, so ist die Konver¬ 
genz gut. Bei Funktionen, die keine reellen Nullstellen haben, kann die Itera¬ 
tion manchmal zwischen zwei Werten hin- und herpendeln und so Konvergenz 
Vortäuschen. Zur Kontrolle wird daher im Programm noch der Funktionswert 
des letzten Iterationsschritts mit ausgedruckt. 

Als Programmbeispiel wird das kubische Polynom 
x3 — x — 1 

behandelt. Beim Startwert 1 ergibt sich in 5 Iterationsschritten die Nullstelle 
x = 1.324718 

auf Maschinengenauigkeit. Der zugehörige Funktionswert ist 
f(x) = 5.960465-IO -7 
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100 REM Muel1 eriteration 
110 : 

120 CLS:PRINT STRING*(25,42) 

130 PRINT" Mueller-Iteration" 

140 PRINT STRINGS(25,42) 

150 : 

160 'Definieren der Funktion 
170 DEF FNF < X)= X A 3-X-1 
1 SO PRINT"f (x) =x•"'•3 -m-1 " 

190 : 

200 INPUT"Startwert";XO 
210 INPUT"Abbruchgenauigkeit";EPS 
220 INPUT " 11 er at. i onsz ah 1 " ; IT 
230 s 

240 'Initialisieren 
250 K=1:D=.5 

260 X3-X0:X1=X3~D:X2=X3+D 
270 PRINT XO 
280 : 

290 'Muel1 er-Interpolation 
300 IF X2--X1 TUEN X2=X2+1E-0B 
310 LI-(X3-X2)/(X2—X1) 

320 D1 = (X 3— XI)/< X 2— X1) 

330 IF K>1 TUEN 350 
340 E1=FNF< X1):E2=FNF < X2) 

350 E3=FNF(X3) 

360 A1 -L1 •"••2*E 1 -D1 A 2*E2+ (L1+D1) *E3 
370 C1=L1*(L1*E1-D1*E2+E3) 

380 B=A 1 "'■2-4*D 1 *C 1 *E3 
390 : 

400 '11 eration 
410 IF EK0 TUEN B=Ö 

420 IF A1O0 TUEN A1 =A 1+SGN (Al) *SG!R (B) 
430 IF ABS(Al)+ABS(B)=0 TUEN A1=4*D1*E3 
440 IF A1=0 TUEN Al-.0000001 
450 L—2*D1*E3/A1: X-X3+L* (X3-X2) 

460 PRINT X 
470 : 

480 'Konvergenztest 

490 IF ABS(X—X3)<EPS*ABS(X3) THEN 560 
500 IF K<IT THEN 510 ELSE 530 


141 



510 

K=K+1: 

: X1—X2: 

X2=X3:X3=X 


520 

El—E2: 

: E2—E3: 

GOTO 300 


530 

PRINT 1 

'Kei ne 

Konvergenz 

nach";IT;"Iterationen" 

540 

PRINT* 

'Letzter "; 


550 

: 




560 

PRINT' 

‘Fun kt i 

onswert®"; 

FNF < X 3) 

570 

END 





*•*•*•**■**•**■**•**■•*•* *•*•**■*•**•#•** 
Muel 1 er-1 terat i on 
***•*•***********•#•#*■*■*•*•#*■** 

■f (x ) =x ‘"'3-x~l 
Startwert? 1 

Abbruchgenauigkeit? 1e—06 
Iterationszah1? 50 
1 

1.31344632 
1.32456943 
1.32471808 
1.32471796 

Funktionswert® 5.43427E-07 
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31. MÜLLER-ITERATION (ZWEIDIMENSIONAL) 


Die Müller-Iteration (Programm 31) läßt sich auch auf Funktionen 

f(x, y) 

zweier Variablen ausdehnen. 

Dabei wird die x-y-Ebene in einfachen Müller-Iterationsschritten rechteckig 
in x- und y-Richtung durchlaufen 

x = X 3 + h (*3 - *2> 
y = Y3 + <2 <Y3 - V2> 

Das Verfahren bricht ab, wenn die Änderung der x- bzw. y-Werte die vorgege¬ 
bene Abbruchgenauigkeit unterschreitet. 

Zum folgenden Programm 

Startwerte und Abbruchgenauigkeit werden mittels INPUT eingegeben. Das 
Programm ist gegen das Verschwinden des Nenners gesichert, jedoch nicht das 
Überschreiten des Wertebereichs (Overflow). 

Als Beispiel wird die Funktion 
f(x, y) = eV - x2 -2 
behandelt. 

Für die Startwerte x=1, y=1 erhält man in 2 Schritten die Lösung 
x = 0.84751501, y = 1 

Analog liefern die Startwerte x=1, y=0 die Lösung 
x = —0.7637183, y = 0.9490546 

Da die Funktion unendlich viele Nullstellen hat, erhält man für verschiedene 
Startwerte wechselnde Lösungen. Sie erfüllen die Gleichung 

y = In (x2 + 2) 

die man erhält, wenn man die gegebene Funktionsgleichung nach y auflöst. 
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100 REM Muelleriteration (2-dimensianal) 

110 : 

120 CLS:PRINT STRING#(39,42) 

130 PRINT" Mueller-Iteration (2~-dim.) 

140 PRINT STRING#(40,42) 

150 : 

160 'Definieren der Funktion 

170 DEF FNF ( X , Y) =EXP (Y ) — X "'2—2 

180 PR I NT: PR I NT " f ( x , y) =exp (y) -x •'-2-2 " 

190 : 

200 PRINT:INPUT "Startwerte";X0,Y0 
210 INPUT" Abbruchgenaui gkei t " ; EF : 'S 
220 INPUT " 11 er at i onsz ah 1 " ; IT; PR I NT 
230 : 

240 'Initialisieren 
250 K= 1: B1 =. 5: B2=E< 1 
260 : 

270 X=X0:Y=Y0 

280 X3=X0:X1=X3-B1:X2=X3+B2 
290 PRINT XO,YO 
300 : 

310 'Muel1er-Interpolation 
320 IF X2—X1 TUEN X2=X2+1E~08 
330 L1=(X3-X2)/(X2-X1) 

340 Dl=(X3-X1)/(X2-X1) 

350 IF K>1 THEN 370 

360 E1-FNF(X1,Y);E2=FNF(X 2,Y) 

370 E3=FNF(X3,Y) 

380 GOSUB 600 

390 B1=L*(X3-X2):X=X3+B1 

410 : 

420 'Konvergenztest 

430 IF ABS (Bl) +ABS (B2) <EF’S*ABS <X3) THEN 570 
440 X0=X:Y3=Y0:Y1=Y3-B2:Y2=Y3+B2 
450 IF Y2=Y1 THEN Y2=Y2+.0000001 
460 Ll=(Y3-Y2)/(Y2-Y1) 

470 Dl=(Y3-Y1)/(Y2-Y1) 

480 E1=FNF(X,Y1):E2=FNF(X,Y2):E3=FNF(X,Y3) 

490 GOSUB 600 

500 B2=L*(Y3-Y2):Y=Y3+B2 
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510 IF ABS(Bl)+ABS(B2)<EPS*ABS(X3) THEN 570 
520 IF K<IT THEN K=K+1:Y0=Y:GOTO 270 
530 : 

540 PRINT"Keine Konvergenz nach";IT;" 

Iterationen" 

550 PRINT"Letzter "; 

560 ; 

570 F'RINT"Funkti onswert = " ; FNF (X3,Y3) 

580 END 
590 : 

600 'Unterprogramm 
610 A1=L1 •"•2*E 1 -D1 •"■•2*E2+ (Ll+Dl) *E3 
620 C1=L1* <L1*E1~D1»E2+E3) 

630 B=A 1 "•■2-4*D 1 *C 1 *E3 
640 IF B<0 THEN B=0 

650 IF AlOO THEN A1 =A 1+SGN (Al) *SQR (B) 

660 IF ABS(Al)+ABS(B)=0 THEN A1=4*D1*E3 
670 IF A1=0 THEN Al=.0000001 
680 L=~2*D1*E3/A1 
690 RETURN 
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************************ ****** *** ****** 
Muel 1 er- I ter at. i on (2-dim.) 
**************************************** 

•f < :•{ , y) -exp <y)-x A 2-2 

Startwerte? 1,1 
Abbruchgenauigkeit? le—06 
11erationsz ah1? 30 

1 1 

.047515091 1 

F unktionswert = 0 


•f ( k , y) =e;< p (y) —x ""2--2 

Startwerte? 1,0 

Ab br uc hgen auigkeit ? 1e-06 

Iterationszahl? 30 


1 

-1 

-0.764252978 
-0.763719359 
-0.763718694 
Funktionswert= 


0 

1.20627398 
1.00921891 
0.949423254 
0.949054804 
5.68107E-07 
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32. NEWTOIM-ITERATION IM KOMPLEXEN 


Die Newton-Iteration im Reellen 
f(*n) 


x n+1 - x n 


f'( x n) 


kann auch auf die komplexe Ebene ausgedehnt werden 
fUn) 


z n+1 - z n _ 


f'( z n) 


Um eine komplexe Arithmetik zu vermeiden, wird die Funktion f und ihre 
Ableitung f in Real-und Imaginärteil zerlegt 

f (z) = u + iv 

f'(z) = U x + iuy 


Damit erhält man 

f U + iv UU X + VUy VU X — UUy 

f' U x + iUy U x 2 + Uy2 U x 2 + Uy2 

durch Erweitern mit dem Konjugiert-Komplexen des Nenners. Dabei stellen 


u x , u y 

die partiellen Ableitungen von u nach x bzw. y dar. Trennt man Real- und 
Imaginärteil, so erhält man die komplexe Newton-Iteration 


x n+1 


UU X + VUy 

u x 2 + Uy2 


Yn+1 = Yn - 


VU X — UUy 

U X 2 + Uy2 


Polynome können mit Hilfe der binomischen Formel in Real- und Imaginär¬ 
teil zerlegt werden. Z.B. gilt 

z3 + 1 = (x + iy)3 + 1 

= x3 + 3ix2y — 3xy2 — iy3 + i 


Dies liefert den Realteil u = x3 — 3xy2 + 1 und Imaginärteil v = 3x2 y—y3. 
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Andere Funktionen können mit Hilfe eines geeigneten Tafelwerks zerlegt wer 
den. Nach [4] gilt z.B. 

sin z = sinx coshy + i cosx sinh y 
Dabei treten die hyperbolischen Funktionen sinh bzw. cosh auf. 


Zum folgenden Programm 

Der Realteil u bzw. Imaginärteil w der Funktion muß im Programm als Un 
terprogramm in den Zeilen 170—180 vereinbart werden. Die partiellen Ablei 
tungen von u nach x und y folgen in den Programmzeilen 230—240 

Die Startwerte, die Abbruchgenauigkeit und die maximale Zahl der Iteratio 
nen wird über INPUT-Anweisungen eingegeben. 

Als Programmbeispiel wird die Gleichung 
z3 + 1 

gelöst. Die Zerlegung in Real- und Imaginärteil ist 
u = x3 — 3xy2 + 1, v = 3x2y — y3 
die partiellen Ableitungen von u nach x bzw. y sind 
u x = 3x2 _ 3y2, u y = —6xy 

Eingabe der Startwerte x=1, y=1 liefert die Nullstelle 
0.5 + i * 0.8660254 

Entsprechend liefern die Startwerte x=1, y=—1 die Nullstelle 
0.5 - i * 0.8660254 
bzw. x=1, y=0 die Nullstelle 
—1 + i * 0. 

Die exakten Nullstellen sind hier 
—1 und 
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100 REM Newton-Verfahren im Komplexen 
110 : 

120 CLS:PRINT STRING*(30,42) 

130 PRINT" Komplexe Newton-Iteration" 

140 PRINT STRING*(30,42) 

150 : 

ISO 'Definieren Real- u.Imaginaertei1 
170 DEF FNU ( X , Y) =X""’3~3*X*Y'"'2+1 
180 DEF FNV(X,Y)=3*X Ä 2*Y-Y A 3 
190 PR INT"Rea 11ei 1 =x A 3-3*x *y-"2+1 " 

200 PRINT"Imaginaertei1=3*x rt 2*y-y A 3" 

210 : 

220 'Def. der part.Ableit. von u(x,y) 

230 DEF FNX ( X , Y) =3*X ""2-3*Y""'2 
240 DEF FNY(X,Y)=-6*X*Y 
250 : 

260 PRINT:INPUT"Startwert Real,Imag.";X0,Y0 

270 INPUT"Genauigkeit";EPS 

280 INPUT"max.IterationszahlIT 

290 : 

300 "An f an g swert e 
310 K—0 

320 PRINT X0,Y0 
330 : 

340 'Iteration 
350 K-K+1 

360 LJX—FNX ( XO , YO) : UY=FNY ( XO , YO) 

370 G=U X 2+UY2 

380 IF G=0 TUEN PRINT" 

Gradient ist Nu11vektorEND 
390 X--X0+ (FNV ( XO , YO) *UY-FNU ( XO, YO) *UX ) /G 
400 Y-Y0-- (FNV ( X0 , YO) *UX+FNU ( X0 , YO > *UY)/G 
410 PRINT X , Y 
420 : 

430 ' l<onverqenz test 

440 IF ( X — X5) 2 + (Y—Y0) ’2 < EF'S*ABS ( X ) THEN 490 
450 IF K<IT THEN XO-X:Y0=Y:GOTO 350 
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460 F'RINT"Kei ne Konvergenz nach ";IT;" 
I ter at i orien " 

470 PRINT"Letzter WertGOTO 500 
480 : 

490 PRINT"Nu11stel1e="; 

500 IF Y<0 THEN 520 
510 PRINT X;"+i*";Y:END 
520 PRINT X;"-i*";ABS<Y) 

530 END 


■*•***■*•** ***••*■*■*•*•*•*•*■*•**•** •**•**•*•** *• 
Komplexe Newton-Iteratian 
■****************************** 

Real tei 1 =>: ■'"•3~3*x*y'"2+1 
Imagi naertei 1 =3*x'" - 2*y-y""3 

Startwert Real,Imag.? 1,1 
Genauigkeit? le-06 
max.Iterationsz ah1? 20 
1 1 

0»666666667 O.833333333 
0.508691916 0.841099874 

0.499329996 0.866269172 

0.499999911 0.866024903 

0.5 0.866025404 

Nullstelle= .5 +i* 0.866025404 
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****************************** 
Komplexe Newton-Iteration 
****************************** 

Real tei 1 — ""'3—3*x*y""2+1 
I mag i n aer t eil =3*x ""•2*y -y'"'3 


Startwert Real,Imag.? 1,-1 
Genauigkeit? le-06 
max.Iterationszahl? 20 


1 

O . 666666667 
0.508691916 
0.499329996 
0.499999911 
0.5 

Null stelle— . 


-1 

-0.833333333 
-0.841099874 
-0.866269172 
-0.866024903 
-0.866025404 
5 -i* 0.866025404 


Start wert Real , Imag. 7' 1,0 
Gen auig keit? 1e—06 
max . 11er at. i orisz ah 1 ? 20 
1 0 

0 . 333333333 0 


”2.77777778 0 
-1.89505185 0 
-1.35618683 0 
-1.08535856 0 
-1.00653708 0 
-1.00004236 0 
-1 0 


Nul1stel1e=—1 +i* 0 
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33. NEVILLE-INTERPOLATION 


Unter Interpolation versteht man Aufsuchen von Zwischenwerten einer tabel¬ 
lierten Funktion. Im Gegensatz dazu wird bei der Extrapolation ein Wert 
außerhalb des tabellierten Bereichs gesucht. 

Sind die Funktionswerte Z(I) an den Stellen X(l) gegeben und soll ander 
Stelle X interpoliert werden, so läßt sich die Neville-Interpolation wie folgt 
beschreiben 

260 FOR K=1 TO N 
270 FOR l=K+1 TO N 
280 IF l=N+1 TH EN 330 

290 Z(I) = (Z(l) * Y(K) - Z(K) * Y(I)) / (Y(K) - Y(l) 

300 NEXT I 
310 NEXT K 

Dabei gilt 

Y(l) = X - X(l) 

N ist die Anzahl der vorgegebenen Funktionswerte, Z(N) der gesuchte, inter¬ 
polierte Funktionswert. 

Die in Programmzeile 340 auftretenden Differenzen heißen dividierte Diffe¬ 
renzen. Mit ihrer Hilfe kann auch das Polynom vom Grad höchstens N—1 be¬ 
stimmen, das durch die vorgegebenen Punkte geht. Ein Programm zur Newton- 
Interpolation findet sich in [17], 


Zum folgenden Programm 

Im Programm wird nach dem Neville-Verfahren im angegebenen Punkt inter¬ 
poliert. 

Folgende Funktionsstellen (xly) seien gegeben: 

(-4 1-26), (-1 14), (0 12), (2 116) 

Gesucht werde der Funktionswert im Punkt x=1. 

Das Programm liefert den Funktionswert 
f(1) = 4 

Das interpolierende Polynom ist hier 
p(x) = x3 + 2x2 — x + 2 
wie man leicht nachrechnet. 
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100 REM Nevi11e—Interpolation 
110 : 

120 CLS:PRINT STRING#(25,42) 

130 PRINT" Neville—Interpolation" 

140 PRINT STRINB#(25,42) 

150 ; 

160 INPUT"Wieviele Stuetzstel1en";N 
170 DIM X(N) ,Y(N) ,Z (N) 

180 : 

190 PRINT:PRINT"Eingabe der Stuetzstellen" 

200 FOR 1=1 TO N 

210 PRINT IJSING" ;•:<##), y (##)"; I , I; : 

INPUT X(I) ,Y (I) 

220 NEXT I 
230 : 

240 PR INT: INF'UT"x - Wert " ; X 

250 FOR 1=1 TO N 

260 Z(I)=Y(I):Y <I)=X — X <I) 

270 NEXT I 
280 : 

290 'Interpolation 
300 FOR K=1 TO N 
310 FOR I=K+1 TO N 
320 IF I=N+1 THEN 380 

330 IF Y(I)=Y(K) THEN PRINT"y-Werte 

muessen verschieden sein":END 
340 Z <I) = (Y <K>*Z(I)-Y(I)*Z (K) )/(Y <K>-Y(I)) 

350 NEXT I 
360 NEXT K 
370 : 

380 PRINT: PRINT" Interpol . Funk.tionswert = " ; Z (N) 
390 END 


154 



***■**•**■*-**-*•* *■***-**■■*■*•*•* ■*•*•* 
Nevi11e—Interpolatian 
#**•*■**•*•*-***•*■*•**•■****■-****** 
Wieviele Stuetzstel1en? 4 

Eingabe der Btuetzstellen 
x< 1 ), y ( 1)? -4,-26 
x( 2),y( 2)? -1,4 
x( 3),y( 3)? 0,2 
x < 4) ,y( 4)? 2,16 

x.Wert? 1 

Interpol.Funktian swert= 4 
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34. GAUSS-JORDAN-VERFAHREN 


Löst man im Gleichungssystem 

an xi +ai 2 x 2 + ....+ai n x n = bi 
a 21 xi + a 22 X 2 + . . . - + a 2 n x n = b 2 


a nix +a n 2 X2 + • • • • + a nn x n = b n 
die k-te Gleichung nach der Unbekannten xk auf, so erhält man 

xk = (—akl xi - ak2 X2 ---akn x n + bk) / akk 

Dies setzt voraus, daß das Diagonalelement akk ^ 0 ist- ist a kk = 0 ober sehr 
klein, so führt man eine entsprechende Spalten- oder Zeilenvertauschung 
durch. Findet sich keine Zeile oder Spalte, die ein Element =£0 liefert, so ist 
die entsprechende Matrix singulär, d.h. die Determinante der Matrix ver¬ 
schwindet und das Gleichungssystem ist nicht eindeutig lösbar. 

Setzt man diese Unbekannte xk in alle anderen Gleichungen ein, so verschwin¬ 
det sie aus diesen Gleichungen, d.h. sie wird eliminiert. Führt man dies für 
alle Unbekannten durch, so bleibt in der 1. Zeile nur noch die xi, in der 2. 
noch X2 und entsprechend in der k-ten Zeile xk übrig. 

Man hat somit ein Gleichungssystem erhalten, in dem die Unbekannten nur 
noch in der Diagonale der Matrix stehen 

dl 1 xi =ci 

d22 X2 = C2 

d33 X3 = C3 


dnn x n - Cn 


Dieses Gleichungssystem kann sehr einfach gelöst werden 


x k = 


Ck 

dkk ' 


d k k^0 
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Lineare Gleichungssysteme finden zahlreiche Anwendungen bei 

— Dimensionierung elektrischer Schaltkreise 

— der Lösung von partiellen Differentialgleichungen 

— Vektorrechnung 

— Input- und Outputmodellen der Wirtschaft 

— Ausgleichsrechnung 

— elastische Verformungen von Fachwerken, Karrosserien 
usw. 

Zur Konstruktion eines Flugzeugflügels treten z.B, Gleichungssysteme mit 
mehr als 20.000 Unbekannten auf. 


Zum folgenden Programm 

Das Programm löst beliebige lineare Gleichungssysteme, bei denen die Zahl 
der Unbekannten mit der Zahl der Gleichungen übereinstimmt. 

Gibt es mehr Unbekannte als Gleichungen, so heißt das System unterbe¬ 
stimmt, im umgekehrten Fall überbestimmt. Zur Lösung von überbestimmten 
Gleichungssystemen, wie sie bei der linearen Regression auftreten, gibt es spe¬ 
zielle Verfahren wie z.B. die Methode der kleinsten Quadrate. 


Als Programmbeispiel wird das Gleichungssystem mit 6 Unbekannten behan¬ 
delt. 

6 xi — 3x2 + 2 x 3 + M — X 5 + xß = 11 

—3xi — 7x2 + 4 x 4 — 2xß + X 7 = —5 

4xi — 3x2 + 6 X 3 — X 4 + 2 x 5 + X 7 = 28 

2 xi + 4x2 + 5 x 3 — 7 x 4 — 3 x 5 + 2 x 7 = —6 

—xi + 5 x 2 — 4 x 3 + 3 x 5 — 2 x 7 = 25 

3xi + 4 x 3 — 2 x 4 + 5 x 5 ~ 6 x 7 = —4 

Die Anzahl der Unbekannten und die Gleichungsmatrix mit rechter Seite wird 
im Programm in Form von DATA-Werten eingelesen. 

Das Programm liefert hier die exakten Lösungen 

x 1 = 1, X 2 = 2, X 3 = 3, X 4 = 4, X 5 = 5 und xß = 6 . 
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100 REM GauBB—Jordan--Ver-fahren 
110 : 

120 CLS:PRINT STRING*<27,42) 

130 PRINT" Gauss—Jordan—Ver-fahren" 

140 PRINT STRING*(27,42) 

150 : 

ISO DEFINT I--N 

170 READ N:'Anzahl der Unbekannten 
180 DIN A(N,N+l),X<N) 

190 : 

200 PRINT"Gegebenes Gl eichungsystem:" 

210 FÜR 1=1 TO N 

220 FÜR J=1 TO N+l 

230 READ A(I,J):PR INT A <I,J); 

240 NEXT JSPRINT 
250 NEXT I SPRINT 
260 s 

270 FÜR K=1 TO N 

280 IF K=N TUEN 370 

290 G=A(K,K)sM=K 

300 FOR I=K+1 TO N 

310 IF ABS(A(I,K))< =G THEN 330 

320 G=ABS(A(K,K))s M=I 

330 NEXT I 

340 FOR J =K TO N+l 

350 H=A(K,J)sA(K,J)=A(M,J):A(N,J)=H 

360 NEXT J 

370 IF A BS(A(K,K) )<.0000001 THE N 

PR I NT"Matri x si ngulaer"s END 
380 FOR 1=1 TO N 
390 IF I=K TUEN 440 

400 F=A<I,K>/A(K,K) 

410 FOR J=K+1 TO N+l 

420 A(I,J)=A(I,J)-A(K,J)*F 

430 NEXT J 

440 NEXT I 

450 NEXT K 

460 s 

470 PRINT"L.oesung: " 

480 FOR 1=1 TO N" 

490 X(I)=A(I,N+l)/A(I,I) 

500 PRINT"x(";I;")=";X(I) 
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510 NEXT I 
520 END 
530 : 

540 DATA 6 

550 DATA 6,-3,2,1,-1,1,11 
560 DATA -3,-7,0,4,-2,1,-5 
570 DATA 4,-3,6,-1,2,1,28 
580 DATA 2,4,5,-7,-3,2,-6 
590 DATA -1,5,-4,0,8,-2,25 
600 DATA 3,0,4,-2,5,-6,-4 


*************************** 
Ganss-Jordan-Verfahren 
*************************** 
Gegebenes Gleichungsystem: 

6-3 2 1 -1 1 11 

-3 -7 0 4-2 1-5 

4 -3 6-1 2 1 28 

245 —7 —3 2 —6 

-1 5 -4 0 8 -2 25 

3 0 4 -2 5 -6 -4 

Loesung: 

x< 1 )= 0.999999998 
;•! ( 2 )- 2.00000002 
:•; ( 3 )= 3.00000002 
( 4 )— 4.00000002 
x < 5 ) = 5 
x ( 6 ) — 6 


160 



35. MAT RIZEN-IN VERSION NACH FADDEJEW 


Eine einfache Methode zur Invertierung einer quadratischen Matrix, die zu¬ 
gleich noch das charakteristische Polynom liefert, ist das Verfahren von 
Faddejew [9]. 

Ist die zu invertierende Matrix Al von der Ordnung n, so wird iterativ die 
Matrix A n in n—1 Schritten über die Gleichung 

Aj+i = Ai (Aj - q E) 

berechnet; dabei ist E die Einheitsmatrix und cj = Spur (Aj)/i. Die Spur einer 
Matrix ist die Summe der Diagonalelemente. 

Die gesuchte inverse Matrix ergibt sich dann aus 
Af 1 = (A n —1 c n —1 E) ~ 


Der Wert c n liefert außerdem noch den Wert der Determinante 
Det(A) = (—1 )n+1 Cf1 

Programme zu anderen Inversionsverfahren finden sich z.B. in [17], 


Zum folgenden Programm 

Das Programm berechnet die inverse Matrix nach dem angegebenen Verfahren. 
Die zu invertierende Matrix wird in Form von DATA-Werten eingelesen. 


Als Beispiel soll die Inverse von 

( 5 7 6 5 \ 

7 10 8 7 \ 

6 8 10 9 I 

5 7 9 10 / 

berechnet werden. Das Programm liefert hier die exakte Inverse 



Durch Matrizen-Multiplikation (Programm 22) kann geprüft werden, daß gilt 
AA“ 1 = A —1 A = E. 
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100 REM Matrizeninversion nach Fadelejew 
110 : 

120 CLSSPRINT STRING*<25,42) 

130 PRINT" Matrizeninversion" 

140 PRINT STRING*(25,42) 

150 : 

160 DEFINT I-N 

170 READ N:'Ordnung der Matrix 

180 IHM A(N,N) ,A1 <N,N) ,B(N,N> ,G(N,N) ,H(N,N) 

190 : 

200 PRINT"Gegebene Matrix:" 

210 FOR 1=1 TO N 
220 FOR J=1 TO N 
230 READ A(I,J)SPRINT A(I,J>; 

240 Al(I,J)=A<I,J) 

250 NEXT JSPRINT 
260 NEXT I SPRINT 
270 : 

280 FOR 1=1 TO N 
290 S=0 

300 FOR J=i TO N 

310 S=S+A(J,J> 

320 NEXT J 

330 OS/I 

340 : 

350 FOR J=1 TO N 

360 FOR K—1 TO N 

370 IF J =K THEN H < J,J)=A(J,J)-C 

ELSE H(J,K)=A(J,K) 

380 NEXT K 

390 NEXT J 

400 : 

410 IF I=N—1 THEN 420 ELSE 480 
420 FOR J=i TO N 

430 FOR K=1 TO N 

440 G(J,K)=H(J,K) 

450 NEXT K 

460 NEXT J 

470 : 

480 FOR J=1 TO N 

490 FOR K=1 TO N 

500 S=0 


162 



510 FOR L=1 TO N 

520 S=S+A1<J,L)*H<L,K) 

530 NEXT L 

540 A(J,K)=S 

550 NEXT K 

560 NEXT J 

570 NEXT I 
580 : 

590 IF C=0 THEN PRINT "Matrix si rigul aer ": END 
600 PRINT"Inverse Matrix:" 

610 FOR 1=1 TO N 
620 FOR J=1 TO N 
630 G<I,J)=G(I,J)/C 
640 PRINT G <I,J); 

650 NEXT JSPRINT 

660 NEXT I SPRINT 

670 D= (-1 ) (N+l) *C 

680 PRINT"Determinante 3 ";D 

690 END 

700 : 

710 DATA 4 
720 DATA 5,7,6,5 
730 DATA 7,10,8,7 
740 DATA 6,8,10,9 
750 DATA 5,7,9,10 
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*********•************•*-**•* 
Matrizeninversion 
************************* 
Gegebene Matrix: 

5 7 6 5 

7 10 8 7 

6 8 10 9 

5 7 9 10 


Inverse 

Matrix: 

68 -41 

-17 10 

-41 25 

10 -6 

-17 10 

5 ~“3 

10 ~6 ■ 

t n 

“ ■Jt X- 

Determinante= 1 
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36. CHARAKTERISTISCHES POLYNOM NACH FADDEJEW 


Das charakteristische Polynom hat genau die Eigenwerte der Matrix als Null¬ 
stellen. 

Eigenwerte X und Eigenvektoren x sind Lösung der Gleichung 
Ax = X x 

Eigenwerte von Matrizen geben an 

— Klassifikation von Flächen n-ter Ordnung 

— Schwingungsfrequenzen eines Systems 

— Kreiseleigenschaften von festen Körpern 

— Stationäres Verhalten von Markow-Ketten 
usw. 

Wie beim vorhergehenden Programm entstehen beim Faddejew-Verfahren [9] 
Ai+i = Ai (Aj — cj E) i = 1, 2, . .., n-1 

die Koeffizienten q aus 
Ci = Spur (Aj) /i 

Versieht man die Werte cj mit alternierenden Vorzeichen, so erhält man die 
Koeffizienten des charakteristischen Polynoms der Matrix. Bestimmt man die 
Nullstellen des Polynoms mit Hilfe der Programme 17, 18, 31,32 oder 33, so 
erhält man die gesuchten Eigenwerte der Matrix. 

Programme zu anderen Eigenwertverfahren, die auch gleichzeitig die Eigen¬ 
vektoren bestimmen, finden sich in [17], 


Zum folgenden Programm 

Analog zum vorhergehenden Programm werden die Koeffizienten cj bestimmt 
und mit den richtigen Vorzeichen versehen. Die Matrix wird wieder in Form 
von DATA-Werten eingelesen. 
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Als Beispiel wird die Matrix 

4 4 1 \ 

6 1 4 \ 

16 4 

4 4 6 / 

behandelt. Das Programm liefert das charakteristische Polynom 
p(x) = x^ — 24x3 + 1 50x2 _ 200x — 375 

dessen Nullstellen bereits in Programm 17 bestimmt worden sind. Die Eigen¬ 
werte von A sind somit 

Xi = —1 ; X 2 = 5 ; X 3 = —5 ; A 4 = 15. 
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100 REM Charakt - Pol ynom nach Faddejew 
110 : 

120 CLS:PRINT STRING*(30,42) 

130 PRINT" Charakt. Pol ynom einer Matrix" 
140 PRINT STRING*(30,42) 

150 : 

160 READ N:'Ordnung der Matrix 
170 DIM A(N,N),Al(N,N),G(N,N),H(N,N),C(N) 
180 : 

190 FORINT "Gegebene Matrix:" 

200 FOR 1=1 TO N 
210 FOR J = 1 TO N 

220 READ A<I,J>:PRINT A(I,J>; 

230 Al(I,J)=A(I,J) 

240 NEXT J:PRINT 
250 NEXT I SPRINT 
260 : 

270 FOR 1=1 TO N 
280 S=0 

290 FOR J=1 TO N 

300 S--S+A (J , J ) 

310 NEXT J 

320 C(I)=S/I 
330 : 

340 FÜR J=1 TO N 

350 FOR K=1 TO N 

360 IF J=k! THEN H (J , J ) =A (J , J ) -C (I) 

ELSE H(J,K)=A(J,K) 

370 NEXT K 

380 NEXT J 

390 : 

400 IF I=N—1 THEN 410 ELSE 470 
410 FOR J=1 TO N 

420 FOR K=1 TO N 

430 G (J , Fs:) =H (J , Fs:) 

440 NEXT K 

450 NEXT J 

460 : 

470 FOR J=1 TO N 

480 FOR K=1 TO N 

490 S=0 

500 FOR L=1 TO N 
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510 S=5+Al(J,L>*H(L,K) 

520 NEXT L 

530 A < J,K)—S 

540 NEXT K 

550 NEXT J 

560 NEXT I 
570 : 

580 PRINT"Koe-f -F . des charakt. Pol ynoms: 
590 IF (N MOD 2)=0 TUEN S=-l ELSE S=1 
600 PRINT-S; 

610 FÜR 1=1 TO N 
620 Cd) =S*C ( I) 

630 PRINT C(I); 

640 NEXT I:PRINT 
650 END 
660 : 

670 DATA 4 

680 DATA 6,4,4,,! 

690 DATA 4,6,1,4 
700 DATA 4,1,6,4 
710 DATA 1,4,4,6 


•****•*•*■** **-***4e* ■*****•*****■•*■■*•*•** 

Charakt. Pol ynom einer Matri;-: 
*-*-*■***■*•#•*•*•**■*■*■*•* •***-*•*•*•* *■*■**••# 
Gegebene Matri : 

6' 4 4 1 

4 6 1 4 

4 1 6 4 

1 4 4 6 

Koe-f -f » des charakt. Polynoms: 

1 -24 150 -200 -375 
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37. NUMERISCHE DIFFERENTIATION 


Obwohl Funktionen — im Gegensatz zur Integration — analytisch differenziert 
werden können, ist dennoch das folgende Programm von Interesse, da es einen 
Überblick über Funktions- und Ableitungswerte im ganzen Intervall ermög¬ 
licht. 


Das Programm benützt die 3-Punkte-Formeln zur Ableitung: 


f'(x) 


—3f (x) + 4f(x+h) — f (x+2h) 
2 h 


f"(x> 


(f(x) — 2f (x+h) + f(x+2h) 
h2 


Da die Formeln nach links "übergreifen", ist zu beachten, daß die Funktion 
auch links vom betrachteten Intervall definiert sein muß. 

Da in BASIC die Verkettung von Funktionen implementiert ist, können die 
beiden Formeln direkt als Funktionen mittels 

DEF FN. . 
vereinbart werden. 

Dazu muß nur in Zeile 170 die zu differenzierende Funktion als 
DEF FNF(x) = . . . 
definiert werden. 

Das Intervall und die Schrittweite h werden über eine INPUT-Anweisung ein¬ 
gegeben. Wegen der unvermeidlichen Rundungsfehler werden alle Funktions¬ 
werte im Programm auf 4 Dezimale gerundet. Dies kann, wenn gewünscht, 
leicht geändert werden. Die Schrittweite zur Ausgabe ist im Programm mit 
der Variablen Fl (Zeile 370) gegeben; diese kann notfalls geändert werden. 


Zum folgenden Programm 

Als Programmbeispiel wird die Funktion 
f(x) = x^e -x + sinx 

gewählt. Das Intervall soll [1,2], die Schrittweite 0.1 sein. Wie man dem Pro¬ 
grammausdruck entnimmt, wächst der Funktionswert bis x = 1.6 an und 
nimmt rechts davon wieder ab. Die Funktion hat somit zwischen 1.6 und 1.7 
ein Maximum, wie es auch am Vorzeichenwechsel der 1. Ableitung erkennbar 
ist. Die 2. Ableitung ist im ganzen Intervall negativ, d.h. der Graph ist rechts¬ 
gekrümmt. 
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100 REM Numerische Di -ff erenti ati on 
1 10 : 

120 CLS:PRINT STRING#(25,42) 


IN 
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Numerische Differentiation 
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38. NUMERISCHE INTEGRATION (TRAPEZREGEL) 


Da viele Funktionen wie 


sin (x2), e~ x ^, 


sinx 

x 


nicht analytisch integriert werden können, ist man hier auf numerische Inte¬ 
gration angewiesen. Es gibt zahlreiche Integrationsmethoden. Ein bekanntes 
Verfahren ist die Trapezregel (auch Sehnentrapez-Regel) genannt: 


b 

J f (x)dx 


= [f(a) + 2f(a+h) + 2f(a+2h) + . . + 2f(b—h) + f(b)] 


_g 

mit der Schrittweite h = ——, n Zahl der Intervallunterteilungen. 


Die Trapezregel eignet sich auch zur Extrapolation (siehe [33], [35], [36]). 
Weitere Programme zur numerischen Integration finden sich in [15] und [17], 


Zum folgenden Programm 

Das Programm berechnet das bestimmte Integral einer Funktion nach der an¬ 
gegebenen Trapezformel. 

Die zu integrierende Funktion ist in Zeile 170 mittels 
DEF FNF(X) = . .. 

zu definieren. Die Intervallgrenzen a, b und Anzahl n der Unterteilungen wer¬ 
den über IN PUT eingegeben. 

Als Beispiel wird das Integral 



0 

mit 50 Unterteilungen berechnet. Das Programm liefert 
3.141526 

Vergleicht man dies mit dem exakten Wert 
Tr = 3.14159265 


so liefert dies den Fehler 
7 • 10-5. 
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100 REM Numerische Integration 
110 : 

120 CLS:PRINT STRING*(25,42) 

130 PRINT" Trapezregel" 

140 PRINT STRINGS(25,42) 

150 : 

160 'Definieren der Integrandenfunktion 
170 DEF FNF (X) =4/(X--2+1 ) 

180 PR INT " f (;s) =4/ ( A 2+1 > " 

190 : 

200 INPUT "Interval1 grenzenA,B 
210 IF A>B THEN H=A:A=B:B=H 

220 INPUT "Anzahl der Interval11ei 1ungenN 

230 : 

240 H=(B~-A)/N: 'Schrittweite 
250 I=FNF(A)+FNF(B) 

260 S=0 

280 FÜR K=1 TO N-l 
290 S=S+FNF(A+K*H) 

300 NEXT K 

310 INTEGRAL=(I+2*S)*H/2 

320 : 

330 PRINT"Integral=";INTEGRAL 
340 END 


•**■***********•*•*■*•*••***■***■#■* 

Trapezregel 

****•******■**•*•*•*•■*■*■*■*•***•*■** 

•f <;•;) =4/ (:■; "-2+1) 

Interval1 grenzen? 0,1 

Anzahl der Interval 1tei1ungen? 50 

Integral= 3.14152599 
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39. HAMMING-VERFAHREN 


Bekannte Verfahren zur Integration von gewöhnlichen Differentialgleichun¬ 
gen sind die Verfahren von Runge-Kutta und Adams-Bashforth (siehe [33], 
[35], [37]. Programme dazu finden sich z.B. in [17]. 

Ein neueres Verfahren stammt von R.W. Hamming (1959) [12]. 

Die Differentialgleichung 

y' = f(x,y) 

mit dem Anfangswert 

y(a) = y0 mit xe [a,b] 

habe die Lösungsfunktion u(x). Im folgenden wird 
U(Xj) = Uj, f(xj, Uj) = fj 

gesetzt. Das Hamming-Verfahren läuft wie folgt in 5 Schritten ab: 

(1) Zunächst wird mit dem Prädiktor der Funktionswert 

vj+i = u;—3 + | h (2f j - fj_i + 2fj_2) 
vorhergesagt. 

(2) Dieser Prädiktorwert wird mittels des vorhergehenden Korrektorwerts zj 
modifiziert 

112 , . 
wi+1 = vj+i - (zi — v;) 

(3) Für diesen modifizierten Wert wird die rechte Seite der Differentialglei¬ 
chung bestimmt 

w'j+i = f(xj+i, Wj+l) 

(4) Der neue Korrektor lautet 

zj+i = 1 (9uj - uj—2 + 3h (wj+i + 2fj - fj_i)) 
o 

(5) Der gesuchte Funktionswert uj+i wird schließlich mit Hilfe des Korrek¬ 
tors verbessert 

g 

Ui + 1 =Zj+1 + (Zj + 1 -Vj+l) 
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Ist der rechte Rand des Integrationsintervalls [a, b] noch nicht erreicht, so 
wird x um die Schrittweite h erhöht und das Hamming-Verfahren erneut 
angewandt. 

Wie man der Prädiktorformel entnimmt, benötigt man bei jedem Schritt die 

Kenntnis der Funktionswerte fj, fj i und fj 2- Man sa 9 l - das Verfahren ist 

dreifach zurückgreifend. Da aber nur der Startwert ug und somit auch fg be¬ 
kannt ist, müssen die fehlenden Funktionswerte f-|, f 2 mit Hilfe eines anderen 
Verfahrens berechnet werden. Man wählt hier meist die Methode von Runge- 
Kutta. 


Zum folgenden Programm 

Als Beispiel soll die gewöhnliche Differentialgleichung 
y' = y - 2x/y 

mit dem Anfangswert y(0) = 1 im Intervall [0,1 ] integriert werden. 

Die Differentialgleichung wird in Zeile 170 in der Form 
170 F = Y - 2 * X/Y 

vereinbart. Das Intervall [a,b], der Anfangswert yg und die Schrittweite h 

werden in From von DATA-Werten eingelesen. Die numerische Lösung kann 

dem Programm-Ausdruck entnommen werden. 

Die exakte Lösung der Differentialgleichung lautet 
y = \/ 2x + 1 

Der relative Fehler am rechten Intervallrand ist somit 
3.6- 10- 6 . 
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100 REM Hamming-Verfahren 
110 : 

120 CLS:PRINT STRINB#(25,42) 

130 PRINT" Hamming-verfahren" 

140 PRINT STRING*(25,42) 

150 : 

160 'Differentialgleichung in der Form f '=f(;;,y) 
170 DEF FNF(X,Y)=Y-2*X/Y 
ISO PRINT TAB (5) " y ' -y-2*:-; / y " 

190 : 

200 INPUT " Interval 1 grenzen a,b";A,E< 

210 IF A>B THEN H=A:A=B:B=H 
220 INPUT "Anfangswert der Funktion";Y 
230 INPUT " Integrat. i ons-Schr i ttwei te" ; H 
240 DIM Y(5),F(5) 

250 : 

260 'Anfangswerte 
270 X=A:Y(1)=Y 

280 PRINT" ;< y" 

290 PRINT STRING*(25,45) 

300 PRINT USING " ##.# #####,######";X,Y 

310 : 

320 'Startprozedur nach Runge-Kutta 
330 F0R 1=1 TG 3 
340 K1=H*FNF(X,Y) 

350 K2=H*FNF(X+H/2,Y+Kl/2) 

360 K3=H*FNF(X+H/2,Y+K2/2) 

370 K4=H*FNF(X+H,Y+K3) 

380 Y=Y+(K1+2*K2+2*K3+K4)/6 

390 X=X+H:Y(1 + 1)=Y 

400 F(I+1)=FNF(X,Y) 

410 PRINT USING " ##.# #####.######";X,Y 

420 NEXT I 
430 : 

440 'Hamming-Verfahren 
450 K=0:P1=0:X = X +H 
460 UJHILE XOB+H/2 
470 'Praediktor 

480 P=Y(1)+4*H*(2*F(2)-F(3)+2*F(4))/3 

490 'Modifikator 

500 M=P—112*(K-P1)/121 
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510 Y=M:M=FNF(X,Y) 

520 'Kor r e k t or 

530 K=<9*Y(4)-Y(2)+3*H*<M+2*F(4)-F(3)))/8 
540 'Schrittfunktion 

550 Y(5)=K+9*(K—P)/121 

560 PRINT USINB " ##.# #####.######";X,Y 

570 F(5)=FNF(X,Y(5)):P1=P 

580 FOR 1=1 TO 4 

590 Y(I)=Y(I+1):F(I)=F(I+1) 

600 NEXT I 
610 X = X+H 

620 WEND 
630 END 


*■**•*•*■*•*•*•**■*•*■*•* *■*■*•* •*«•*• •#•*•■** 
Hamming-verfahren 
**********■#••*•*■*■*•*•**■*■*■*•**** 

y'=y-2*x /y 

Interval 1 grenzen a,b? 0,1 
Anfangswert der Funktion? 1 
Integrations-Schrittweite? .1 

y 


0 . 0 

1.000000 

0 . 1 

1.095446 

0.2 

1.183217 

0. 3 

1.264912 

0 . 4 

1.341564 

0.5 

1.414103 

0.6 

1.483179 

0.7 

1.549159 

0.8 

1.612438 

0.9 

1.673320 

1 . 0 

1.732057 
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ANHANG A 


HISTORISCHE RECHENAUFGABEN 

Altbabylonische Keilschrifttafeln (ca. 2000 v. Chr.) 

1) Länge und Breite habe ich multipliziert und so die Fläche gemacht. Wie¬ 
derum was die Länge über die Breite hinausgeht zur Fläche, habe ich addiert 
und (es ergibt) 183. Wiederum Länge und Breite addiert (gibt) 27. Länge, 
Breite und Fläche (ist) was?. Lösung: (15, 12, 180) (14, 13, 182) 

2) 1/4 Breite und Länge sind 7 Handbreiten, Länge und Breite zusammen 
sind 10 Handbreiten. Lösung (6,4) 

3) Die Fläche und die Seite des Quadrates habe ich addiert und 0,75 ist es. 
Lösung: 1/2 

4) Ein Balken, 0,5 lang (steht senkrecht) an einer Wand. Oben ist er um 0,1 
herabgekommen. Von unten (wie weit hat er sich entfernt)? Lösung: 3/10 

5) Finde, wie lange es dauert, bis eine Geldsumme sich bei 1/5 Jahreszins 
verdoppelt. Lösung: 3.801784 


Papyrus Rhind (ca. 1700 v. Chr.) 

6 ) Haufen: sein Ganzes, 2/3 hinzu, (von der Summe) 1/3 hinweg, bleibt 30. 
Lösung: 27 

7) 7 Häuser, in jedem 7 Katzen, jede frißt 7 Mäuse, jede von diesen 7 Ähren, 
von jeder Ähre könnte man 7 Scheffel Korn ernten, wieviele Dinge sind es ins¬ 
gesamt)? Lösung: 19607 

8 ) Zu verteilen 100 Brote unter 5 Personen in arithmetischer Folge, daß 1/7 
der Summe der drei ersten Anteile gleich der Summe der letzten drei Anteile 
ist. Lösung: (38 1/3, 29 1/6, 20, 10 5/6, 1 2/3) 

9) Ein Quadrat und ein weiteres, dessen Seitenlänge 3/4 des ersten ist, bil¬ 
den die Fläche 100. Wie lang ist die Seite des ersten? Lösung: 8 

10) Haufen, 2/3 davon, 1/2 davon, 1/7 davon, zusammen genommen ergeben 
33. Lösung: 14 28/97 


Anthologia Graeca (ca. 500 n. Chr.) 

11) Wieviele Äpfel werden gebraucht, wenn 4 von 6 Personen 1/3, 1/8, 1/4 
und 1/5 jeweils vom Ganzen erhalten, während die 5. Person 10 Äpfel erhält 
und genau ein Apfel für die 6. Person übrig bleibt. Lösung: 120 
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12) Demarchos verlebte 1/4 seines Lebens als Junge, 1/5 als Heranwachsen¬ 
der, 1/3 als Mann und verbrachte 13 Jahre im Alter. Wie alt ist er geworden? 
Lösung: 60 

13) Ziegelmacher, ich muß dringend an meinem Haus weiterbauen. Es ist 
schön heute und ich brauche nicht mehr als die üblichen 300 Ziegel. Du allein 
kannst soviele am Tage machen, dein Sohn 200 und dein Schwiegerson 250. 
Wenn ihr alle drei miteinander arbeitet, in welcher Zeit könnt ihr die ge- 
ten Ziegel hersteilen? Lösung: 2/5 

14) Ich bin ein Löwe aus Erz; aus meinen Augen, aus meinem Mund und aus 
meiner rechten Fußsohle springen Fontänen hervor. Mein rechtes Auge füllt 
das Becken in 2 Tagen, das linke in 3 und mein Fuß in 4 Tagen. Wie lange 
dauerts, wenn alles vereint? Lösung: 333/37 


Algebra des al-Hwarizmi (ca. 825 n. Chr.) 

15) Ein Quadrat und zehn seiner Wurzeln sind gleich 39. Lösung 3, —13 

16) Ein Quadrat und 21 sind gleich zehn seiner Wurzeln. Lösung 3,7 

17) Ich habe 10 in zwei Teile geteilt. Dann habe ich den einen Teil durch den 
anderen geteilt und der Quotient war 4. Lösung (8,2) 

18) Ich habe 10 in zwei Teile geteilt. Dann habe ich jeden Teil mit sich selbst 
multipliziert, das kleinere vom größeren abgezogen und der Rest war 40. 
Lösung (7,3) 


"Aufgaben zur Schärfung des Verstandes" von Alkuin (ca. 790) 

19) Wenn 100 Büschel Korn so auf 100 Leute verteilt werden, daß jeder 
Mann 3 Büschel, jede Frau 2 Büschel und jedes Kind 1/2 Büschel erhält, ist 
die Frage, wieviel Männer, Frauen und Kinder waren es? 

Lösung: (17, 5, 78) (14, 10, 76) (11, 15, 14) (8, 20, 72) (5, 25, 70) 

20) 30 Körbe, 10 davon voll, 10 halbvoll und 10 leer, sollen so unter drei 
Söhne verteilt werden, daß die Körbe und Inhalte gleich verteilt sind. Wie 
kann dies geschehen? Lösung: z.B. 10halbvolle 

21) Ein Hund jagt ein Kaninchen, das 150 Fuß Vorsprung hat. Der Hund 
springt immer 9 Fuß weit, wenn das Kaninchen 7 Fuß springt. Nach wieviel 
Sprüngen hat der Hund das Kaninchen eingeholt? Lösung: 75 

22) Ein Wolf, eine Ziege und ein Kohlkopf müssen in einer Fähre einen Fluß 
überqueren, die außer dem Fährmann noch einen weiteren Gegenstand beför¬ 
dern kann. Wie muß der Fährmann sie über den Fluß bringen, damit weder 
der Wolf die Ziege noch die Ziege den Kohl frißt? 
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Lilavati des Bhaskara II (1150) 

23) Das Quadrat aus 1/8 einer Herde Affen tummelt sich im Walde, während 
die übrigen zwölf auf dem Gipfel eines Hügels brüllen. Lösung: 16, 48 

24) Jemand hat 300 Rupien und 6 Pferde. Ein anderer hat 10 Pferde vom 
gleichen Wert, aber eine Schuld von 100 Rupien. Beide haben dasselbe Ver¬ 
mögen. Was ist der Preis eines Pferdes? Lösung: 57 1/7 

25) Schönes Mädchen mit den glitzernden Augen, sage mir, so du die richtige 
Kunst der Umkehrung verstehst: Welches ist die Zahl, die mit 3 vervielfacht, 
sodann um 3/4 des Produkts vermehrt, durch 7 geteilt, um ein Drittel des 
Quotienten vermindert, durch Ausziehung der Quadratwurzel, Addition von 
8 und Division durch 10 die Zahl 2 hervorbringt? Lösung: 28 

26) Von einem Schwarm Bienen läßt sich 1/5 auf einer Kadambablüte, 1/3 
auf der Silindhablume nieder. Der dreifache Unterschied der beiden Zahlen 
flog nach den Blüten einer Kutuja, eine Biene blieb übrig, die in der Luft hin 
und her schwebt, gleichzeitig angezogen von dem lieblichen Duft einer Jasmine 
und eines Pandamus. Sage mir, reizendes Weib, die Anzahl der Bienen! 

Lösung: 15 


Liber abaci des Leonardo von Pisa (1202),- 

27) Jemand reist in Geschäften nach Lucca, von da nach Florenz und von da 
heim nach Pisa. In jeder Stadt verdient er so viel, daß sich sein Geld jeweils 
verdoppelt; in jeder Stadt gibt er aber auch 12 Denare aus. Nach der Heim¬ 
kehr hat er nichts mehr. Lösung: 10 1/2 

28) Erhält ein Mann vom zweiten 7 Denare, so hat er fünfmal so viel wie der 
zweite; erhält der zweite Mann vom ersten fünf Denare, so hat er siebenmal so 
viel wie der erste. Wieviel hatte jeder ursprünglich? 

Lösung: (12 7/17, 46 2/17) 

29) Ein König sandte 30 Mann in seinen Obstgarten um Bäume zu pflanzen. 
Wenn sie 1000 Bäume in 9 Tagen pflanzen können, wie lange würden 36 
Mann für 4400 Bäume benötigen? Lösung: 33 

30) Ein Mann hinterläßt seinem ältesten Sohn 1 Bezant und 1/7 des Verblei¬ 
benden. Von dem Rest gab er dem Zweitältesten 2 Bezant und wieder 1/7 des 
verbleibenden Rests. Diese Teilung setzte er fort, er gab dem nächstjüngeren 
1 Bezant mehr als dem vorhergehenden und je 1/7 des verbleibenden Rests. 
Wieviele Söhne waren es, wenn der letzte Sohn den Rest erhält, aber alle Söhne 
gleichviel? Lösung: 6 
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Algorismus Ratisbonensis (1450) 

31) Item ain fraw hat veigen vnd hat auch kinder vnd sy gibt iglichen kind 12 
veigen, so pleibt ir 37 feigen. Nu nympt sy dij veigen widerumb von den 
kindern vund gibt ander wais hyn vnd gibt iglichen kind 15 feigen, so zw rint 
(fehlen) ir 44 veigen. Nu frag ich, wije vil sind der feigen vnd der kinder gewe¬ 
sen? Lösung: (361,27) 

32) Item: diuidatur 8 in tales 2 partes, ut veniant in diuisione 3. (Teile 8 so 
in zwei Summanden, daß ihr Quotient 4 ist). Lösung: (6 2/5, 1 3/5) 

33) Si tune vixisses quantum vixisti et iterum tantum et dimidium et dimi- 
dium dimidij tanti, 100 annos compleuisses. (Wenn du noch solange leben 
würdest, wie du schon gelebt hast und nocheinmal so viel, und noch halb mal 
soviel und ein viertelmal soviel, so würdest du 100 Jahre alt werden). 

Lösung: 26 2/3 

34) Item ein thurnn (Turm) ist 1/4 im ertrich (Erdreich) vnd 10 schuch im 
wasser vnd 3/5 ym luft. Nu frag ich, wye lang der thurnn seij vnd uil schuch 
im wasser sey vnd vieuil im luft. Lösung: 66 2/3 


Bamberger Rechenbuch (1483) 

35) Item ist ein Faß, das hat drei Zapfen. Wenn man den ersten zieht, geht es 
aus in 2 Tagen. Mit dem anderen Zapfen geht es aus in 3 Tagen. Mit dem drit¬ 
ten Zapfen geht es aus in 4 Tagen. Und wenn man sie alle drei zieht, wielange 
muß es ausrinnen? Lösung: 12/13 

36) Es lauf ein Has gegen ein Holz und ein Hund lauf ihm hinten nach. Und 
wenn der Has 12 Sprünge tut, so tut der Hund 15, und der Has ist vor dem 
Hund 100 Schritt. Nun ist die Frag, wann der Hund den Hasen erläuft, in wie¬ 
viel Schritten? Lösung: 33 1/3 

37) Es sind zwei Gesellen, die gehen nach Rom. Einer geht alle Tage 6 Meilen. 
Der ander geht am ersten Tag 1 Meile, an dem anderen Tag 2, etc. Um alle 
Tag eine Meile mehr denn vorher. Nun willst du wissen, in wieviel Tagen einer 
soviel gegangen ist als der andere? Lösung: 11 

38) Item ein Turm gebaut nach solchen Sitten: 1/4 des Turmes ist im Erd¬ 
reich, 1/5 im Wasser und 100 Schuh in der Luft. Nun frag ich, wieviel Schuh 
sind im Wasser und Erdreich und wieviel Schuh sind an dem ganzen Turm? 
Lösung: 1819/11 
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"Titelblatt von Adam Rieses Rechenbuch (2. Aufl. 1529). Es zeigt einen 
Wettbewerb zwischen einem Zahlen- und einem Linienrechner'' 


Rechenbuch des Adam Riese (1526) 

39) Item ein kauffmann zeucht (verleiht) hinweg mit gelt/gewinnt ein drittheil 
seines Hauptguts (Kapital)/vnd 4 fl. mehr/legt an hauptgut vnd gewinn/ge¬ 
winnt den vierdten theil /bringt zusammen 40. Lösung: 21 

40) Item einer spricht zum andern/wann ich noch so vil/ein drittheil/vnd ein 
vierdtheil so viel hett /so wer meines gelts über 100 fl. so viel als jetzt darun¬ 
der. Lösung: 55 35/43 

41) Item ein fuhrmann fehrt von Leiptzig ghen Nürenberg in 6 tagen/vnd ein 
ander fuhrman fehrt desselbigen tags aus von Nürenberg /kompt in 8 tagen 
ghen Leiptzig / in wie vile tagen körnen sie zusamen? Lösung: 3 3/7 

42) Item 21 Personen / Männer und Frawen / haben vertrunken 81 d / ein 
Mann sol geben 5 d vnd eine Fraw 3 d. Nu frag ich wie viel jeglicher in Sonder¬ 
heit gewesen sind? Lösung: (9, 12) 
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43) Item einer hat zween silberne Becher / vnd ein oberlied (Deckel) / so das- 
selbig auff edn ersten gesetzt wirdt / heit er viermal des andern gewicht. Wirdt 
es aber auff den andern gesetzt / so ist er dreymal schwerer dann der erste / 
vnd das oberlied wigt 16 loth / wieviel wigt ein jeglicher Becher in Sonderheit? 
Lösung: (7 3/11, 5 9/11) 


Algebra von Christopher Clavius (1608) 

44) Um seinen Sohn zum Rechnen zu ermuntern, zahlt ein Vater seinem 
Sohn 8 Pfennige für die Lösung einer Aufgabe, wenn sie richtig gelöst wurde 
und zieht ihm 5 Pfennige ab, wenn sie falsch war. Nach 26 Aufgaben schuldet 
weder der Vater dem Sohn etwas noch umgekehrt. Wieviele Aufgaben hat der 
Sohn richtig gelöst? Lösung: 10 

45) Wenn ich jedem Bettler vor der Tür 7 Pfennige geben würde, würden mir 
24 Pfennige verbleiben. 32 Pfennige würden fehlen um jedem 9 Pfennige ge¬ 
ben zu können. Wieviele Bettler sind es und wieviel Geld habe ich? 

Lösung: (28, 220) 

46) Einem Diener werden 100 Denare und ein Mantel als Jahresgehalt ver¬ 
sprochen. Nach 7 Monaten verläßt der Diener das Haus und erhält 20 Denare 
und den Mantel. Wieviel ist der Mantel wert? Lösung: 92 


Vollständige Anleitung zur Algebra von Euler (1770) 

47) Ein Vater hinterläßt seinen 3 Söhnen ein Vermögen von 1600 Reichs¬ 
talern. Nach seinem Testament soll der älteste Sohn 200 Rtl. mehr haben als 
der zweite, der zweite aber 100 Rtl. mehr als der dritte; wieviel bekommt 
jeder? Lösung: (700, 500, 400) 

48) Suche ein Zahl von der Beschaffenheit, daß, wenn man sie mit 5 multipli¬ 
ziert, das Produkt so viel unter 40 liegt, wie die Zahl selbst unter 12. 

Lösung: 7 

49) Zwey Bäuerinnen haben zusammen 100 Eyer, die erste spricht: wann ich 
die meinigen zu je 8 überzähle, so bleiben 7 übrig; die andere spricht: wann 
ich die meinigen zu je 10 überzähle, so bleiben mir auch 7 übrig; wieviel hat 
jede Eyer gehabt? Lösung: (63, 37) (23, 77) 

50) Eine Gesellschaft von Männern und Weibern sind in einem Wirtshaus; ein 
Mann verzehrt 25 Cop., ein Weib aber 16 Cop. und es findet sich, daß die Wei¬ 
ber insgesamt einen Cop. mehr verzehrt haben als die Männer; wieviel Männer 
und Weiber sind es gewesen? Lösung: (7 + 16 t, 11 + 25 t) t ganzzahlig 
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ANHANG B 


ZITATE 

Wenn jemand die Mathematiker fragt: Ihr Wunderlichen — von was für Zahlen 
sprecht Ihr eigentlich — Was werden sie wohl antworten? daß sie von Zahlen 
reden, die man nur denken kann, unmöglich aber auf irgendeine andere Art 
handhaben. 

Plato 

Geometrie ist die Erkenntnis von Gegenständen ewigen Seines. 

Plato 

Wer die Geometrie begreift, vermag alles auf der Welt zu verstehen. 

Galilei 

Ceux qui ne sont pas mathematiciens sont portes a considerer les mathema- 
tiques comme une Science inhumaine. 

E. Borei 

Gott existiert, wie die Mathematik widerspruchsfrei ist, der Teufel, weil wir es 
nicht beweisen können. 

A. Weil 


Die Zahl ist das Wesen aller Dinge. 

Pythagoras 

Ich hörte mich anklagen, als sei ich ein Feind der Mathematik überhaupt, die 
doch niemand höher schätzen kann als ich, da sie gerade das leistet, was mir 
zu bewirken völlig versagt blieb. 

Goethe 

Die Mathematiker sind eine Art Franzosen: Redet man mit ihnen, so überset¬ 
zen sie es in ihre Sprache, und dann ist es etwas ganz anderes. 

Goethe 

Ein Mathematiker, der nicht irgendwie ein Dichter ist, wird nie ein vollkom¬ 
mener Mathematiker sein. 

Weierstraß 

Die Wissenschaft der reinen Mathematik in ihrer modernen Entwicklung darf 
von sich behaupten, die ureigenste Schöpfung des menschlichen Geistes zu 
sein. 

Whitehead 

Es gibt jedoch noch einen anderen Grund für die hohe Wertschätzung der 
Mathematik: sie allein bietet den exakten Naturwissenschaften ein gewisses 
Maß an Sicherheit, das ohne Mathematik nicht denkbar wäre. 

Einstein 
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Zur Mathematik führt kein Königsweg. 

Menachmos 

Es gibt nichts Unpopuläreres als die moderne Mathematik, und auch darin 
steckt ein Stück Symbolik der unendlichen Ferne, der Distanz. Alle großen 
Werke des Abendlandes von Dante bis Parsifal sind unpopulär. 

Spengler 

Nicht etwa, daß bei größerer Verbreitung des Einblicks in die Methode der 
Mathematik notwendigerweise viel mehr Kluges gesagt werden würde als heu¬ 
te, aber es würde sicher viel weniger Unkluges gesagt. 

K. Menger 


Die ganzen Zahlen hat der liebe Gott gemacht — alles andere ist Menschen¬ 


werk. 


Kronecker 


Die Zahlen sind freie Schöpfungen des menschlichen Geistes, sie dienen als 
Mittel, um die Verschiedenheit der Dinge leichter und schärfer zu fassen. 

Dedekind 


Der Begriff der Mathematik ist der Begriff der Wissenschaft überhaupt. Alle 
Wissenschaften sollen daher Mathematik werden. 

Novalis 

Wer ein mathematisches Buch nicht mit Andacht ergreift und es wie Gottes 
Wort liest, versteht es nicht. 

Novalis 


Die Mathematik ist wie die Gottseligkeit, zu allen Dingen nütze, aber wie die¬ 
se nicht jedermanns Sache. 

J. Kraus 


Ich glaube . ., daß es, im strengsten Verstand, für den Menschen nur eine ein¬ 
zige Wissenschaft gibt, und diese ist die reine Mathematik. 

Lichtenberg 


Die sogenannten Mathematiker von Profession haben sich, auf djg Unmündig¬ 
keit der übrigen Menschen gestützt, einen Kredit von Tiefsinn erworben, der 
viel Ähnlichkeit mit dem vom Heiligen hat, den die Theologen für sich haben. 

Lichtenberg 


Die Mathematik ist eine Art Spielzeug, welches die Natur uns zuwarf zum 
Tröste und zur Unterhaltung in der Finsternis. 

d'Alembert 


So kann also die Mathematik definiert werden als diejenige Wissenschaft, in 
der wir niemalsdas kennen, worüber wir sprechen, und niemals wissen, ob das, 
was wir sagen, wahr ist. 

B. Russell 
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Stets sind sie eilig, nur zu messen und zu rechnen, halten es für die Hauptsa¬ 
che, und le calcul! le calcul! ist ihr Feldgeschrei. Aber ich sage: ou le calcul 
commence, l'intelligence des phenomenes cesse: während Einer bloße Zahlen 
und Zeichen im Kopf hat, kann er nicht dem Kausalzusammenhang auf die 
Spur kommen. 

Schopenhauer 

Das Wesen der Mathematik liegt gerade in ihrer Freiheit. 

G. Cantor 

Die Geometrie ist einzig und ewig, ein Widerschein aus dem Geiste Gottes, 
daß die Menschen an ihr teilhaben, ist eine Ursache dafür, daß der Mensch 
Ebenbild Gottes ist. 

Kepler 

Die Mathematik ist eine gar herrliche Wissenschaft, aber die Mathematiker 
taugen oft den Henker nicht. 

Lichtenberg 

Wer die Sicherheit der Mathematik verachtet, stürzt sich in das Chaos der Ge¬ 
danken. 

da Vinci 

Befaßt Euch zuerst mit Mathematik, denn die Mathematik ist die Zahl, und 
die Zahl ist Gott. 

J.H. Fahre 

Die Natur ist mathematisch, die Meisterwerke der Kunst sind im Einklang mit 
der Natur; sie drücken die Naturgesetze aus und bedienen sich ihrer. Folglich 
ist das Kunstwerk mathematisch, und der Wissenschaftler kann das strengste 
Urteil anlegen und unfehlbare Formeln anwenden. 

Le Corbusier 
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Für den Anwender des Schneider 
CPC bietet dieses Buch eine Reihe 
von Programmen aus den Bereichen 
Finanzmathematik, Unternehmens¬ 
forschung und Betriebswirtschaft, 
Sie finden z. B. Programme für 
Zins-, Rendite-, Renten-, Tilgungs-, 
Optimierungs-, Abschreibungs- 
Berechnung. 
lnVorb,2.Q.'85.216Seiten. 

Kart. DM 44.-/Fr. 44- 
ISBN 3-88322-153-8 


Das Buch führt den Leser mit vielen 
Beispielen aus der Praxis in das 
Datenbanksystem dBase II ein. Es 
sind keine Vorkenntnisse von dBase 
erforderlich. Das Werk ist ähnlich 
einer (PU)Programmierten Unter¬ 
weisung aufgebaut und ist daher 
sehr gut verständlich. 

1985. Ca. 300 Seiten. 

Kart. DM 48,-/Fr. 48.- 
ISBN 3-88322 154 6 



Dieses Buch wendet sich an Super- 
Calc-Anwender. Anfänger führt es 
schnell zur sicheren Handhabung, 
Fortgeschrittene finden ausgereifte 
Anwendungen. Beispiele mit aus¬ 
führlicher Beschreibung machen 
das Buch auch für VisiCalc-Anwen- 
der interessant. Die Unterschiede 
zwischen VisiCalc und SuperCalc 
sind in einem Kapitel dargestellt, 
ln Vorb. 2. Q. '85. Ca. 200 Seiten. 

Kart. DM 48,-/Fr. 48- 
ISBN 3-88322-155-4 



Teil 1 gibt eine Einführung in die 
Grundlagen; Sprach-Interpreter und 
Compiler werden erläutert. Teil 2 ent¬ 
hält eine umfassende Beschreibung 
aller Anweisungen, Befehle und 
Funktionen, abgerundet mit zahlrei¬ 
chen Beispielen. Teil 3 beschreibt die 
Anwendung des Compilers und sei¬ 
ne Unterschiede zum Interpreter. 
1983 306 Seiten Spiralh. DM 56,-/ 
Fr. 56- 

ISBN 3-88322-024-8 



Dieses Buch behandelt alle grafi¬ 
schen Fähigkeiten des Schneider 
CPC und führt in die Grafikprogram¬ 
mierung ein. Von der elementaren 
Programmierung bis hin zur komple¬ 
xen mehrfarbigen Abbildung im 
hochaufgelösten Grafikmodus wer¬ 
den zahlreiche Beispiele erläutert 
und entsprechend illustriert. 

ln Vorb. 2.Q, '85.288 Seiten. 

Kart. DM 48.- / Fr. 48- 
ISBN 3-88322-147-3 
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CP M für die Praxis 1 
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Das Buch ist in drei Teile gegliedert: 
Teil 1 führt in die Eigenschaften von 
Computern und CP/M im besonde¬ 
ren ein. Aufbauend auf diesen Kennt¬ 
nissen werden im 2. Teil die zentra¬ 
len CP/M-Hilfsprogramme vorge¬ 
stellt. Der 3. Teil geht-nach einer Ein¬ 
führung in die Funktionsweise des 
8080-Prozessors - auf die CP/M-Sy- 
stembesonderheiten ein. 

1982. 386 S Mit zahlt, prakt. Beispielen 
Geb. DM 48,-/Fr. 48.- 
ISBN 3-88322-004-3 



















Dieses Buch beschreibt den Lei¬ 
stungsumfang von FORTRAN-80 für 
Personalcomputer, wie es im Relea¬ 
se 3.44 von Microsoft festgelegt 
wurde. Alle Sprachelemente sind in 
detaillierter Form abgehandelt und 
in Beispielen dargestellt. Das Buch 
kann in dieser Form als Program¬ 
mierhandbuch und auch als Nach¬ 
schlagewerk genutzt werden. 

1984.344 Seiten. 

Splralh. DM 56.-/Fr.56.- 
ISBN 3-88322-079-5 


Dieses Buch ermöglicht die schnelle 
Einarbeitung in die Programmier¬ 
sprache CBasic. Alle Eigenschaften 
werden mit Beispielen erklärt. Syn- 
taxdidgramme zeigen grafisch die 
Einsatzmöglichkeiten. Die Funktio¬ 
nen des Compilers, des Linkers und 
des Bibliotheksprogrammes LIB 
werden erläutert. 

1983.190 Seiten. 

Splralh. DM 56.-/Fr. 56.- 
ISBN 3-88322-034-5 



Dieses Buch behandelt die Proble¬ 
matik der Schnittstellen V 24 und 
dem lEC-Bus. Die Übertragung zu pe¬ 
ripheren Geräten (Drucker, Plotter, 
Terminal u. a.) ist gut verständlich 
aufgezeigt. Hardwareaufbau und 
Beispielanschlüsse machen das 
Buch zu einem guten Nachschlage¬ 
werk. 

In Vorb. I.Q.'85.200 Seiten. 

Geb. Ca. DM 48,-/ ca. Fr. 48- 
ISBN 3-88322-094-9 



Wer hat nicht bereits verzweifelt ver¬ 
sucht, das »Computerchineseng- 
lisch< zu verstehen? Hier hilft das 
Wörterbuch der Computerei mit sei¬ 
nen über tausend Begriffen. Außer¬ 
dem sind die wichtigsten Begriffe er¬ 
klärt. Ein handliches Nachschlage¬ 
werk für jeden, der sich mit Compute¬ 
rei beschäftigt. 

1983. 2., erw.Aufl.144 Seiten. 

Kart. DM 32,-/Fr. 32- 
ISBN 3-88322-026-4 



Die Programmiersprache >C< ist be¬ 
sonders zur Erstellung schneller, ma¬ 
schinennaher Programme geeignet, 
weist jedoch gegenüber der Assem¬ 
blerprogrammierung wesentliche 
Vorteile auf, wie z. B. lokale und glo¬ 
bale Variable, Prozeduren, Funktio¬ 
nen usw. Diese Einführung setzt 
beim Leser keine umfassenden 
Kenntnisse voraus. 



Dieses Buch führt Sie zur sicheren 
Handhabung der Tabellenkalkula¬ 
tion. Zahlreiche Anwendungen mit 
ausführlicher Beschreibung schlie¬ 
ßen sich an. Unter anderem: Werbe¬ 
planung, Reisekosten, Kapazitäts¬ 
auslastung, Baufinanzierung, Reak¬ 
tionszeiten- Auch für die englische 
Version nutzbar. 

ln Vorb. Juni 1984. Ca. 250 Seiten. Spiralh. 
Ca. DM 56,-/ca. Fr. 56.- 
ISBN 3-88322-074-4 



1984.256 Seiten. 

Geb. DM 56,-/Fr. 56.- 
ISBN 3-88322-041-8 


















D. Herrmann, G. Schnellhardt 


Schneider CPC Band 1: 
Mathematik 

Dieses Buch enthält 39 mathematische Programme aus den Berei¬ 
chen: — Mehr-Register-Arithmetik — Zahlentheorie — Kombinato¬ 
rik — Algebra — Geometrie — Numerische Mathematik. 

Die hier neu vorgelegte Langzahl-Arithmetik gestattet die Grundre¬ 
chenarten für Zahlen bis 255 Stellen. Zahlreiche Anwendungen fin 
den auch die angegebenen kombinatorischen Prozeduren. 

Sie finden auch das Rechnen mit Polynomen, Matrizen und kom¬ 
plexen Zahlen: Neben dem komplexen Hornerschema werden ins¬ 
besondere Algorithmen zur Polynomdivision und Matrizeninversion 
gegeben. 

Anwendungsbezogen sind die Programme der Numerischen Mathe¬ 
matik. Hier werden Verfahren zur Lösung von linearen und nichtli¬ 
nearen Gleichungen gegeben, zum Eigenwertproblem von Matrizen 
und zur numerischen Integration und Differentiation. 


ISBN 3-88322-152-X 
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